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PREFACE 


The present volume has evolved out of the lectures delivered 
by me to the Post-Graduate students of Pure Mathematics in 
the University of Calcutta. It aims at giving a fair, comprehensive 
and connected account of the subject of Finite Differences. 


‘The subject-matter has been presented in а form which will 
be easily understood, and the different problems have been 
domonstrated in as simple a manner as the topic under discussion 
would allow, In many places new and original proofs have been 
given. ‘The needs for beginners in the subject have been kept in 
view во that the book will be found suitable for beginners as also 
for advanced students and researchers. 

A comprehensive account of symbolic operators, of finite and 
approximate summations, and all the various propositions of 
finite differences bave been given. Interpolation bas been 
discussed at greater length. Difference equation bas not been 
treated in thin volume, as the same is intended to find place in 
а separate volume. 

Every branch of scientific enquiry has got a practical side 
but unfortunately in most cases little or no stress is given 
in that direction so that the learner knows the theory but cannot 
make practical applications of bis knowledge. The theories of 
Finite Differences admit of wider applications. They are 
required by scholars in many different fields—Astronomy, 
Actuary, Engineering, Statistics, Physics, Meteorology.......08 
well us Pure Matbomatics. An attempt bas been made here 

X to combine theory with practice by bringing the ignored practical 
side into some prominence. The subject-matter covers an 
extensive feld snd has been amply illustrated by numerical 
‘examples and applied to data taken from observations—actuarial, 
‘astronomical......s0 that the worker may learn how to apply his 
‘theories to diverse problems. 























x PREFACE 


The present volume has been designed to meet the require- 
ments of students of Pure Mathematics as also of students 
Preparing for actuarial, statistical and other natur сө 
examinations. Я 

In the preparation of this work I Һай to consult the existing 
treatises and various memoirs and research papers on the subject 
and my thanks are due to their respective authors. My own 
researches also have been incorporated in proper places in the 
body of the book. ^ 

Any suggestion for improvement will be thankfully received. 

In conclusion, I wish to express my thanks to the officers and 
the staff of the Calcutta University Press for the pains they ha 
taken in printing off the book. 








UNIVERSITY or CALGUTTA: 
} P. N. MITRA. 


March, 1939. 





CHAPTER I 
DIFFERENCE OPERATORS 


1.1. Introduction. Calculus of Finite Differences, a» its 
name implies, deals with finite increments of the variable or 
argument aud the corresponding increments of the function 
dependent on the variable. 

Let f(z) be a function of the variable x, and let zo, жү, ту, ... 
be the successive values of x, such that the intervals 


بو “و2 .20—21 .21—20 


are finite. These intervals may or may not be equal. When 
these are equal, to Ar say, the successive values of x can be 
put as 








z, reAr r4DAz, 2434s, .. 








and the corresponding values of the fonction as 


Қа), fir Az), [(e+2d2), f(r48Az), .. 








"The difference of f(x) from its successive value, 
a) = + Az) – f(z) 
called the first difference of f(z). We shall be mainly con- 
med with such differences. In most cases the interval Ar 
will be taken to be unity. 
When the values of x increase by unequal finite quantities, 
the intervals of the variable x are also unequal. We have such 
cases in Divided Differences, Reciprocal Differences, .. 
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If the intervals of z become infinitesimally small, we -come 
to Differential Calculus. Thus from (1) we have 


flz+Ax)—-f(2) E 
ar 


which in the limit as Ar—>0 gives the differential coeficient 
ot f(z). Thus 
Pig) = Lim 2ھ‎ fix). 
ere A 
. 1.3. Difference Notation. Let Ax be the increment of the 
variable z, then the corresponding value of the function is 
fle+Az). The difference 
Har Az)- f) 
in denoted by Af(z) and is called the Firat Difference of f(r) 
with respect to tho increment Ar. Here we shall be concerned 
with finite constant increment. 
Hence denoting A by w, we have 
А+) e ]( + w) — f(x) 
as the first difference of. f(z). Similarly 
Af) m Afr ы) Afla) 
` fs called the Second Difference of f(x), and 
AM (2) = Дауа + w) – АЗ) 








m im the Third Difference of jis). Proceeding in this manner we 


сап form the fourth, fifth,.,...nth differences of f(z). 
Consider for example f(z)sz*, then 
As) = + Bru? eut, 
ArI = ru + Gu, 
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Argument Function д a a at 
a fi) - 
2 Аа) 
ate Kat) Atla) 

Afiat w) Atay 

4+9 fla+2u) Af =) A*f(o) 
ма +) Afla же) 

ae — fla+3w) AS fa +2) 
Afía Se) 


44e — f(a4u) 


and similarly for differences of orders higher than the fourth, 
‘Tho first entry f(a) is called the leading term, /(a + u), fla+2u), 
aro called the successive values of the function, and the 
differences Afla), A*fía, A*f(a), „are called the leading 
ifferences, Each difference in the Table is obtained by sub- 
tracting the upper entry immediately to the left from the next 
lower entry in the same vertical column. It may be noted that 
to form a fourth difference, five consecutive entries are required; 
and generally to form a difference of the nth order, (n1) 
consecutive entries are necessary. 


Етатріс. Difference table of f(z)— x, w=1. 











Li tin) А AR ar ae 
1 1 
7 
2 8 12 
19 в 
a 27 18 0 
n эт 6 
4 (4 24 ° 
„ о [] 
5 125 30 
өт 





o 
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4.3. Operator E. This operator denotes the opiration of 
increasing the argument by w. Thus 
Еу) = fia +»). z 
‘The operation can be repeated any number of times. Thus 
E? f(z)  Ef(x + w) = f(x +2), 
E? f(z) = Effz + w) = f(x + Bw), 








E*^f(z) = E" f(z + o) = für n), 
where" is an integer. We bave by $ 1-2 
oA Alla) = fle + d — f(z) = Eft) — fix) 
4 = (Е – уба). 
` Awo Низы) а) + Afta). 


‘Tt is thus evident that the operators E and A are connected 
by the relations 


a &еЕ-1,‚ or Em1«A. 
Again by Taylor's Theorem, we bave 
Ера) = fis) 


| 
| 
i 


ы? d 


LUE + r + а 








DIFFERENCH OPERATORS 5 
1.4: Properties of the Operators. ‘The operators A, E and D 
follow tbe fundamental laws of Algebra. We bave 
AG) + gr) sues s 
= für + w) + gru) nes = fin) — gir) — 
= 2/02) + Ag) + .... 
Thus A satixfes the Distributive Law, 








Again. A"A"J()=(AA...m timps)(AA...m times) fix) 
=(AS.....m + п times) f(z) 
=۵". 

‘Thus A satisfies the Index Law. 

Finally Дауа) = аре + ы) — afi) aafe), 
where а is constant, ‘Thin result is also true if a is replaced by 
a periodic function a (х) of period w. Thus 

(у) = ala + fir + e) ау 
= (е) + w) а) +) 
=н). 

Let A, operate оп z and A, on y, and let and ы? be the 
respective increments of z and y, then the partial differences of 
Jz, y) are given by 

Ale = fir ө, y) — fir, De 
Auf, y) für, y +) — für, у). 
Agen — A,A, fir, y) m Ay +, у) v) 
fete, yu) far, fix, у + ы) + für, у), 
. А,&,/(х, у) = Au (fin, у n) — Кг, У} 
m fin, y +e) Јн, у) fer, y) e fn, у). 

Hence ۵.۵, = AVAL 

Thus the operator A satisfies the Commutative Law with 
regard to constants and periodie functions of the seme period 
and to the operators not affecting the argument on which А 
operates, 


CALCULUS оғ VIXITR DIFVEARXCHA 
I moy be shown in like manner that E and D alid satiety 


e the fundamenta! laws df Algebra subject to the same condi- 
le’ tious, Thus'im many respects the oporstors A, E and D bebave 
| 





“Hike Algebraic symbols and may be combined ая such. 
1.5. Differences of Elementary Functions. 
0) Afü)gun = fle + ugi + e) — Naga) " 





id mjia + og + o) = flzlgle + wj [бедра + w) — fang) 
| == + Afi) + Аида). 
Putting «1, we bave 
Ану = gor + ПАЙ) + fingit. 


Agan Aflergieh= (fori + Afta Hola) + Agi} — feign 
farte) + уе) + )اود )د‎ 


К moon а-а 








ч N Ep + Ag) 

| Eur . 

FA sinu, ена (v, Au, ) uin wy 

jt EIL eos (u, میگ‎ ( 

| 7 ما د“‎ Spe rious єр), 

| د‎ sonno mri ма Ste عر + ,)مہ‎ ). ^ 


" A ens nhan ч, tam eum "Дел. 3279 








«dia 
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Лы, maz +b, and 1, we bave 


Asin عم‎ 2 sin р دا‎ (gr ene SS" } 
Similarly 29% sin (az b) m (2 sin ф” sin (ев eati), 
* к 7 ate 
and generally A* sin (ar eje (2 win $) (овен). 
2 me et 
Similarly Д cos (5+) =2 sin pon (pr ttt) 
* = 3 et 
and а" oos (ax +b)= d ма g)" con (orden tz). 


sina 
tan ру ome LTTE 


A tan"? (are bm ten"! aad 





0 Astmete (a — lja*. 
Hence Aat (а= = та", 
and generally. 2" m (em =1) "a". 





Agin Aunt e (a — Lan, 
Ata = (am —l)*atn. 
UI. Abg ule da. — log us 


mese rr) Rr ry 








СЕКАДЕ Monee d) m log d 


pon 











8 CALCULUS OF FINITE DIFFERENCES 
1.6. Differences of a Polynomial. Consider the polynomial 
f(x) e ао” каут" casa 79 + ......+0n 

Now Аа) Йе) — f) : 
SAo{(e+ 0)" —2"} * a (( + w)*71 DII 








+ 





ва. 





RES [nr ненац 





of degree n — 1 in x and therefore the first 


This in a polynomi 
ial of degree 


difference of » polynomial represents another polyno 
less by one unit, Repeating this process, we bave 
A? f(z) s ag n(n —1)2*7 32 +... 
a polynomial of degree » —2 and soon. Itr <n, we obtain on 
differencing т times 
Ar fl2)=agn(n—1)(n—2).. 











тит" 














Finally A*f(z) вом" — 1)n —9).......2.1,*. 
Putting w=1, we have A*f(z) = вот" —1)(n —3) 
In particular A^ =1,2.3. 





"Thus the nth difference of в polynomial of degree т № 
constant and the (m+1)th and bigber differences are all zeroes, 


1.7, Differences of Factorial Expressions. Producte of the 


Na)ftz— f(z — 2o)... ..f(z ты + 9) 
and fixie + f(z + 24)... fer ть) 
where m is а positive integer, nre called Factorial expressions, 
the first being а descending factorial and the second an ascending 


factorial expression. 
‘The simplest ond of special interest are the following 


factorial 








arr — Ца 2). 





1 

EEFE +i 

If m=0, the value of these expressions is unity, 
Mate,” 
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Now Ar™= (e+ Пет 





(xm +2) 


—#(z—1)fz—2). 






maim eDim). 


Again Ar em A i7 D m (m —1) 2-2) = m gin 


and generally Ati = m(m —1)......(m —n +I- юата), 





Tt may be noted that the sum of prefixes of m and z is always 
equal tn m and thus the number of factors in m and x is 
nhwnys m. 

Again Arm 


1 
vm) Жеан] 





sey 
(2+ 100 +2) 


SHS iy 3. m) 
в-он), 
Hence by repeated operations af A we hate 


Atre —m(- m1). m т A Dt? 





= 


= (1 mim + Dm +2)......(т +n = Пя. 





Here the number of factors in the denominator exceeds the 
number of factors in the numerator by m. 


The two oases can symbolically be put in the single form 
A tattoo ہے‎ eaim) 


where m may be а positive or negative integer. 
For the general form of fuctorials we have 
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И u,=or+b, we ha 
(а= bj. 








m (az + b)i, 





A (ax € b) = a? m(m — (о + Б), 





И т=н, Ата +b) m atm 1 

and if m <n, A*(az + b)" 0. 
Again Alor +Ьу-® = ama 4 b) ==, 
Araz + bj = ( — 1) m(m + 1)......(m n — 100+ 6) 





For the interval w we bave 


am )چ‎ u)(r-me)es(r mete) а 





ЕВЕ Е 
GA hills SEF oa Do), 


Hence Ar“ = muzin, 
Atm mim —1 m —n + 1ана ce, 
Atte ایر‎ 
AAEM = (=m) (т 1). (т от + отне сты), 


1.8. Difference Quotient A. We introduce here Nórlund'& 
‘operator A which is defined by the relation 


lteter‏ ررد 





e+ meu) 








We call Ai(x) the first difference quotient of fiz). The 


Operation can be repeated. Thus 
f(x w)— Afis) 
aja] - علد‎ 





аза) 





+ Qu) — 





DIFFERENCE OPERATORS n 
Апа generally for the nth difference quotient 
ania [279]. 
We obtain from this the relation. 
Alera) m eA" firj + Afin) 


The operators А and A а 





connected by the relation 
At mut. AT. 
When ө = 1, both the operators become the same. 
When «—>0, we have 
Lim A /( = Dji) 


=o 
and Lim A*firj-D"f(r), where р-#. 


1.9. Gregory-Newton Formula, Expression of u,... in 
terms of u, and its successive differenc 


We have ty. qu = Ем, = (1+ Д)" 








1+ + Das, 10-2 aay. уд, Ju 





Hence 


nin 








а) маме =u, emu, + Atu, +... AT. 





2T 


Ezample. Express x* in descending and ascending factorial. 
We have by Gregory-Newton formula 





=ч, ды, + ج‎ зи, +...+ Atu, 
Put u,=2", then 
(ee ر‎ = ر٣‎ адут TEL Day 
Batting y=0, we have 
ато" + до, + 670дан + 
giving z* in descending factorials. 





+ary". 





+270" 
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To express x" in ascending faetorials put —2 for £ and ч 
let w,-z**. Then ` 


(yat tet алута + Dae x 





ziz +1)(2+2) 4s 
as d ammi Any 
Putting y=0, we have 


(a)i = —>Ао* 





2+1) asosi 25 +1148) Aspect 
+ T3 atomi ae 0 +... 


orate (—1)"{ Aone — t$ aso, EHNE Daora 





Expression for A*u, in terms of u, and ite succeesi 





-(E-D*u, 





- [met -Up 0-0)" Ju 


a ninar + PET, rr 


When =1, we have А 

Ашчы Pann +" 1), аса. + CHI. * 
ейди, жшше i e hare 

A дле (кта en те ас: аът 2) (тат, 


=u, oe + (=1) thay 
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1.14. Theorem. If F(4 be a ratio 
with coefficient 


il integral function of t 
which are independent-of t, then will 








n F(Eja*$íz) = a*F(a* E)o(z), 
Let Fit) Ҳа, 
then F(Eja* (x) = Xa, tela) 


= Хана" eir ne) 











In particular when «1, we have 





@) F(R)a*o(z) Р(аЕ)Ф(®). 
From this we easily have 
NS Ра" P (aEja-*olz). 


1.12. Expansion of a Polynomial in Factorals. Let а(х) 
be а polynomial of degree. m. To express it in factoriale assume 


(1) pla) ag + apa + азии — ш) + азы) — 2e) +. 
| 7 чать) — 2ш)...(хк— me =) 


which is а polynomial of degree m with m1 constante, 
| Differencing it we have 


Agr) ma ү» + 2а or + Ba ушл(ж — e) +... + maur س ا‎ 








Apr) a, + 2a ә + Bag 2 +... + magn Te, 





14 CALCULUS OF FINITE DIFFERENCES 
Now putting z—0 in the above equations we obtain «f£ 
POs ApiO)=ay, A%(0)=2ag, ... 


Ie = 2790) =г1 а, 





000) = т! a... 


Substituting these values of do, ву, ау, ... in (1) we have 








@ $i) 90) +2910) + 7; 2.4960) + E 


+ [ro = 





аео). 


| When w—>0, we get Maclourin’s Theorem in Differential 
j Caloulus 


4 (3) ole) = 900) +2090) + De) +... 2 = юч). 
When w=1, we bave 


(8) 0) 60) адн) у два). 





Р 
эө. 





_ $43. Secondary Form of Мас!ашп'в Theorem. То expand 


g(a) in the form 
08.28 а 


Ке 


orgs, 


eI ts 





727 establish this result we shall first prove the | 
ond фа) be two rational 




















DIFFKKRENCE OPERATORS as 





Then s мо seas ( 4 уле 
Еу 2 
Now (E ymo itm 


=0, if men and 2=0, 





=n!, it men, 


Е 


‘Thus the corresponding terms of the two series are equal 
and hence 


а 9 is = д o. 


Now Maolaurin’s theorem (3) of $ 1.12 can be put as 





a 
@ Одо 060) 57+... 


Е И (1) we get the secondary 
form of Maclaurin 


e c 





о ое. 
- ЕО, О 


а: 








a6 CALCULUS OF YiMiTE DIFVEMEXCES 
4.18. Herschel's Theorem. If fiz) be a rational integral 
function of x, then 
| е") e fio: egy ++ fur 97 





Asume ег) ао ае! base?! +age +... 
| This cun be pub 
a eint DI + age + 
M Hence f(e*) = ((ao +0, E e eg E? +... 
TB haso, 
= (Bee! 
зву, 
Now ГЕЛ = f().1. ык 


Ie =f) + (E)0.t + fc) 95 








where x=0. 






















" 1.18. Deduction of Maclaurin's Theorem from Herschel's. 
Theorem, and conversely. 


We have 
(t) ellog е“) ee(log E)c?**, by Herschel's theorem 
=ollog eP)e9*! ee(D)e9*, 





which is the secondary form of Maclaurin’s theorem. 
‘Tu establish the converse cuse, put /(е')= ө). 


ties gena ($ )-Ц ив. 
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Ezample 1. Develop in ascending powers of t 








Q il. 149 nh. 
@ Here. i= тк = dy ey Herschel's theorem 
"ата «(1+2 ) 
a[- 










orat 
T, 





о”. 





Hence f(e!)== $ [-2 +p r+ (oe St } 


20°, A?*0 


te. т = Ba { 0° А+ oz" -. +(e IE: 


zi 


ün - 








+ )-1( 25- Je 


y= 240" } ee 


safo- + 537 + ао Mu 





and calculate tbe first four terms of the expression. » 
„19-0. М.А-4 M.Sc. 1938] 
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We have 


Ке) me! +A, 60-1, by Herschel's theorem 


2e 0t 








RETE A 0. 
efe 1+ sotte © 


А. 





[+20 tee 
giving the result of the first portion. e 


° Now Кет) ее (1 0. t ое. 





uu 9 
ore os 


E. ДП 








Putting х= 0, 1, 2 and 3, we have 
e mes ots et Jet? +... 
1.46. Lagrange's Theorem. То prove that 
Se pie ie сы 
а^ +0! * 
. 
where n in а positive integer. 
We bave Au m (B-1) u, = (e "n, 


a 








| 
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w by Herschel's theorem. 





1)" = (E-1)*e9:( ater! 


AS ors 
a [ro 3t 





But A*0"=0 if n>m, aod =n! if nem. Hence 


=4 870", 
RESETS 








Replacing by D we have 








а) Au, e+ 








w+) МТИ 


147. Generalised Lagrange’s Theorem 
We have by Taylor's theorem 


маем, {Dug + D*u, + m : 


Operating by (4) on both sides (A operating on ! only) 
we have 


(at. бие + р. OE +...‏ + ھاو = .دام 
‘Thy first term. 9(S)u, =4,9(d).1=u,9(0).1.‏ 
Now putting 1=0, we have‏ 


оды, = olO)u, +9500. Че 00-07. e 





+ Putting #(2) = 2", we get Lagrange's theorem, $1.16(1). 
4.48. Expressions for Ати... 






t = Bu, Er, мис, 


© 
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where E operates оп u, alone and E! on v, alone. Then 
Ам „т, = (FE! uaz 
Hence A*u,v, = (EE! — 1)*u,v, . 


=(-1)" [orones 


n(n 
ES 
Agnin to express it in ascending differences of ve, we bave 


Aus, = (EF) -1)*u,0, m {A+ EAP use, 





- [aren -1 BAM + ae A" EAP + dee ы 





Аби», nA" a, a AD, Mpp latu, vae, +.. 


(we havo dropped the accents of A! as no longer necessary). 
Putting ч, 


*, we have 

















а" ©“ fi -1)"e, та 1)" Заде, 
+" ОН aon Atoa +... |. 
1.19. Expressions for Du, and ¢(D)u, in terms of 
differences of и,. 
We have «<P =1+A. hence D=log (1+۵). 
25 Du, log (1 + Aju, 
С, A 5 
(А+ 
C 
5 








ig —n for n we hi 


555.099", = (log 1+ A))-*u,. 


Again ug. =u, ttan, + дво. 
Operating by o( 3, ) on both aides we have 


OCG) mere e fi) ам, 


а уел) As, 
(d + uai 


Putting ¢=0 and noting that 
[eG es] e eG 


we have 


*( & Jusu ooto Ж, Jo. au, +9( 2). a^ 








‘of which the above are particular cases. 

1.20. Differences of Nothing. The differences of the various 
“positive integral powers of =, in which the increment of x is 
taken to be unity and x is put equalto zero after evaluation, 





‘They nro represented by the symhol AO". To 
x^ is differeuced m times and then x js put equal to zero. 

They are of frequent occurrence apd play an important part 
in the Calculus of Finite Differences. We establish here a 
formula connecting differenees of mothine, from which we are 


127576 


N 
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able to cunstruct а table of such differences. We shall first 
И establish the formula 


(Е Ю* = Еф(Е)0"-1. E 
Leta,E^0* and Era, E'-!0*-! be the corresponding terms 


[ of the two sides. Then, since each is equivalent to a,.7^, the 
f above identity is established. 


Again since B=1+A, we bave 
90 + A)o* = Ее + Ajo! 
or Као"  Efr(agon-t. 
Potting /Д) =", we have at once the recurrence formule 





ы в) Ate Бут дл-ЗО" 1 = m(A"-10"-1 + A051), 
Ii is to be noted that A90! =0 and A!0! —1 

IL. and ATO" =з! for т=п 

: -0 for mL. 


"The above formula (1) can also be obtained from the expansion 
| ef Az" of $1.10. We have 


(27 тт Шш а айй таш бат) 
et C Ttg. 


(Putting z=0, we get 
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We give below n table of differences of nothing: 


п до, Atos  A*o* дч AM aton 
1 2 : 

2 1 2 

з 1 в 6 

a 1 14 36 24 

5 1 30 150 240 120 

6 1 вз 540 . 1560 1800 720 
m 1 126 1806 8400 16800 15120 


To find the value of а particular difference AMO", add two 
numbers in the preceding row which are situated i 
column and in the preceding colum», and multiply by m, 


A106 = 4(240 + 150) = 1560. 


1.21. Difference Quotients of Zero. When the interval is w. 
tho expression 2"0* is called the Difference Quotient of Zero, 











Now Az" = Anr- zm ASII man Mr ы)*-1, by $1.18. 
.. Дт" = moAT- 1 EO! = maf ATO) + A-1071). 


Непсе we have the recurrence formula for difference quotients 
а) ATO" тым"! + тд" 





Putting w=1, we have A0" = m(A70*-1 + д-10*-1). 


1.22. Differences of Powers of Unity. We establish here 
the recurrence formula for caleulating the differences of powers 
of unity: 





AMP mne 1)A* 1^ mt 
We bave A*17*! = A*EO**! аде + 4)01 
sd Har ort aroni 
- - = (n+ 1A 10" + A707} 
+n{4"0* + 4^-0}, by 51.20 (1) 
= (n+ IA EO" + nA*- EOF. 
in Otte нА". q 














рб = 
|4 i Ед 
Goi сотта or rima буаз еке 
| pa 
D Example. Shew that agar” 271^! « A^ "0^ — 
ned 


x | We have Ath +A 11^ (n JA*I^« nA" INSANE 
(падт arty 

=(n+ ar ЕЛ” 

en 1)A* "M Bir"), 0, 

| = (+a. 





H رھ کروی‎ 
E n+l 7 





^ (1.23. Generating Functions. If 4í(f) can be «реле in 
. the form of а series Xu,t*, then 9(¢) is called the Generating 
Function of u, вой i» represented in the form 





= Guy. 
(0 For example, we have tmz =0. үр. 
S) Зак ет 
^4 p Mar 





e(t Xu,t 
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Again Qu, a 7 Gu, (1-1 4(4 


or СЄ -1 oto. 
GA*u, = GA(u,.,—,) 


(s, a uus Me) 
=( 4-1 ee. 


(2) дати, = ( 1-1 Yow. 


Qv. Уно (а) "wo 
=o +m( 2-1) oy + "OFI 1-1 Js 


=Gu,+nGAu,+ nin aanu, tate 


| 


Hence we have the expression for м,,„ alrendy obtained 


-G [pens onan, +. 

















in $1.9. 
Also e(t) = Хи." = Gu, 
tolt) = Хи," = би, 
(8) pN = Xu, „= Gun 
Gu, = (0 —09(0, 
" GAw, у= (1-090, 
* алла,.,=@ч,-9бОч,.у+@ч,-а 
-(0-—10t9(0. 
d XE GA*u, a= ü — ^el). 
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4.28. Theorem. If += бы, and (t) = Xo,t*, then ill 
оноо д у. ч 


and ets (1 )- ava, 
Tho coefficient of (= in (бу) 









дум +95м, + Oye) Fagg Fags + << 


The coaficient of t^ im +01) 





=... AL) Maca +аом, dius, ү Hay 


=ИЕ)ч,. 
iA ашы „= 














DIFFERENCE UPERATORS ат 
Putting w(E)-A"-(E-D*. ¥ H )«( 1-1 M 
(au: 
Putting YE=( $ )"=(1-2 )".v¥( )=а-0- 


and hence GAU,- a= (1 —0*e(0 
Example 1. Shew that 


and hence * балч 





(s Y 
маечки + О Даи, aee DO дэм, „у... 
Wo have by §1.28(1) 


Gu,. 





о-оо 






онно а зро + MOF?) 0+. 
Hence by 81.2808) we have 


Gu, m Gu, e nGAu, +" Naaru, x 


EDED awe, аж... 
whence the result. 


Example 2. Shew that 






малечу тои, + Vary, 


е + meer) {n+8r—2) vs, 


J We have by 61.2303) (7900) - Gua... 
Hence by $1.24 


e( il )«o-( i-i Jon "em Au, e 








8 P 
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Similarly "(2 —1)"9)=Ga%,.5,, amd во on. 


Here we shall have to expand Qw,.. in a series whose 
terms are of the form 


| fe( 1-1 X} ш). k-1, 9, э... 


| Pot e( + -1 у=, eas і -1 





or i-e. or у=1+ғу" where s=} 


д Ву Lagrange’s theorem, if у = u + apy), then we have 
А 

E ipfe) кенг)» уу (eren) +... 
+ 

g Taking fly=y*. yuj=u” and «=1, weave 





ees saro, ala +B (attr ang., 
я 


Substituting the values for y and 2 we have 
dem (а) КЕБУ 


eemper e(t 








P» tires (nar Dg 4 ( 
GU... - GU, «nGAU,., + ддар, а, 
ybi = ee IEA 






DIFFERENCE OPERATORS 


MIsceLtaNnous ExaurLES 
Example 1. Sum to infinity 


> 15 
аталат 1-8 дз 
дат" + LIA 





[C. U. M.A. & М.е. 1938.) 
АЕ ас мы 

The series ah да+ کا‎ arp Stare |= 

=д{а+д% уа” 


=a te 


=". 
Example 2. Shew that 


f 


п and a being positi 





ти уте de Ге) Ап", 





quantities and A operating on n alone. 
(C. U. M.A. & M.Bc. 1946, 1935.] 














CALCULUS OF FINITE DIFFERENCES 


Hence je rite" dr AT T(a)jn-* 





* -Ta)JA*S 





Example 8. И == 


4 ү° До» d 20° , d* 220° з. 45 
(a) et rt mene 


[С. U. М.А. 4 M.Sc. 1933, 1935.] 


‚ establish the formula 








аз 
тм ыт” ы de =) 








1 Xie) jo 
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Example 4. Show that 
n! E 


б-т)! mst 








Arone nator? Met YAa-07-24. 


Now 





As =e ca ener! SOLD aca н=л) Dp „у, 


Hence A0"—n0*-14. mero -a min. О-о im 


A1 10911 (m +1) (A"0* + A7*107). 


‘The last term =—_"! | "2(n—1)(n—2)... (m1) Aeg. 
in-m)i n— т}! Ыы 


Also 





>. ‘The given expression 


= В)» Но 





manot +а {тоз +" То, 





a min 1) (n 3). (m+ on} 


nom) 
= 4۵*0" + А”Ао* 
= д"0" + дт. 10" 
ПЕВ 











ExAwrLES 
1. Find tbe first differences of the following functions :— 
() tan 20. (i) cot 20, (ai) asin P. 


(wi) eot 2*8. 
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в 
(2) RY би 
(е “ (a£ 








log (2 sin 270) | (zio) очаат, 





rest : "T 
Ge) are tis шш tan шр. 


2. Prove the following relations :— 
(0 0*—A0**A*0*—...£ A*0*— 








Ало” 
A +... 





d) S 
8. Prove the following identities: 
© asor =" 1 DA noe, 


(ii) (n+ 1)A*0" - 2(A7-10" + 470}, 
A. Eatablish 
HO д)" =н + д)0". 
р 5. Shew that if m be less than r, 
{1 + log E)" = r(r—1) (672)... ir +1). 

| [C. U. M.A. & M.Sc. 1936.] 

6. Shew that 

{log (1 + A) ] "f()0* = и — 1) (22)... (e n +100400", 


7. Prove that 
{log Е} "0" —-0, 


s unless mem, when it in eqonl to mt. " 
8. Shew that 


nao" =+ n {Блу 0°. s 


[C. U. М.А. ФМ. 
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9... Prove that 


фана por 


marah 


е deduce the expression 
A0?* _ A202" , Aso?" as 
9 — 3. t p tee 




















10. Assuming ra cae {їкї кама. prové that 





in all cases except whore z—1. 
11. Prove that 





(P,A3 Руд + Ps3 — ТО, it та. 
{с. U. М.А. ё M.Sc. 1036] 
oaa I PAO”, show that РТР: 


+nP? 





2 
Фф чожши+ utem وس + و‎ e "Ано +. 
8 


sre. 
* 


















+ se Matus Р, 





Чектен ды а 


Go) wo tzAu, B A 





> mu, жадай, 1+ atu a+ 
16. Find w, from the following relations 











Gu,= 1-7 ¥1 


5 + Gu, fte 





16. Show that Glog LE yas = арм, 





= ‹соаб—# 
Spon; ut FEM 
Show that Вт +1, п) (тат, where m is а 
integer. 
Shew that 
рии Anu a SAP, уж. AT Tus „Аби, 


| и. A 4 М.50. 1938.] 
ЕЯ Sum to т terma 


1.2A2* 2,0 A72 8.4 Дз" — -4. ا‎ 
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39.. Shew that Afar + b)" = am (ae + bme о, 





Alex +b)" = —ain(az +Ьу-==- m. 


29, Sew that . 


fei dia 
Ld 27 Га + 1jcos- Ж" y 


Whero 2n > m > a, all being positive and A operating on т alone. 


Hence shew that 
n 2nz siria 
= ioe 


is constant for all values of т between qd =. 














24. Itu, a r tage али? +......ай infinitum, 
show that it can also be expressed in the form 
ye Ana? , Atat 


“mi-a 0-23" =a) 





A operating on a^» only. 
- [C. U. М.А. & M.Sc. 1939.) 


* 25. Show that the differences of zero 

4*0*, A*0**1, A*0**9, 

form » recurring series, and find its scale of relation. 
26. Shew that 





E TU EC » 
(a= ‘Sta dee 27 
where m is a positive integer greater than a, 2k the even integer 

‘next greater than а + 1, and A operates on n alone. 


„ 397. Prove that 
1*4 9*2 8927 + ...ad infinitum 


XAR {го «E am za Anne. 
aaa Tz ® 











CHAPTER П 
FINITE INTEGRATION AND SUMMATION OF SERIES 


2.1. Summation Problem., The proble 
Summation consists in determining a funo! 





of Integration or 
m ọla) which will 








Apiz)=u., 

where u, is a given function of =. Thus ф(х) i» such a 
funotion of = that, when it ix differenced, it becomes equivalent 
to us. We are to find the value of the unknown function 9(z) 
whose difference i» known. If we operate by the inverse of A 
оп both sides wo got 


(2) (= Aue. 

This operation is known as Summation, for, 
presently wee, it consists in finding the sum of terms like ty, 
Tt in analogous to integration and for that reason it is also called 
Integration, We shall bo concerned mostly with Finite Summa- 
tion of series. 

From (1) we have 9+9) =u, 

Hence 9(2)—9le-1) ч.л, 

Ф001) — 9-2) =u 

















ei . 2 
giro n*1)- ez n) u, s, 
Summing these we get Е 
pz) ема + uua t Mena es 
a = ч 
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"his sum is called the symbolic solution of (2) to the left. 
If the series converge uniformly in пу region, it represents 
an analytic solution of (1). 


‘The equation (1) is also satisfied И we stop st a particular 
term us, а<=—1. For taking 


Pa) =U ey ture t-te 
we have Ди, eu, asa m (uri Mara cM) 





Tbus (s) m Su, = S'u, 





It iv usual to represent the series 





u, by the symbol Xp. 


Thun in Finite Calculus 

Su, mu, g t urat 

the series beginning with the term u, ,. Hence 
в Ф) = A7 Tu, Xu, 

=. 


u.a tua- 


Thus A and X are inverse operators. Sometimes tbe 
symbol Sy is used for the inverse operation of A. Thus 





ole) А-а, e S ue. 
Again from the general relation (8) we have 
er) — ea) usg Melgar 
or eia) e Au, = Su, +96). 


„ It may be thought that corresponding to different integral 
“yolues of a, we shall get different values of A-'w,, eo that 
the numberof expressions denoted by A-! ie infinite. But 
that is not но, ая i» evident from the fact that g(a) representa 
the sum of the terms 


Many + tag May 
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The value of ф(а) does not change when = is increased by unity 
хо that Ae(0)—0. Hence we may regard it as a constant’ and 
denote it by C. Thus 

(6) ¢) =A Xu, +0. 





‘This is independent of апу particular value of 2, and represents 
the sum of all terms of the series beginning with м The 
difforence between the forms (5) and (6) is analogous to the 





difference between the definite integral ("Кох and the indefi- 


nite integral (/(2)dz in Infinitesimal Calculus 
The expression for Xu, can slo be symbolically obtained 
thus: 
Хм, = А-и, = (6-1), = E! (1B), 
m(E-' + B-2+ E 4 


I m 





Here X stands for the general sum of all tbe terms. It corros- 
Ponds to a=—a. То deduce the general form (6) we must 
add a constant C and making Ce —Xu, we have the sum for 
the selected group of terms wu, y иаа. Man 











Another mode of obtaining the sum is as follows, -We have 
ple * 1) 7 (z) =u. 
Hence +#@)—ф(=+1)= —u,, 
$(z*1)—o(2*2)— -м,.1, 
m 92+ 2)—9(2+3)=—uy sg, 


glen) фі=+н+1)= —u, 








Adding we have 
(2) = —u,—wu,. мв 





© Eee 
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‘This sum is termed the symbolic solution fo the right or the 
Principal summation. by Nórlond. И this series converge 
uniformly in any region, it represents an analytic solution of (1), 

2,2. Ihtegrable Forms. We shall not “be able to find the 
sum for any and every form of u,. Only for certain special 
forms of us, it will be possible to obtain finite expression for 
the sum of the series Su,, and we give bere the most important 
types of such series. In other cases, the sum will be expressed 
in the form of infinite series (Chay. IIT). 

1. Factorial Expressions of the form =" and (ox + b). 

We havo by $1.7 

Arm = (m+ 1)2™ 


gen 
mei 








Hence 
С Ima 





+0, 





ie, Xr(z-1)(2-9)..(e— m4 1) £61) (Е—)..(е-— т) +0. 
m 

In this indefinite summation C may be any arbitrary constant. 

То determine C we shall take X within certain proscribed limits. 


Thus taking the limit between m and n (n > m) we have 
m Y). (n — m) (n — m +1)...( —2) (0 — 1) 


«o p unte ct cel. 








Or we can take the upper limit only and determine С from 
the first term. Thus putting n for x we hi А 





завот D (n= 2)... (n-m) e, 


mt 
and putting n=m +1, we get the first term 


_ (m+ Dm...2.1 


H 1.9.8... m CI 





so that in this case C=0. 
Again ” A(arib) m P m a(m +1) (az +b) 


@ = Sox b) = ع‎ + bm tm 


ama) t0 
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‘The operations can be repeated. > Thus 


(az + by 





s om 
ө) Жа e by (mT) rT) * C1 Ce 
for ХС, = Cis. 

Tn general 
a) Stax + by = 








(m +1) (m + 3)...(m n) 

+з Сун" ++ C n On. 

Example. Sum to т terms the serios 
1.8.5 +8.5.7 + 5.7.9 " 

Here the nth torm м, = (2n —1) (2n +1) (2n +8). 
‘Hence the required sum = Хм... у 

=%(@н +5) (9n 8) (2n +1) 

m (2 +5) (2n +8) (ди + 1) (n 

" 8 

"To determine С, put n=1, then 








155-7 94.0, or C=, 





) (2n + 1) (2n +8) (2n +5) +15 
үн З 


- 2.8. II. Factorial Expressions of the form =‘ and (ох + b)i) 
We bave by $1.7 
Aziat (m+ 1), 


aS 
dito ^ 
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Again (ax + bY = a(—m +1) (a+ 0) 


2 зову Et SD 
e ^ eene ay жа 


Repenting the operation we have 

(ах +b) 
Siam 
То general 8 


(az + by-7* 


Share ион я) 
4 Cyan + Cur 





(3) — xX*(ax by jets 








+...+6. 12 Ca 


In order that the summation may be possible we must have 
т" +1. Otherwise the denominator in the result will vanish 
giving no finite expression. Formel. we must have mzs2, wo 
that there must be at least two factors in the denominator of 
the function to be integrated once. 

Tt may be noted that there results are all obtained from 
the results in $2.2 by putting —m for m. 

Brample. Sum tow term» the series 








нк» the nth term = TZ 





5 
Hence S.— Xo ту (8% 4) n 7) 


1 
C-ST +1) Gn 74) 
* Putting n1, we have 
1 
170-547 °"' С 





а 
30 San +1) (8 +4) ° 
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Working Rule for finding the sum to n terms of S(ax+b)™, 
™ being a positive or negative integer (Boole) : 

Write down the nth term of the series with its factors in 
ascending order of magnitude, add ono factor at the end (take 
away one factor at the beginning), divide by the number of 
factors now remaining and by the coefficient of x in each factor, 
and add to (subtract from) a constant. 

‘The bracketed portions in the above paragraph are to be takon 
when m is ре; е. . 

2,4. IIL Polynomials. 14 the rth term u, of the series 
be а polynomial in z, we sball have to transform u, in a serios 
of factoriais by $112 and then apply the summation formula of 
f22. Ву 91.124) we ha 








2 
u, = uo + жАмо + aT Аш, +... 





w > wetem e ыы Aug + ЖҮ Atuo +... 


3T 
‘The number of terms in (1) will bo finite, for u, is 'a rational 
and integral function of ж. 
Example 1, Find the sum of n terme of the serios 


134224 92... 





Here Elm + 1) 19 +29 + 88+ +n, 
and Alm +1)# = (n +2) — (n+ 1)* 9n 3, 
A*(n e 1) =2. 


о Sm + =0 tn hDg, "О-В 


acy "O+ n +1) 
shen Ges) | 


Putting n —1, we get C—0. 
A direct proof can be given thus: 
Хы, =мо Huy tugt. tue, > 
=а+Е+Е*# +... + Buy 
ажат, 
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"Therefore. 


E. 
o зч, f Ed E 

Here thé undetermined constant has beco avoided and the 
operations performed on the first term мо nre direot, We have 
only to construct the table of differences for tho first term. 
Thus for the above example, we bave the table 


Auo „Atuo Aup 





1 мон 
з 
2 ма -2 
5 0 
8 шо 2 
т 
4 «$4716 


Substituting for the values of ug, Aug and Aw, in (2) we 
get the same result as already obtained. 
Example 2. Find the sum of n terms of the series whose 
nth torm is n? 7n +4. 
Aw, Auo Аз 








uo 12 
t м 
u,=26 12 
26 в 
ча=52 18 
44 
из=98 
14n(n — 1) , 12nín — 1) ( 
Hence S,-12n4 7 S. + э 
+ n(n =1) (n—2) (9—3) 
at 
В = Lin +2) (nt +15). 
2.8. IV. ny Rational Fraction of the form - a 





where ч„=о®+Ь, and g(z) a polynomial in ж of a degree lower 
(by at least two unities than the degree of the denominator. 
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"To sum the series we express ¢(z) in the form 


92) = ay tarus + 





РРС 





where the constant д, a}. аз. ... are to be determinél. Then 






fx) 








Cy pa 

Each of the fractions can now be soparately integrated. 

Again it may happen that some intermediate factors’ are 
absent in the denominator. In that case we can supply the 
missing factors both in the numerator and in the denominator 
and then the factorial character in the denominator being 
complete we can apply the above method to obtain the sum. 


Ezample 1. Sum to n terms and to infinity the series 





14 18 
28458456 * 





4+6 


Here the nth termo СЕ TEY 


= a * 5 5 
WI) + 3j +8) wine) 305) 


: 2 6 
n iam e. (% +2) (» +8)” ЗП) 


Е 





10 229. = 
Potting n=1, we have ао =€ = g g" whence С 





and 5.=2 





Y 
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Example 2. Sum to п term» the series 





(C. U, M.A, & М.Бе. 1939.] 





+2 


Here the nth term = "+? 
n(n +1) +8) 


= 59. 
nin +1) (n2) (08) 
= "++ +4 
тт +1){% +) n +5) 


1 8 
22 (x8) * (13) 02) 08). 


4 
Чи 


ма ЕЕК... 
Hence Sn = 0 Зы ЗОО 


Gn? +27п +20. 
бз +1) (+2) (һ+%)` 


6+27+ 








Putting n=1, we have 





Functions of the form а” and a*9(x), where ¢(z) 


is а polynomial. х 
Aa” = (a — о". 





. 
In general 
ea хе gr Orate ae +6.=+6.. 
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| Again — Xa*e(z)— A-ha" p(x) $ 
(E 1)7!a*e(z) . 
Е — 1)7! ez), by $1.11(2) » 
=a" (o(1- A) 1)-19(2) 


аә 











"Thus 
(8) хатаа) 
_ Repeating tbe operations, we obtain am expression for X*. 
"Thus, 

X*a*e(z) = (E — 1)77a*(z) = 


328600) og +0. 





etz)-—. کے‎ * 














“өз ГА} нө 


© cg uj “үзмө +" (Sa Ганно} 
4 C 771 + Cua 4. + Du ya One 

Ла particular putting а= —1, we have 

C077 fet Тао Taten - c. 











20-0") 
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Applying (3) we bave 
8,-8*(2—1) + Xi*(22—1) 





“=8*(22—1) +S ios-n-3 2} +с 
=з'(в=—в)+ €. 
Potting х= 1, we get C3. Thus 
8,28 (5—1)49 
S, is also obtained by evalunting Sus. ү  X9** Mr 1). 
Example 2. Find the value of S*8** 1 (2r +1). 
Applying (4) we bave 








X +ع‎ (= 





(0+1) 7 2] 





+ OE e Cant +... Cu an C, 





(2r 1—8n) € Ce") +... +04 m C, 


where the C's are to be evalunted. 
Putting n=1, we obtain the value of 8,. 


2.7. VI. Functions of the form — e(z)/(z), 
where p(z) is а polynomial and (2) isa function of z admitting 
‘of repeated integration X* for all integral values of n. 
We have 
рее) = (RE! 1-1) = (A + BAN- (оу) 
where E and A operate on p(z) only and E’ and A’ оп yiz) ов у. 
Then 


Soll) = Еду (1+ к J eve) 











“Er f» 1-357 * gua free 
Е чы — els 2)А-542) 
жаз 8)а (а) – 





dropping the accents of A ав no longer necessary) 
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0) =e(e-NEyz)—Ap(e—-2)E* (x) 


+ A"p(a- BY) oe C. 





We need not add arbitrary constants for each sdmmation X 
оп the right side of (1). It will sufico to add one constant С 
in the result. Аз p(z) i» а polynomial the series (1) will 
terminate. 


Again 








Solel) m gelo (1+ gy ) о 





a 
1 a a(n + 1j ^з 
== f$ +З AT TUS jean 
(2) ен") -nAe(z п —1)E"* ys) 
» + ED Assam WE" fr) e 
Суат Оззи +... Cu cit Oa. 
A slightly different form proceeding by successive differences 
of (7) can be obtained as follows 
Ta ent 





-2{ " pom mv 





(8) а) — Аа) 9 (r1) 





Similarly, 
(4j Se) 900) уа) - nA (r) X" * Ve 1) 


= 22.0 аза)" еж)... 


4 O + Cau aso Ca¥ Ca 
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Example. Sum ton terms the series whose nth term is 
(8n + 5) sin (an + Б). -* 
Неге 8. = (Bn + 5) sin (an +b) + X(Sn + 5) sin (an +b) 
= (Bn + 5) sin (an + b) + (3(n — 1) +5} sin (ап + b) 
+ A(3(n —2) + 5)X* sin (am b) C 
ne 
(=) 
bir) 


_ 3 sin (an b a s) 
(2 а 2 +O. 
2.8. Partial Summatlon. To prove that 

Ж(и„Ао„) =ч, Fws., Au), 

Wo bave Аир. m uu AD, + Ursa. 

Hence operating by X on both sides we got 
7X(u,Av,) + X(v. 
©, X(u,Av,) = u,v, — Хе. 
Example. Find the value of Xza*, 
We have Aa* =о^(а—1), 


Hence  Xzo*(s—1)-—2a* — Xa*'! = za". 








= (9n 5) sin (an +b) + (Bn +2) 














2.9. Miscellaneous Forms. When the nth term of the 
senes proposed for integration does not fall under any of the 
given forms, it is often possible to conjecture the form of the 
sum by inspection from в general knowledge of the effects of 
the operation A on it. Thus if ¢(z) be а rational function of =, 
then ¥(z), defined by the relation Аа“ (=) = a(x), will also be « 
rational function, 

412868. 
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Ezample. Sum when possible the series. 
зз: 
Mri 








atat nat 


xs 87 ey Ferar” 


As the summation has reference to m, we can legitimately 
assume 


sa a 
(+1) (n+) 


Operating by A we get 
nts” atu an +0+6 ап+Ь 
deeem {= "+ } 


2 n+l 
B Jn? + (2a + b) (к—1)п + (a + b)z—2b 
| о RFD 
‘Equating the coefficients on both sides, we must bave 
a(e—1)=1, (Qa«b)r-1)-0, (a+b)z-2=0; | 





an+ 
nai 


















‘Thus the proposed series can be summed И ғ= 4, and then 
naar ^ »-2 


STD 


E 
ET М 





+0 * 








Putting n=1, we get 0=. Hence A 


mel 





үч s.=4 ва. 2. 


а. 
№ 
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2.10. Condition of integrability of the function 





The above function is not alwaye integrable in finite terms, 
and here we try to find out the condition satisfied by p, q and r 
so that it may be integrable. 

We know that 


Atan-t 








ч а 
[Fab 685099 + Dade xata 
1 
ал 
тетт rz” 
Hence comparing the above expressions we have 
1+ab+b*= —ap, a+2b=~q, а= 











or q?—r?=4 (pr-1), 
which is the condition to be satisfied that the given function 
may be integrable. Thus 


а Xtan 





Example 1. 


1 
i. 
a sy 





Here the coefficients satisfy the condition of integrability, 
Hence by (1) 
e 1 
вети 








‘The condition of integrability is satisfied bere. Hence by (1) 
= 1 
Som ee" улу eg n 
=! 2 E 
ты 


Putting 21, we bave tan7! 1 = C—tan- j, whence Co tant). 






Hence Spaten i “ھا‎ Samtang 


24i. Another method of summing such series can be 
_ ebteined thus: 


VIMITR ISTRORATION AND MCMMAYION OF «anme LJ 
Beample 1. Evaluate Xtan t Ta 


atts ды 
"rapper и EE 
=O + tan” ele — 1). 


Bromple 2. Bum the series whose sth tern ie 


Now Iter 


зе ПЕ Ы 





= 1 
Bere Sive qo i-i edi 








Xa a eT or e" 
Hence 8, C „маг! (100-7). 
Potting z1, we have C tan-! ( 8) — tap! (S) tan" (7). 
Thus 8,» tan 7 (7) tan" (10s — 7). 


"The same result can also be obtained from (1) of $ 2-10. 
212 Sum of the series 
ay 1e ay 1 + +e 1 2*2 09) +... 
-..2- ++ 





when the sum of tbe series 
NO mas анат +в,” +... 

de known. 

* Here the general term ea, (1*4 2* +... +a)" 

+5*+..+5*}0* 





8 9 
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98 (Et —1)07 
Ы = fare- antr on ^ 
Hence the required sum 
A =0} on 


Ezample 1. Sum the serios 





194-95 1849548", 184894. tet | 
вт ec чар: 


F 179^... 
Here the general tern Fe 
=4 Egon 


(the operators having reference to 0) 
О + 8 + ..(-)1+ + ait Jo 
= a ue E xi 
Dad 
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Here the general term = (1^ +2"  8* +... + z*)t* 
EN (gs quas B gage Ре 
EU 1)0 АЕ t*).0*. 
Hence the required sum 


E 
“Хх faciens... 





ажене. 





p 
-Е {аву о0о 


a MEI е Жк НЕР 
Еа е 70 





АТ 
a 





MISCELLANEOUS EXAMPLES 


Example 1. If X^, denote the sum of the first n terms of 
the series uo, tig, t4, 





E (Was—wo). Boole.) 


XI us m tg tg tenet pug 
2 Xu, uo tty Fg +... Mn a ча, 









Xue) + (но ми +p 
* Now uo =u +ug—... + Mang sa1 


--i0-$ кез r ман — мо). 
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Hence Yu, =} {2-3 +ô- = deeem. 


Example 2. И Ди, =... —u, and A= Shew tbat 





ч, FAMAM, HA A up +... HAA U 
=a (a7 — хаты, «Aat Anus у. 
C. U., M.A. & M.Sc, 1936.] 





We bave 


u. +АДи, + АЗ Зы, +... FATA =}. 





= (I-AA tu, -A* * 11 -AA]-1A**1u, 
= (a* —1)(a* E—1)-!w, + A*a* (a* E—1)-1A* 
=а- а= (Е) а (Е 0) Зандан, 
=a * (at хати, +л"Жа' = дити 











EsampleS. Sum to terms the series 
£21 ТЕ че E de HP 
sin 6 sin 26  sin26 win 36 sin 36 wm 46 








1 
ee Sem Sura Ха FT win EF " 


sin @ 
Per 2000 = аб ain FFI 


2 Aet (e+ D0 a E+ sip EFA 


"Hence вже. | 





1 cot 20 
таана 21, we have C= ag уыл 


gueule, 


Ee 
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Ezample 4. Sum to terms and to infinity the series 


1 лы 9 
ton 6+ z tang + 4 tan T 








Putting z=1, we get C= —2 cot 26. 


Hence S} tan в 12049, 





‘Sum to infinity is 3-2 cot 26. 
Example 5. Sum to terms and to infinity the series 
м1 EY ‚А 
tani} +2 tans + 2° tan? ggg +- 
The series can be put tbus 


ml а e 29tan-71—— 
tan +2 ва”! gigy treno raa +... 





4x8 ai 
- ZI 
Hence tee m tanc See eae 
and 8: = Зи, = 2" tan! 
Ал” tan) Laar lan! cu ТИНИ 
H s^ 4.843.857 





= 1 
S,=C+2* tanga: 
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Potting 2—1, we have C—tan-17 —2 tan = — = 
хо 8,=2* tan 1, -1-£. 
a 27 tan“! à 8.-1-1 
Ezample 6. Som the series 
Ae d 1 
1 PU e Z 
tan rs #2 еп”! ^° te раг А 
ача 
Hore صا 2= .ن‎ ты, 
Now ар-11 tan? gk =Фбав-!,Ь. 
n dn? + 3" ЕЛ 
т = = 1 ail 
s. tan m 2 tan tancii, 








1 





Hene 8,93«,., = 32° tan ost gary 





=32" а امھ‎ 
эз“ (2 “مه‎ Se tan zz] 


E 





22*tn-1 5 +0. 






-11 
Putting z=1, we get C= — tan? $ 
Ei Lock. 3 що 
Thos) B,-2' tàn йал. — Sen баве 
Езотые 7. Sum the series 


tan uy tan ууу ted 
where /°(n)=n, f'(mm4n?3n, /%(п)=4{(я)}* *Sf(n), eto. 
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= = 1 
Here aea m (D) lan ау 


^ 1 1 Ы 1 
zl 1 — = i 
Now we have tanii аа, =2 tan 


Hence putting /*(») for т and multiplying throughout by 
(—1)*, we get 








(=1)* tan" zg + (—1)* ten“? 





f) 

AA M iE Ta > 
oF A(- 1*3 ЕО 
whence Bue =F (е ten + O- 


Putting z=1, we got C=} tan-! T 


‘Thus 8, 





д. (e +(=1)**1 tan“? 





Ps 


eem 
8, * we. 





ExaMPLES 
1. Sum to z terms the series 
@ 1.8.5.7: 8.5.7.0 5,7.9,11 +... 


[2] 





CALCULUS OF FINITE DIFFERENCES 
10 12 14 
© stas ruso 
(vii) 1.2.442.3.5+3.4.6+... А 
8 Га. СРЕ, Д 
7.13 


+... 


(ей) ї 





(=) 


з 5 а 
258'5811 Emma h 
(zi) 1.8.5 cos 6+ 8.5.7 cos 20+5,7.9 сов 86 +... 
(zi) 1+2a cos 0+ Sa? cos 204+ 4a? cos 364... 
2. Sum to terms the series 
(0) 19,24 22.22 482,95 «42 264... 


(ii) 1.82 9,52 +5. 








8. 
© 
1 5.1 > 
з 1 4 б 
TIT + 2573 * 34 35*- . 


6 1 Y x 
UD Ba ава вав" 


| 
| 1 1 С Rim z 
| (©) a+ V+ a+ Des VD + ауу *+-- 





Й 
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4. Sum to x terms and to infinity the series where u, is 
equal to i 


ed: 
69 ten ser eg 








Gi) tan су (iv) tan~ 


1 
-—12:-18 " 2627 — = 4" 


5. Sum tox terms and to infinity the series 








(ii) sim. 8 ( 


“= 


(e) (sin өе +4( sing Ks “( A ys 








6. Sum to z terms the series 





« 7. Sum to x terms the series 
(i) sin 6+sin 30+ sin 50+... 


* (ii) сов 6 +008 30+ сов 50+... 
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+: 





OUS ы кл ry 1 
Sos 6 сов 26 * роз 28 сов B * cos 89 cos 40 


OD ее аә 





i 1 Pie | 1 
bed CLA EE] 


tiin (aa) (Ea): TE (ыы 


B. Sum to = terms the series where м, is equal to 








( 3295. 


6) 2° +78. z+) 


: (e+)? 
te ® Seer 09 DES rm 
9. Find the sum to z terms of the series 

G) 2.24 7.44148 +23,16 + 84.82 +... 


id 
EU 








6) rm 


10, Expand Xeí(z) cos ox in a series of differences of (z). 
11. Shew that 


Mg 4, cos 29+ мысов 40+... o o шу + Aug, enin о 






жыр con 26— A win BOF... - 
19. И ооо + tari yr) +... shew that 
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ава И. фо) = фора +0924... shew that 
DOES 





¥ чова + аде) Aug T7 ра) Ado +... 


‚ [Guderman, Crelle, УП, 800.] 
18, Sum to infinity the series 


or 417.5. ча 2-1 


Hanc 
z 5 


) 

14. If 1542743" e... m" m8, and mm +1) p, sbew that 
S, p*f(p) or (2m +1)pf(p) according ая п is odd or even, f(p) 
being a polynomial. 














[Nouvelles Annales, X, 199.] 
15. Shew that 
1 1 1 1 12 
n'est erg "m 


ЗИ t ЕЕ) 
10. И бз) meg +z Far? +..., shew that 
ам, г" +a 


ими e SEI ее а, cantly cant” * ° +... 





-1 [steer exorto) amen.) 


whore a iw an nth root of unity- 
17. Prove that 





9680 —х^Ли„) 


+ sito 27a) e 
18. Prove tbat 
a) =a) 
And hence sum to n terms the series 
148.214 7.81 € 13.41+21.51 +... 





CHAPTER III 


APPROXIMATE SUMMATION OF SERIES : 
BERNOULLI'S AND EULER'S NUMBERS 


3.1. In the previous chapter severo! methods have been 
given for finite summation. But when such finite expression 
cannot be obtained, we can find out an approximate expression 
for the sum. In many cases the difference botween the actual 
value and the approximate value so obtained will be very small. 

3.2. Euler-Maclaurin Sum Formula. Expression for Xu 
in terma of the differential coefficients of us. 

We haye Xu,= Au. Ce (60 —1)7! u, C. 








Here -4 is the only term involving odd power of t. For 





t do. E ا‎ E a 
we have 5 жа ear tie тю? 





and these expressions remain unchanged when f is changed 
into —t, во that we are sure that the right-hand expression 
contains on expension only even powers of t. 

It is usual to express (1) in the form 





٤ م 8 م کے‎ Ва ہم‎ В 
e атлы a aie pt 





«=== Bannann + s 


The quantities B,, Bs. Bg, ete, are called Bernoulli's 
mumbers.* Their values are р 
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‘These numbers begin to increase after В, and these ulti- 
mately form а divergent series. и 
Dividing both sides of (2) by f, we have 





Ba a, Bags. 


e 4! et 





Hence 


ву du, _ By du, 


0 uy = 0+ fusa Lus ик 4ш E ЖЫ 


1 1 du, 1 d*u,, 1 dw, 
он fuda- ен Е т 223 оно def 
"The value of the constant C is to be determined from the 
particular conditions of the problem. 
If и, be a rational integral function of z, the integral is аб 


once obtained. 
‘This theorem will be very useful when finite summation is 


impracticable. 
3.8. Bernoulli's Numbers. The numerical values of 
Bernoulli's numbers may be obtained thus: 














Let f()—— {| then by Maclaurin's theorem 


a) оло +(4),t+ (RE unt unte 
But by $8.2(2) 


fü-1-i mg 












e L 
Bap + Bog, 





жо that all the pdd differential coefficients in (1) except the 


(eet ›„=% {ос n=1, 2, З... 


6—1996B. - 


first vanish, 








=0, when n2-1, 





‚ when nel. 
= Putting t=0 and n=], 2, 3, 
n=1, /0)=1, 





in (2) we have 


n=2, 1+29 =0, or ( 3 3. =x 


т=з, 16535 107 





5 
ie. 1-3 





-%=(4 au 


Proceeding in this manner we get 





Example 1. Find the value of Sz? =19 422+... +(2—1)%. ^ 


We bave Xi9—€4 безас рана з 


292.2 
=045-3+5- 


Putting ==2, we get C=0. Hence 
дез ew ЕП. ЕЯ - b: 
6 


We have obtained here a finite expression for the sum. 
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Brample 2. Somtheseries Xi. , 


We have by $3.2(4) 





To determine the value of С put z=% , then 


el. .Rol1 1 ТРЕК 
E = +g + +... ad infinitum =o, 


But it is known that 


TOL а S 
{я * аз * з +... ad infinitum = Tr, 





If x be large, we need only take а few terms on the right- 
hand «ide. 

When the sum to infinity is unknown or infinite, we 
oan determine the value of C approximately by giving to z some 
such value which renders terms after the first few terms on the 
right side of 88.204) practically negligible. 


Thus taking the series 1+ 2 +1 











з 
1 1 
we have жа 
1 ee Se E 
rg + C+ utque ~~ 





© y 
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То determine C, ceo 
1.4. d LA 
1+3 +4 + REC + log 1045 — + 
Here we can neglect higher terms as their values are infini- 
tesimally small. Simplifying and giving to log 10 its value 
2:302585, we get C= 577215 о Hence 


8= 817218 + log «+ 3. 





1 
2 "uat 1202* 


3.5. Stirling's Theorem." Approximate value of 1.2.3...=. — 


Pot м,=1.2.8...2 
then log м, =log 1+log 2+log 3+... «log = 
log 2+ log =, 
| Hence by §8.2(4) we have 


| log w, log #+С + flog xdz— 3 log x 








а) 


To determine C, suppose z to be infinitely large, then. 
e log (1.2.3 ... а)= C (e+) log ==. 
(8) Therefore — 129..2ee 77,57 +4, 
а Hence 1.2.8... 2e 40 25 (got +4, 
м _ Again multiplying (8) by 2* we get 
в з.4.6...2==2°‹©—** +4, 
К. ‘Dividing (4) by (5) we get 
© 1.4.5.7 ...(22—1) 9 9" * 3e 7 * P. = 
_ Again dividing (5) by (6) we have 


ia Te | CA МЫ: 2 
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But by Wallis’ theorem, * since z is infinitely large, we have 


24.0..(27 9) (2x) мя. 


- уз” 





Hence from (7) we has 

y (22) y e e ve 

or «6 = y (22). 
хо C=log 2s. 





Hence 


log ч, log V (2a) log a” M ±+ а 





„a= y (ur) = 


zy 1 1 139 
@ = vee (= ) (1 * цы gases вис * 7 )- 
For infinitely large values of z we 





0) 1.2.8... 27 y (2==) ( = Ja 


Ezample. Find an approximate value of 





when « is large but not infinite. 
Aw з їв large but does not tend to infinity, assume 


123..5- (ar (= (1+ 
meglecting tbe higher powers of T 3 
. 


+ Wallis-— Arithmetica Infeitorem (1024). 
“+ Anelementary proof of (10) was given by Glaisher in Quarterly Journal 
of Mathematics, 1878, See niso Raabe, Crelles Journal, XXV. 
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Then 128...305 Ме. Z J (1 i. ) 
and 24.0..24—2* vem. ( (1*3 


sotbat 8.5.7...(2z—1)= v2 ( = yc le de ) 


As veal) as) 






(22-1) (2z+1) _ 9541 ( 1 
m FI Вг. 





8.6. The cbief use of Euler-Maclaurin Sum formula §8.2(4) 
i» to convert a series into another which converges more rapidly 
than the given serios. Tho series so formed is convergent for 
the first fow terms and then diverges, but by keeping only the 
convergent portion we obtain an approximate value of the 
given series. It may be necessary in some cases to find the 
sum of the first few terms in ordinary manner before we apply 
this formula. This gives а very close approximation to the 
value of the function. 

As an illustration take Ех. 2 of §3.8. We have 


В, B 
рот Ват 


Taking the series between the limits co and 1, we have — 


+++ ad infinitum =1+ 5 + В Bs + Вь | 


As the Bernoulli's numbers are convergent up to Bs and. 
then diverge, let us stop at Bs. Then 











Here the convergence is so rapid that if we stop at By 
we have 


1.305 ,1, 1 , By 104494, 
хатта 5 oF bo - 
which differs very slightly from the teue value 1164403. To 


get а similar approximation from the original series, thousands 
of terms would have to be taken. 


ЗЛ. Ifthe terms of the series be alternately positive and 
negative, we can make use of the formula for X(—1)*e(z) to 
convert it into в rapidly convergent one. Thus taking the series 


$(0) — e) e(2 — 
we have this equivalent to 
ü-E«E- 


. ad infinitum 





0900) 





=а+ку-уш= ( 1+ 








@ xcswe-i-2 a esum. 


The formula (1) transforms the series into a rapidly oon: 






1 1 
S-isxi6s3 GHEE) 77 





| 
| 
| 
i 
А 
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Applying (1) we have 





Dur $e exu t Ae {= +2) E+) 





1 84 
ао EFS EF ^ 
which converges more rapidly than original series S. 
In particular, putting 2—4, we have 


ае В 
4725678 * 


38. Expression for Xa*e(z) in terms of differential 


coefficients of у(х). 


We have Xa*o(z)- (E —1)7!a*e(s) =a" (a — 1-14 (2) 
=a" (aP —1)-19(2). 
Now by Herschel's theorem, we bave 





ве” = GE -enm [1+ Ape 
where A operates on 0 only. 
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= (—1)*-! (ео —1)7! —2 (020 —1)7?)e(2) 


-142)* - By Bs ps 
-(-1) (5-543 1-р--Җа-Рэ+... 


1. „А Анн 55H, 
-a(z 50-46 wp? +... )} em. 
Thus 





® seneo [1—4 en 


B = 3- 
+ 0*-00 





where B,, Ву, eto., are Bernoulli's numbers. 


3.9. Expression for Xu, 
of uy. 


Since ="=(e?-1)-* 
we shall expand (e^ —1)-" in ascending powers of t. 





terms of differential coefficients 





Put vaml 177 ; 
AN de. —mne* 
— 
"Hence ма (т) 





Similarly | (n-Dvz- -( dti) 





2). = * Р 2)... 





1e --( 4 +h) " 
“Multiplying these together, we have 


ден n (а, а) (ана) +1 Ye 
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zo (e 1) Ea $2) (2+1 
3 at "X at 
e 1-1 led ta Pss | 
61 





Putting D for ¢ in (2) we get the expression for X" in terms 
of Bernoulli's numbers. 





Example. Find an expression for X?u, in terms of D. 
We have from (2) 


e-i 2X3. Xie Be Baers... 


(ананан per) 


аазы! ‚э ВВ _8.2Byt_ 8.3.Bys | 
zi atii ent D» e. 


Mes T 52 8u- 2 a 
Boate e 1340 160" * 








Hence 


X*- Ума» рч” + фе i 


19 du, _ 1 due, 
240 dr ~ 160 dz* 


8.10. Expressions for Bernoulli's Numbers P5... 




















We have from §8.2(2) 
a) 1-1 Bre Batts Bete. 
Я + +: 
Again by Herschel’s theorem 2 
e wit e i 
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Hence comparing the coefficients of t?" in (1) and (2), we 
have 





—1)** Bas 
cante 






As the series (1) does not inelude terms involving 
(п = 1), its coefficients are zero. „Hence 


(4) p- EAN CANSA, -an онно, 












1 фу by $1.90, 


Giving to m the values 1, 2, 3, .,. and substituting in (8) we 
have 


Basan {(4-4 + 2r t Ir yr 
-(4 a 


But since ( a--+ 








= Jo 
=log (1 + Ajo = DO?" —0, 
we have 


(5 — Bas-1=(-0)" 





yen 
3.11. Another elegant expression for Bg,—, № obtained from 


the expansion of ; 7. We have 








03.12 
ЕП 


E a Вен ae Bp ene 


у Же eS, € 


Hence, comparing tho coefficients of ¢3* 








ene S mpi 












=en; 1 o-i 
e Ваа" aut “ga 7, 





а i-i Get Gre fom 


8.19. Expressions for tan 6, cot 6, eto., in terms of Baami. 
Wo have from $3.2 
еа pe tee е 
4 
_ Putting t=26 in it, we have 


hs son ecl. дане – pee о а 3 


"ga. 








Cay Bae 





i 9 
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Expressions for tan 0 and cosec 6 are obtained from the 
identities 


tan 6=cot 9—2 cot 26 


coscc $c 6 соё 6. 
Thos 


(8) tan s 27-9 B,64 2105 gs Р ыз РИА 


(EP p, rectus 


and 


(4) еоесб= у 1 «20-1738 «299 1) Ha 09 «as nete 


Py sett} 
Gn) 





Bgng +... 
Similarly 








(6) cosech ө= 1. ETD, + cmm вз - 205-0 p 





348. Agein since sin 6=6(1-% (175513 )-— 


wd heve log sin 010g ¢+10g (1—03) нот) 
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ИСТЕ ТЕЧЕТ Ари 








Aue +} 


Comparing the coefficient of 93" in (1) with $3.12(2) we have 
E] 


2 (олу 1 
Q io, Basi args fee 








"T 





iv {i+ ge tae +} 








Let 8, 





Then gis Seem i 





Hence (1-31, )Sa.=1+ gle + gar t gis teen 
rers dil tho terms whose denaminators are multiplas of 2 have. 
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1f we take the next prime number 3 and multiply 
1 А 1 
( - ge )s- by 1 ie. 


we shall doprive the series of all terms involving multiples of 8, 
and soon. Thus we get 


(1- he) lea ) (de ent die 


where 2, 8, 5, 














p are the succession of natural primes up to P, 


and q is the prime next to р. Taking q^ to be very large, “= 


p 


becomes very small. Hence 


э: [Gir Xo cte ig) 
where the succession of primes continues to infin 
= е {сус түр» 
А ZU (1m ate t= ate аа 
For large values of n this formula will be of great use. 

3.18. We have from Algebra 


е 1 range eae 
тз gae + 529 + aet 




















Hence from $8.18(2) we have 


Byss) +3) Qn + DS acs 
(2«)98,. 


Bee 
£4 
А ея ар асе) 


an ст] 


(п) —1 
ро а 
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Hence B». , > Bsn- provided 4n? —1 > 42% or n> (92 +), 
ic.,ifn>816. Now B; <B,. Hence 
<в,<В„<... 

Aguin since Sa tends to 1 for n—>œ, we have Baw-i 
tending to @ forn—>a. Thus 

Bernoulli's numbers are all positive and they increase after 








‘The value of E ultimately approaches to equality with a 


geometrical series having for its common ratio ae 


3.15. Again from §3.12(2) we have 





B -Bi gm — Bs забег 
Ө cos 6. sin oft 3i? 215 3 





zB. pr 
р" | 


Expanding cos @ and sin 6 and equating the coeficiente of 
3^) we have 


ET (-0"* 











m)! ^ + Di 
Бона ооа Baigta 
*в1 =й! e ^. 
‘This i» a recurrence formula connecting the Bernoulli's 
numbers. 





Euler's Numbers. Another serios of numbers closely 
analogous to those of Bernoulli was given by Euler." They 
appear in the coefficients in the expansion of sec б. Thun 

Е, 


а) seems Epon + Baars + Bar rns E 
Again we have 
کے کروی‎ 
у = 


* qnit, Cal. Dif., $224. 
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"o9 (ра: 


verita а Gastos 











where ot} iri 





Hence the coefficient of 02" in (2) в % > 1 





a ме 
BIET n? А, 


Comparing the above with the coefficient of of" in (1) 
we have 


O ван A (1 





1 
-1)°-1 а > 
+ ут 
This can easily be transformed into the form* _ ү, 
ELT А d МЕ 
те meee, т xr т 
= 


3.16. Recurrence Formula for Euler's Numbers. ` 


2a 
Since x sec б=1+ PE 


Ega?" 
(аи) j’ 











we have 1= cone [tes 






rj! 








‹ [r$ Ba өз ром Bas 





© Eulee—Inired. in Anal. Inf., 394. 
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АН Da "Gr =o, 





. acu 





n "C Es. +S * 
к *(-1*7* °C, Eg  (— » 

* тые equation enables us to calculate th 
s Thus 











values of E's, 





Y B,=5, E, 01. E, = 1885, Eo 700521, ete, 
PS Again since 


1 1 1 
11-е tg --21- ід. 





EXIT TEJ ж-з SA 
= $ “per pia] А 


En p 








gir]. 
constantly increases ns n increases and is greater than 1 
- forn=1. Therefore Bg <E, <E g.. 





Abo Lim fi gary tye =1, when n= 


Coefficients. 
ри usually termed a» | 


cients. The numerical value ot ihe eoetficiens c 
by С. Thus 
=O — C3. mH 
Ea 











T M «x xx 


~ 7 ы APPROXIMATE йомйлтюн “4 sa” 


By Maolaurin’ theore: 
0 42D ое кыш 













tho actual ovaluation of which involves laborious caleulations. — › 

We can however obtain a recurrence formula, from which 
the values of C**! аге easily obtained from those of с". ory 
we have 








aD (ene, à 
во that 1 
(3) of" +n. 

‘The values of C} are 1,0, 0,..., so that those of Cf are 1, | 








1, 0, Oss, for ke 1, 2, 8... 
- 8.18. Again representing by Cz" the numerical value of the 
coefficient of à in the expansion of zf" io negative powers of 2, 


wo have 





р SF +... 2 





X 2 (-g*71(0-1) 127m, 
z 





Ко ee ге) ы; 
EC. Ш 


E. ATi 
а C= — Ey 
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аа 


RTT ae ветена, 





Again wince A*710*-1— (ı_1){4*-30*-2  A*-104-2), 





from those of ©,- 
Again since 


W^ аточ адо" + —Vasor + 





we have by $3.18(1) 
| w ES [20:2 526 пог ud 





Similarly from $ 1.19 we bavo 


р-н 
IT 


127-4... by 1. 





+... 





o 








Ехлмріва 


1. Find an approximate expression for the series 
Lt he tyne etue 


(Go Find an approximate expression for 22, and also the 
value of P f 





Ittt 
о places: of decimals. 








APPRORIMATE: SUMMATION. в 
3. Find an expression for 


1*3 





obtain an approximate value for the sum ad infinitum, 
4. Show how you would convert the series 


an 








1 1 
+++ 


to в more rapidly converging series, and thus evaluate ite sum 


correct to 5 places of decimals. 
[С U., М.А. MEy ву" 
5. ‘Transform the series y= lst f=... into a rapidly 
convergent one. . 


6. И мочит + мых? +... o jiz), show that 


мониро оз + шаю». +... = [(1)е„ +РО)Аю, + Өля», ‚= 





and apply this theorem to transform tho series 
سي‎ iet Dom ae rm 


to one proceeding by factorials only. 
7 Find the sum to n terms апа to infinity of the series 








[C. U., M.A, 4 M.Bc., 1938] 


XA 1 
sate +H 
1 af к 


9. Prove that if n be an odd integer 
Bis в, + aen 0-2) 0-5) Bac. 


о-в, +... 





gon 


to (r — 1) terms. 
~ 





* 

> 

^ 86 CALCULUS OF FINITE DIFFERENCES 
10. Bhew that 


Bisa 





Sa jam 
11. Shew that 


1 1 1 a " 
++ eec infinitum 





»Ü-bgs,1 BP, B, LA. 
ж = Riga Ret kale oo 


[C. U., M.A. а M. 
° «ва apply this to the summation of Lambert's series 





when zis L pearly. 





(Zeitschrift, v1, 407] 


in terme of successive differences of 





12. Expand XC ру 
log (az+b) and deduce 


хой == {ток sin 2— 2 tog sin 22° tog sin 2 









L. 13. Prove that the coefficient of 8?" in the expinsion of 


еа) ‘equal to BT в» 


a ie 16. per Bernoulli's numbers or otherwise that 








Е 8 
Serer er 94 ncs 
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10. Find an expression for Взл-1 as а definite integral from 
the identity 





17. Prove that 
m. Dose (71). 
nj igi P ent авт Diet 84-9 temm T 


[Grunert, xxvn, 455] 


=0. 





18. Ws that 
i n 





" Вун? _Вут*_ 
` E 1—57 ROPA 
19, Shew that 








ы “22. Show that Cp-C3C3—..—0. 
and' В,0}-В,01+ BsCZ~...= 








CHAPTER IV 
INTERPOLATION AND MECHANICAL QUADRATURE 


| 8.1. Interpolation. In mathematical and statistical 
researches it frequently happens that certain valves of a function 
аге found by experiment, observation or calculation and it is 
Fequired to find other values of the function from them. The 
Process of obtaining unknown values of the function, or the 
approximate form of tbe function, from the given values, is 
known аз Interpolation. 
Although differing fundamentally in the character of their 
respective functions, all mathematical tables are alike in giving 
| the numerical values of the functions for certain assumed values 
of the argument so chosen that the intermediate values of the 
function may be readily derived by the process of interpolation. 
As ро series of particular values can determine a law, the prob- 
lem of interpolation Îs an indeterminate опе. The general form 
of the function is not known; only a series of isolated points 
| on the curve representing the function i» given. Between any 
| two consecutive points we can draw an infinite number of curves 
| and thus we may get an infinite number of curves passing 
through the given series of points. Consequently the analytical 
- expressions of the curves are infinite in number. It becomes 
| therefore necessary to make some assumption in order to bat 
| a problem with definite solution. The assumption usually made 
is that the function can be represented by а rational integral 
function of = of degree n —1, when n points on the curve orn 
‘other relations are known. ‘Then the nth and all higher differ- 
ences will vanish. Tt is with such functions that the actuary 
| and the statistician are mainly concerned. « 
Extrapolation is sometimes defined as tbe process of obtain- 


ie E emas а а serica tn орой ts керей widak 
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4.2. Gregory-Newton Formula of Interpolation. When 
^» consecutive equidistant values of the function f(z) are known, 
the general expression of the function is obtained with the belp 
of Gregory-Newton Formula. Let f(a), f(a +w), .... f(a + n — 1») be 
the given values, then by $1.9 we have 


а) fao) (а) + zAf(a) + = Юдзо) * 





267 1) 2)... +2) As 


amiji чө). 


As n entries аго given. we can find tlie differences up to A*-"f(a) 
from the table and A*/(a) and higher differences all vanish by 
assumption, во that /(а + zw) i» a rational integral function of 
zot degree n This represents an approximate expression 
of the function. 








Taking a=0 and w=1, we have 





(2) (2) = (0) +2 f(0) + € aepo) +... 





: „E-n tD e ago), 


Example 1. Values of z* ure given for 22, 4, 6, 8, 10 and 
12; find the value of (2:8)*. 
Completing the table we have 


= zs a as at at д» 
2 16 
240 
4 256 800 
i206 beri mo ® эзш 
в x 
2800 1844 о 
в 4086, 8104 884 $ 
5004 1728 
10 10000 4832 
» 
12 20786 3 № 

















"o Кельсоков оғ FINITE DivFERENCES 
Here a=2, 02, 2572:8—2— 8, so that £= 74; 


co (8-164 4 х 2404 SD 





71064 90—96--61:44—15:9744 
614050. 
‘This result i ly verified and found exact to the last figure, 
Example 2. ‘The Northampton 8 per cent. annuities for 
the ages 21, 25, 29, 33 and 37 are respectively 184708, 178144, 
171070, 168432 and 155154. Find the annuity for the age 30. 





a a at 
21 184708 
25 178144 — 
‚м - — 0054 
29 171070 E 
А — 10076 
33 1678482 
— 8278 i 
87 155154 ; 
Here a—21, «—4 and zw —9, so that 22:25. ы, 
7. 80) = 18'4708~2°25 x +004 — 2251325. оро 





_ 225 × 1-25 25 „соб 4 2°25 × 195 x "25 x 75 «0022 
41 





LU pi 

109216. Ц 

Exomple 3. From the following table find Sun's Right 
Ascension for April 20, 1898. а 

Date Бип» К. A. A ae at A 


1898 m. м. м. в. в. в. в, 


April16 1 38 1959 














. quantities. > 
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Here а = April 16, w=5 and z= 8. 
n 


м. в. п.м. s 
-. Sun's В.А. for April 20— 1 38 1959= 1 38 1950 
+080 х 18 3695 +0 14 59:000 
-008 x 1099 — 0879 
+0-032 x 105 + 0084. 
— 00176 х ооз — 6000 


+ = 14591176 
- Sun's Н. А. for April 20= ` 1 58 11-75, 


which is the value given їп the American Ephemeris for 1898. 

4.3. Missing Terms. Suppose the given values of the 
function do not form в complete series of equidistant values in 
consequence of the absence of one or more intermediate terme. 
Let it be required to supply the missing terme, 

Let f(a),f(a--w), f(a--20), ... denote the complete series, of 
which m terms are supposed given. If one term is missing, the 
number of terms in the complete series № п+1, and to find 
the missing term put A"f(a)=0. Then by 61.10 we have 


Q) f(a + nw) -nf(a + ZT) MID + (=a) =0 


f(a + — Зы) 
from which any missing term is expressed in terms of the given 
n values of f(z). Thus to interpolate aterm midway between 
two others we have 
а +2) —2 f(a + =) + fla) 0, 
so that Га v) Ма) + fla + 24)). 
"To supply a middle term in a series of five, we have 
f(a) — 4f(a + ө) + б/(а +24) — 4f(a + 3w) + f(a + 4w) =0, 
во that  f(e+2u)=4[4{fla +) +/(а +8} — (f(a) + f(a + 49))].. 
Iftwo terms are missing and п terms are given we must also put 
A"f(a- 9) 0, F 
" Mag A Yu) пата) да + = кй 
e CTI) fla a) =0, 
veda e un oS гт элш unknown. 


























vs 
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In general, if r terms are missing and » terms are given, we 
require the following r equations to evaluate them : 
аа +т—1ш)=0, 721, 2. 8, 
Example 1. Given log, o121 — 20827854, log, o194 20934917, 
log, ө122 = 20803598, log; o125— 20069100, 
find the value of log, 9123. 
Putting A* f(121) =0, we have 
log, 0123 = MA(log 9122 + tog 124) — (log, o121 + log, 9125)) 
= }{16°7191260—4°1796954} 
"0809051. 











Example 2. Complete the following table 
z-20 23 22 23 24 25 26 
Nz)= 185 au 100. 1082 074 


‘Here five entries are known and hence the equations are 
f(2:5)—5/(4) + 10f(278) — 10/0272) + 5f(2:1) —7(2°0) =0, 
1676) —5f(2:5) + 10/(2:4) — 10f(2°B) + 5f(2"2)—f(2'1) = 0. 

Substituting the valves of f(s) and solving we get, 
К 1071) = 7123 and (274) = 090. 

4. Subdivision of Intervals. A frequent problem in 
actuarial work № the interpolation for values of а function f(z) 
between every two consecutive terms of a given series. For 
example, premiums of policies maturing at every fifth year or 
tenth year being given, it may be required to find the premiums 
for every consecutive year. ‘This operation is known as 
aubtabulation. А 

Let the given series be f(a), /(а +w), f(a +2)... and suppose 
that it is required to divide the interval » into m subintervals 
so that (m —1) intermediate values are to be interpolated between 
‘every two given consecutive entries. 

We have by Gregory-Newton formula 


p(er )= rhe as ano 





m 
} f =H + kaosa aCi AH ; 


5 Y =. 


м ae ] 
р. 2 H 
[1 E а: л A м. 

















i 
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and putting k=1, 2, 3,......, m —1 in succession we get the m —1 * 
entries between f(a) and fía--«). The.values of Af(a), A*f(a). 

ete., аге obtained from the given series of values. Next taking 

К from ml to 2m—1, we get m—1 terms betweén fla + w) 

and f(a + 24), and so on. 

But this process is laborious for calculation of a series of 
values. It is better first to obtain the differences of the new 
sequence of values of the function, and then to'ealeulate them 
from those differences. 


Let A, be the difference operator for the new sequence 
of terms 


@ fla) (а), 1 (a+), 


then since (1+ дуа) fla +0) = (1 А)"у(а), 
the two operators Д and A, are connected by the relation 








, fla +)... 





e 








For practical purposes we can ignore higher differences 
beyond A*. 

Putting r=1, 2, 3, 4,....... we obtain the differences of the 
new sequence in terms of the differences of the old sequence. 
Thus 





x TED: DLS —2т) да 


| Сы, 
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1 iim ла т) 
war at + де Gam uma, 





149+ 8 
т 








When the differences are thus calculated, sequence (2) is easily 
obtained by simple addition., 
Suppose for example that m5 and differences up to A* 
аге retained. * 
Wo ensily find that 
A, = "2A —"08A? ¥ OBAS — 033044, 





Example. Annual premiums for policies maturing at a. fixed 
ago are given for the following ages :— 





, Age 20 25 зо 85 40 
* Prémium 305 — 345 40 4775 — 0025. 
Calculate the premiums for all ages from 20 to 40, + 
Forming the difference table for quinquennial entries, we have 
A at „Аз at 
20 90:5 
| a vem 
* 95 945 15 
55 5 
eo 7" 40 225 15 
T5 = 15 
85 4775 B75 
115 
40 50-25 
Hence il v v 
А }(@0) = 72 x 4— 708 x 1:6 7048 x 775—030 x ^75 
6908 n 
А г{(@0) = 04 x 1:5— 082 x -75 + 0256 x 75. 
m =0552 . Ы 
2527020) = 5008 x "75 — 0096 x "75 
: = 0012 
` 4, (20) = 0018 x -76= 0012. É 
j 4 = 
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We assume fourth differences to be constant. 
Hence completing the table of differences we have 





a a, A at 
20 305 
6908 
21 31-1908 -0552 
7460 — 0012 м 
22 319808 '0540 70012 
"8000 0 
28 32-7868 . 0540 "0022 
"B540 ‘0012 
24 83°5908 70562 ‘0012 
"9002 10024 
25 345000 70576 
“0668 
20 9574608. 
and so on. 


4.5. Method of checking the numerical accuracy of the 
‘Differences. 


The algebraic sum of any r consecutive values of А” is equal 
to thefast minua the first of the r+1 consecutive values of A71 
used in forming the r values of A". 


Lot the tabulated values for А”-1 aid А" be as follows 
4i: № Rs 
ats d ky 
Then ky=ha—hy, 
mh, hue 
Hence by addition we have 
Key elg sec Ea a tk =h, mh 
For illustration take the four values of Д, 
(Ex.) and the first and last values of 4, *; then 
— 70012 + 0 + "0012 + 70024 = 70576 — 0552 = 0024. 
Бичу | 0089-10040:« 0805 0002-0076 
= '9668— 6908 = 2760. 
4.6. Inversion ofa Series of Functions. It will be some- 
necessary to invert the series of tabulated values of a 




















me sens Ө 


CALCULUS OF FINITE DIFFERENCES 


to bring the last value into the first position, the 
into the second position, and so on. The effect of 
such inversion on the tabulated columns of differences will be 
that each column of differences will be inverted, like the column 
of functions itself, and that the odd order differences will change 
signs throughout, while the even order differences will remain 
unaltered in sign. 

4. Lagrange’s Interpolation Formula. Given т values of 
a function which are not cónsecutive and equidistant, to find the 
general expression for the function. 

Let the values of f(z) bo given for z=), za. 

Assuming f(z) to be a rational integral function of x of degree 
n —1, it сап be put in the form 
f(x) Ay (2-25) (2—23)......(2— 2.) + Age z) (7 >). 
+ Ag(e—2)) (2-29 
each of the n terms wanting one of the factors z—2,, 2—2 
z-—z, and cach is multiplied by ар unknown constant. 
aluate the A's, put z=z,, then 
fe) miim яа) (z1 —®э)- 











әк. 





ая) 





B 





n аа-а) 









sedem re) 





a Ww, ee 
or Ai Goad oe 








‘Similar expressions are obtained for Ag, As.. 
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Putting plays а) (а)ба), 
it can be put as - 


Hayne Lee, Ме, rci, 2 sum, 





Example 1. From tbe following data find the expression. for 
f(z) and also the value of /(10) x 
я 5 0, 75:9: 1 
Ke): 5, 7. AT, 480, 

Applying Lagrange’s formula, we have 
1—1) (2—8) (== -8- 
IAEN OSOA 1129027" 





z(r—1)(—7) 
*s8-1)0—7 
` MORET 5. 
Putting #=10, wo have /(10)=1195. 
"The value of /(10) is alvo obtained by putting #=10 in $4.70), 
Example 2. Two observations of а certain quantity wore 
made at the interval of в day from each other, the first gave 
for its value а and the second —5. When was its value zero? 
Here we are given /(0)=a and /(1)=—b, and it is required 
to find ж for which f(x) =0. 


„ 22-1) (2-8) 
*9(7=1) (73) | 

















“We have = 





Hence O=a(r—1) +be 





which is the time in fractions of a day from the first observation 
_ Ю the moment required. 
The result is also obtained from $4.2(2) 
Jr) m f00) 2400). 
EU nes are only two observations we take A?/(0) etes 4 
E 3 


TONS "Thrée observations of а quantity are made near 
of maximum or minimum. Determine the time of its 
1 mum or minimum. 4 
| 7—1220B. 
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Let a, b, c represent the билең of observation and f(x) the 
magnitude of the quantity at the time x, so that f(a), /(b) and 
Ale) are the observed values, Ву Lagrange's formulo we bave , 


foe ЕЕ a fia) + dU jas есче еў. 


As this function of x is to be a maximuin or mipimum, Ave 
equate to zero the differential coefficient of ffr) with respect 
toe, Thus we get > " 








= $t — c9)f(a) +(e? в?) f(b) + (a9 — U*)f(e) * 
HO — e)f(a) + (e — 0)/ b) + ве); 


the value of x at which f(r) isa maximum or minimum. Ву 
this formula the meridian altitude of the sun or of estar may 
be found when an observation precisely on шы arit cannot 


bo made 52. 


. 48. Irregularities in Differences. If, instead wf taking: 
the actual numerical values in the table, spproximafé values аге 
taken Ло certain given places of deciendls, the result will be that 
the differences will exhibit а amall degree of irregularity owing 
to the omission of decimals in the approximate vulues of the 
functions employed. No such irregularity will appear if асіп) 
values are taken, so that the differenees then will be perfectly 
smooth or regular. "Thus taking а function of the nth degree 
we know that A" will be constant, but when approximate values 
аге taken, the values of д" for different arguments will Wlightly 
vary. Thus n slight deviation from the true valve ofw tabular 
function will manifest itself by moans of irregularities in differ. 
ences. 


2 49. Accidental Errors. It hawever one уйше of m 
is in error by an appreciable qı 
~ nd its magnitude from an on of the table. 


> Let the valdas: ck Aha: бано Hie) be 
of za, а +1,...... and suppose that [f(a +@) is in error. 
quantity « во that f(o40)*« is tabulated in place of tl 
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value /(а + 6). Тһе ergóra n the table will be manifest-thus т 
а a* arte a as 
fa^ 
SE е V^ aeri RE 
sD- AY 
D PC к Aha > 
+ anatey 
Jass Af, ^em aun Ds 
пажа алш +з) arpata + 
ме» ° „ж.ш, asa ^Ш 
fases alte д, Bitte e A fiae) Me 
PME E Qa. Ае +3 аза эуел1, 
(e Mieten 
ML MEME Oat Se 
дењо asm TL Key ООВ 
` favor г aasa" ahon 
faim. 2% Ate) x atatt 
meni) Ae дааа enm 
$ А+) j 
tasin A oe 


From “an inspection of the table we seo that diflerendes 
lying along and inside the lines sloping upwards and downwurds 
from the incorrect tabulated value are also in error. To the 
regular values are alternately added and subtracted terms in « 
‘The law of progression and ipcrenwe in the coefficients of « 
along the successive orders of differences is that of binomial 
coefficients with alternate signs, We have only to carry on 
differencing to that order ot which the differences of the 

- correct functions would vanish or sensibly so. The location and 
magnitude of the error will then be clearly shewn by a succession 
о! positive and negative terms with binomial coefficients. 

‘Thus if the values of A‘ vanish in the correct table, the 
fourth differences in the incorrect table will be 0, +«,—d4e, + Ow, 
+4, +e, 0. The initial term +é is the error sought. The 
/'imeorrect funcgion ia them obtained by tracing backward along 

- the lins-of downward slope. Thus in the table we come to 
SL f( +B) +e: ‘The correction is the negative of the error «. Hence 
E value is f(a + 6). А 


















| 
| 100 CALCULUS OF FINITE DIFFERENCES 
| 


Example 1, Find the error in tbe following table of 





| i=. x 
| = fe МОЕ А A* с Atte 
1 1 
7 
UM s. 12 
19 6 
| з s 18 0 0 
aT 6 
Н 4 ы м. -10 «10 o 
o1 -4 
5 125 20 +4 -40 0 
в +36 
| 6 206 56 +0 o 
J 137 -u 
7 зз 82 -340 o 
100 ۴ 16 
B 52 48 -10 wO о 
в 211 6 
“9 7m ы 
22 d». 


© 


1000 
“Hore the true value of A* js zoro but in the table wo have 
а complete alternation of signs and iti» easily seen that the 
tabulated values for A* are the binomial coefficients of the 
fourth order multiplied by 10. Further ns the first incorrect 
entry is — 10, the correction ia + 10, and tracing back downwards 
from —10 we come to the required incorrect entry 206. Hence 
‘the true value is 210. The same is also obtained by tracing 
back upwards from the last incorrect entry—10 in the column 
fora: 
















a 

In the columo C are given corrections to A* due to correction | 

0 to /(6) and the last column gives true values of At, 

2 410. When two or more consecutive and nei 

n are in error, the detection becomes difficult and com- 

and in some cases even impossible. ‘The entries да 
tain systematic errors which do not affect the 
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Example. Correct the errors in the following table of 
Қа) 2? —5a +12. E 











= 1/9 Aa д» ‚ А*+С{ бу “A240, +0 
1 в 
m 12 
e 6 è 
3 2 18 
14 в. 6 6 
4 36 за 
04 -18 4 6 6 
5 120 n 
в > |-и "D 
в 108 м 
из æ | в -20 6 
зо 02 
114 -u 30 в 
8 484 38 
212 16 | 16 -10 в 
006 54 | 
206 
10 02 | = 


Here the third difference i» couvtant (=6) and the first 
incorrect entry is 14. Hence in the column C, we put the 
correction —8, and following the law of binomial coefficienta 
of the third order with alternate signs, the corrections are 
—8, 94, —94, 8. The column A%+C, gives the corrected 
series for the error 8, but here also the first incorrect entry ia 
—4, which suggests a further correction of 10. The column 
Cg gives corrections for the error —10. The pal corrected 
column A3 +0, + Су gives the truc value 6 of the third differ- 
ences. The incorrect entries are traced to be /(5)=120 and 
{()=810, which should bé ((5)=120-8=112 and f(7)=810+ 10 
= 320. 

“BAA. If the error begins Ín the column of earlier differences, 
higher differences of some order will indicate corrections 
correspcnding ‘to the binomial coefficients of в lower order 
than that of the differences in question, Thus if ДЭ requires a 
correction of the order 1, 3, 8, 1, it is clear that an error 














exists in A7, of which Д? in tho third difference. In general, E 








© 
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"when A* requires a correction of the rth order, an error will exist 
in Дт”. 
„a y When cither the frst оғ the lost function of a series is in 
"reor, only the first or the last term jn each order of differences 
will be in error by the same ambunt. In such cases, methods 
explained above are of little assistance in correcting the errors. 
4.12. Mechanical Quadrature. One of the chief applications 
of interpolation formulm is the evaluation of definite integrals, 
e., the determination of the areas of curves bounded by the two 
extreme ordinates, the -axis and the curve, When the equation 
of the curve is known and integrable, the result is obtained by 
direst integration. But in physics! and natural: phenomena 
the above conditions aro herdly satisfied. "Thus astronomical 
` observations, retes of mortality, ete.. are obtained from actual 
observations and the expression for the function corresponding 
to these data is not known. The process by which such arcus 
are calculated without a knowledge of the form of function is 
known às Numerical Integration or Mechanical Quadrature. 
„ Suppose n +1 equidistant ordinates /(0), f(1), /(2),. „ 1) are 


eed: “ ia required to find an approximate value ot SK Jede. 
? p 1 d dee Newton Я of interpolation, | e 
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neglected: Тһе expression. (1) giwés an approxunate valua of. 
the integral Е A 
4.13. Particular Cases. , М. 
(1) Simpson's Rule. Suppose that iria entries f(0), | abi 
102) are only given. Them. retaining term: far аз the second. 
difference we have 





$ "ley =) + 24/00) + A 23/0). т, 
Substituting the values of A and А? we have "m s 
коме oaa + о. 


M the third differences, of 70) aro constant, Simpson. nie 
gives tho eact value. 

.. (2) The Threc-cighths Rule. ` Ii four зарана orfinates аге. 
given. wo have 








* рома St +ай жайа jy: euh. 


(8) Weddle's Rule.* Suppose that soven equidistant ordinatos. ` 
are given. Then 


feuz 182/(0) + 27/0) - 24A ^f(0) + A auo) » 


EON anf + ARON 
ARSE the A's in Глена of the CLER ‘it will bo 
found that ratos of the terms aro largo. “This is due to. 
‘the awkward eqe Т, TZ o 9t Af дее" Math ا‎ d 
Ws ve. Ё : ^ 
are small, the. ror involyed in replacing “qq T 


E be Regligibley but the formula. za be much simplified. 
‘The new formula Will involve ân error of = agam and will be 


strictly correc up to the fifth differences m Абау correct, °. 
the D y On redyotion it becomes 








42 
Y do" ien by 5 
















ا 


Im T. 2 


9 
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4.14. Change of Unit and of Origin. The formule of $4.12 
and $4.13 have been obtained by taking the interval to be unity. 
At however the interval be w, tho given values of f(a) will be 
КО), fied. f(2w),...... Hence taking the new variable to bo =’ 
we bave 





a'-er and г. deed. 
"he new limita of integration will be О and мө. 
Thos Simpson" 





pula will transform into 


L T enam [no eas t 

or chinging the variable to x we bave 

- Ede 

aq NL [ran нез) 2 
Similarly Weddle's rule for the interval » will be 

Xr eae у, [rto i tao) + оь) + 40) + бы) + fom |. ° 
Any formula of approximate integration for unit interval can 

thus be transformed to one for interval w. 

А change of origin from ==0 to к= а is obtained by substi: 


tubing f(o). Даны), f(4 + 2w)....... in place of f(0), flu), f(2u), ete 
in the formula. bus in Simpson's formula, it the orig ie i n 








ТЫ» formula wap obtained by James Gregory * in 1008 und. 


; ся à 
n С If the origin is transferred (оли v and the same aren as that - 


given by (1) is required, the transformed formula will be 
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Similarly-the Three-eighths rule becomes 
= °" зах = 5 [ros З/(а + w) + Sf(a + 24) + fla +3) Я 
4.15. Extension of Simpson's Rule: ' 
We have бе {ro + p 


Changing the origin 


Sout i [rne e]. 





Similarly. Ге [феде 


^ í idm 1 | [rem ann ят]. 
Henco by addition 








IN í O) + аа) + 2/2) + 4/08) +34) +... 





This formula gives better result than the formule for the 
interval п, - 


fr rae {I+ af) +120}, 


and is extensively used in practical mathomat 
If tho origin is transferred to а and the interval be w, wo get 
11) in tbe form. - 


2) т (одан у [№ КЫЛЫ дозы T 





чае E f fiae) |. 


"Example. Kind the value of Seca cos 6 40 from the following 
гама: 


coso в cos ð 
0-939093 36° 0800017 
0-913545 40* 0700044 
0-882948 44° 07719840 
0*848048 А - 


чш}. 
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| (2) Taking only three values for 20°, 82° and 44* and applying 
| Simpson's rule, «e Бате 
at L 

so 





cos 0 d=" [eon + soon 32° + con | - 


= {0-ss0009 «o co] T 1ں‎ 





А + e: 
5 < 0051225. 009810 х 5051225 
85904117. bi. 

| . (2) Generalised Simpson's rule gives 











E con вв 1 [eoram eem 24* + 2000328" + 4 соко" 
‹ Caen +4 con 40" + con 44>} 
І р 
' COSS 
= swn s Р 
ы =-852638180...... . ne. 
E (8) The Tbrec-eighths rule with interval of diffe в" 
== ў : Ж 
= coed 8 9 [oon 30° + 3 con 28 cos 30° + con dd ] 
m =3° x 6/724928 = 352640.. E- 
" 
(0 Weddle rule gives 3 


5 d ave 20° + 5 сов 24° + оов 28°+6 con 82° _ 





+ cos 30° + oot 40° + cond} 





=4° x8 10:887081 - 
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Thus the values obtained by Generalised Simpson's rule and 
Weddle's mule are correct up to 6 plac 
given by Simpson's rule and Three-cighths rule are in slight error 
in the Sth or 0th place of decimals, being slightly in excess of 
the true value. 

5.16. О. Е. Hardy's Formule. Suppose the fifth differences of 
f(z) are constant, and values of f(x) for x0, 1, 8, 6, б are given. 
fla) can evidently. be represented „by а polynomial of degree 5, 
Put 





Ma) = ао az 0923 аж? жа a aga 
Changing the origin to —3, 


fuente becomes §” rine 








with the values of f(x) for z= ~0, —2, 0,2, 3. 
Now 10) & ao, [(—2) + fi2) = 2a, + Bag + 3204, 
К-8) + 3) = 20, * 18a, + 16204, 
and s j 
п) Vrbis m bag + 1805 +07 


Now let us seo if (1) can be expressed in the form 
E] MIO) + wf 2) (00) + v{/( 8) + /(8)}. 
Equatíng tbe caeffiéienta in (1) and (2), we bave 
A222, 11702 and у= 728, 





зо that 
f? доага 2/0) + 1'02(/(— 2) + [(2) + B=) + f(9)), 


8) or. fam 280) +100) +1 62/0) К} + 22/09, 


which is Hardy's Formula* (37). 
И the interval of differencing bo w, thia becomes 


Gy Hapa e 2840) + бы) + 176209) + 69) 2 2/88]. 
o 


om o 
пз» r s 
i Ay Е адаа ma 2889-29) 1 atre efte) * 39/099]. 
and soon. 


ао. King- Lie Contingencies 
Б, A 5.3 x SH 
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Hence vi ia 
(4) шем 28{/(0) + 2/(бь) + 2/(18%) + ...} 
' i = + 102 (f(9) + (бе) + J(u) +... 


$224 /(Bw) + (Ow)  [(156) +...) ]. 
"his is Hardy's Formula (38) 
In problems of mortality table, w is во chosen that 7u is near 
about 100 ; (4) is then put iñ the modified form 


- fein oC 28/0) +1°62/(ь) + 2^2/ (Bu) + 162/50) 


} Y Bf (6o) + 1762/(79)] -. 
р ‘This is Hardy's Formula (30а). 





Example, Evaluate g — 
© viz 
А a = i 
ч Here fla) С. Teens 5 


Dias Er gg 28 21547008 (1:02 x 21080285 
- oyi- 12 
+2 2  х1:-0затобву=0 "5235085. 








slight error in the seventh place E decimals. E 
ы 447. Laplace's Formula for Mechanical Quadrature.* 2 
Since Ee cP-1«A л D= log (1+ A). 
(Hence a= PTET" б . \ 


3 EES ES |. 





‘The true value of the integral DH or 075235087, » there ін в 


Й 
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we get, on taking the limits 1 and 0 and putting f(z) = Aple), 


уон 15 + Ts prm 





-/0)«/0) AO азо) _ 
2 12 24 








ilarly 
а ауа = К) + AT fü) , азу) _ 
боа T" u2 ! 3 


and so on. 
Hence on addition of н such expressions we get 





w ffol ep) fine iso fm P 





wi 4 [ал-ал 


+h [ar -A fl 
for А {о + A" fU) а eA fin Y] 
ACT ARD) + А/@) *...... + А/ї#—1)} 
=A Mo) - f00)- * 
Formula (1) labours under tho disadvantage of containing 
differences of f(n), which cannot be calculated from the values 
of fO), f(D, css, fi). а t 
А slightly modified result can be obtained tn 











j See EA а M ON po. 
ЕТ for 4$ (2) and removing equal terme we have 
с MET fe n em 
ES eR- 
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Substituting the above in (1) we get 
2) fte = 1 era) +102)+.. 





-ya 
ftn d 


- 3 fare -n-aon 





em fone -» кэч} 





Неге ай the differen 
This i 


41 


can be aalculated from tho data given. 

simple form of Gregory's formula £4.19(2). 
Euler-Maclaurin Formula." We have from їз. 2(4) 

знае + (fas — fe, E +... 

‘Taking this between the limits 0 and n, we get 

Ella) = /(0) + /(1) + (2)... fi 71) 


= pene} фидо + бемо} 





m рн) mol Ries 


Henee 


емет р) + rays у ши ) ID 


-4 ipm =O} «d des то} — 


The above is « sünple form of Euler-Maclaurin expansion > 
for unit interval. 


Changing the origin to а and taking « as the interval of ж 
k we obtain a more general formula : 


@) ip fr)dz = ito + fce ы) + fla + 20) +... 








C tetni) + farnu- pos =P 
| яо" = an 2 freno ref du 


* Eules— Com, Acad. Sci, Tap, Рейтер. M eI T. 
Maelaurin—T vestire of Plustons (742), p. 672. > 
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Него the terms involving derivatives and the integral are 
affected with w (see 14.14). This i the basic formula for 
expressing a definite integral in terms’ of given ordinates. If 
the form of the fonction js Кооп and differentiable and the 
integral can be evaluated by some method, the above formule 
бап be used ín finding out the sum of a given вез When 
the ordinates are given but the form of the function is not 
known, the derivatives сап be expressed in terms of differences 
(see $4.19) and the integral cun thon be obtained, 

If fl») be а polynomial, the Maclaurin series terminates 
and ig exact. 











Erample 1 Evaluate f? "е correct to five places of decimals, 


In Kuler-Maelaurin formula /4.18(3), put a=1, n=5, w=02 
and ft) = 1. Then 


айт 





i gig ]-- 


10 + 755506 + "25000. 
= 01250 + 0000625 — „, 





N? 42 = 008148 6095, 
ax 


agreeing with the true value of log,2 correct to five places of 
decimals. 
Example 2, Calculate correctly to ten places of deci 
the sum of E 
109? * toa" * 1049 * ^ 2002 
„Find also the sum of the series to MG 
„< неё 4100, v=2, n—50 and f{(x)= 
^ f Maclaurin expansion we have 


zh 

1 (200 4=_ 1 1 1 41 

100 339 271009 * 1083 * 7 * 9-30 
1 





als 





Hence from Euler- 











AEA INE ЖҮРЕ ТҮ у 
ни} si 20057 * ios] G 











spp эшк e 
з [100% 2005| ~ 15 (1005 200 i 
= "0025 + 0000025 + -00000029106006 — 0000000000258 
= "00250279164083 
which is the «um required. 
Again to find the sum to infinity, we have 
moi + тет * im 8 ad infinitum 
=: ee ИКЕ O E ЖАУ 
2100 ° 21002" 37100* 15 1005 
= 1005 + ‘00005 + 0000003333333 — 0000000000960 
0050508332007 
which is the aum required. 
419. Gregory's Formula of Numerical Integration. Here 
tho terms involving derivatives in Koler-Maelaurin expansion 
are replaced by differences. 


Sin D „1+4, 
(1) heme „рее {а-а = 
Thus — ef(o)- Afío) —AA*f(a)  4A*f(a)— 184 








wD=log (1+ А), 











ма + no)  A?f(a + 
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wf" (a +nw) = A*f(a-- n it лая de) _ 


and so on. 
Substituting these values of the derivates in Euler-Maelaurin 


expansion we obtain 





1 


(2) 


"This i» Gregory's Formula, discovered by him in 1670. 


° edem Mia) fa +o) +. 


1 п 
ü [aree 


кл. В {ал+ 8759) -A*f)] mo 


60480 


+ бажн 


* jA^[(a +R 5u) +... 








›-эле)} 


- fare +N Tu) жатка i 
و‎ fom +R) Ala) k 


-ii y fate жи4) + A*f(2) } 


35) + A?[(a + п — 8) + HA *f(a — Su) + 
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1o) + Mía € ne) 


те 


does not require data beyond f(a + ны) to form the differences. 
Example. Obtain the value of i сов 6 40 from tho data 


given in Ex. on p.105. The difference tabl 


o сов б A 
20" 0939603 м 

— 7026148 

24° 0913545 
— 030507 

28° 0882948 
— 84000 

32° 0848048 
* Ж — 89081 

36° 0809017 
: — 42078 

40" 0766044 
^ — 46704 

44° — 0719940 

| 8—1226B. 


— "004449 
— 4808 
- 4181 
— 92 
— зз 





17 
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Here а=30°, w=4° 





5' 
Непое substituting these values in $4*19(2), we have 


=} 





"сов 0 dé "cos 20° + соз 24° +... + eos 44° = +5'878635 


d (con 20° + cos 44°) —0°820516.5 

- 2 (A cos 40*—A cos 20°) + зо 

Д (A? con 36° +A" cos 20) + 8408 

— 19 (дэ oon 32°-A? con 20") — 17 
тю 

-ipp (Af con 28° +A con 20°) — 0.9 


ZF 5051110 





са Ge cox 6 dêm ™ x 5051170  -00081817 x 5051170 
20" 45 
= 252038 
which agrees with the true value given on р. 106. 


4,90. Lubbock's Formula of Summation. Suppose that cach 
interval in Eulor-Maclaurin formula is subdivided into m parts, 
then the new formuls is 


m E Gi rade fp ef) fot D) fom) 1 


- фе) ж/а tne) 1 -= freno ral 





"isa a +m) - f) i 


Eliminating the integral between (1) and the original formula 
$4.18(2) we get 
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в Hla) + a+) + fla +24) +... + Дажа) 


=" (feste fna] meL [rema] 





[re +n) — 1a) | 





тә 
зу [enel Га 





This formula is due to Woolhouse. 





values іа terms of differ- 


+ fla +n) 





0) e+ e+ )+/(а+%= }+ 


- m reef + ss. ены} - لے‎ 





fra + n) +f) 
= [arr + »-afo)] 
С 90-1) лэди + өы дәш} + 

Formula (3) involves advancing differences of /(а+ты). 


If backward differences are employed we geta slightly modified 
result due to De Morgan, t 


w 9+/(а+ m )+1( 04%) + 


т? -1 


{are + ны) – Дуба) ] * 





+ fla+nu) 





=m [ron 





faena] 





= fra + ne) + fla) } 





T [ает] 





m (atte ato] 
б 





Lubbock—Camb, Phil, Trans., 3 р. ma. 
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3F {дэле +n Boj дэле 





{ Aja +» =4ш) + Atja) } =. 


‘These formule enable us to obtain the sum of a large number 
of equidistant terms froma smaller number of terms taken at 
intervala m times greater. 

Еготріе |. Find the sum of 

105 +115 +129 + ...... +505. 

Hore 4-10 and f(x)=2*. Lot w=10 and m=10, then 

forming the difference table, we have 


* a a a At 
103= 1000 
7000 
20= 8000 12000 
19000 0000 
2032 27000 18000 о 
87000 6000 
403 = 64000 24000 
61000 
505 = 125000 


Hence from Lubbook's formula $4.20(4)* 


10 {упо +120) + во) +140) + reo} -$ freos шо} 





-3 {arao—arao} = farao ao) 


90x1890 ( 

- ° [^ 12o- 2.24010} 
= 2250000 — 567000 — 44550 — 14850 
= 1623600. 

ExaMPLES. 


1. Find a rational and integral function Sf = of as low а 
degree as possible that shall assume the values З, 12, 15 and 
—91 when z ie equal to 3, 2, 1 and —1 respectively. 

i M [C. U., М.А. & M,Se., 1936] 
КЫ а < М 
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2. Given the table of values 
2 -8 -2 -r o 1 
y 16 7 4 .1 -8 
Find by means of Gregory-Newton formula an expression for 
y as а function of =. 
[C, U., М.А. & M.Bc., 1934] 
44 an expression for f(x) from the following table : 
4" po ^a А 
Ju. 8; X, 08 es 
[C. U., М.А. & М.Ве., 1039] 
Find an expression for f(z) from the following table ; 
= ФА ET ЫА 5 
f) -6 о о о `зм 18%4 
Hence obtain /(4) to the nearest whole number. 
5, Find an expression for f(x) from the following table: 
и” 5 ГУ 8 7 
fz) 5 205 57 497 эз 
6. Complete the series for f(z) from z=0 to ж=12 from the 
following data: 
= о 4 8 12 
Кг) 271 2888 3087 9175 
7. Find the value of /(60) from the given table: 
= 40 48 56 ва 
fe) 1497 1581 1772 10596 
8. Correct the errors in the following table by the method 
of differences : 


8. Fi 





1 
= 1= 7 z fe 
on 4702 81 5:226 
25 4348 "3з 3080 
-25 4:000 "95 2857 
"e 3704 "aT 2708 
“20 3405 EJ 27504 
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9. The following table giving the values of the moon's 
altitude at equal intervals of time, contains an error in one and 








only one entr ^ 
300" 807 5876 30043" w 2775 
345 32 4453 9 41 586 
t 351 30 30 `9 15 ы 6'3 
857 30 a 8 22 11 30%5 


Utilize the theory of finite differences to detect the position 
of the error and correot it. 





[C. U., М.А. 4 M.Sc., 1035] 
10. Correct the error in the following ephemeris of Moon's 
Intitude : 
Date Moon's lat. 

Mayes — -1* 59 52 спо 
90 -1 9 «4 2 15 
95  -0 4 27:0 120 
100 -0 6 45:8 125 





105 +0 32 39 9 
11. Correct the error in the following table of Mean refraction 
corresponding to the apparent altitude of tbe star: 


alt, mr. 

10° 192 20 2 388 
12 4 mc 22 2 983 
14 з 405 24 2 102 
16 3 1874 26 1 5879 


18° y ots 
12. Correct tho errors in the following tabulation ot 
Қа) = 203 — 252-40: 
т AR LOE 39. & $5 mig eu 
(=) —19. ~6, —17, —40, —63, — 79, —61, —4, 85, 242, 471, 784. 
18. Compute the value log sin 9*22' from the following table 5 


log sin 6 ө log sin ð 
10° . 92800 
ve. in 949800. А — 


9-019023. 
908589 
. 914856. 
9719425 
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14. Calculate the value of sin 33° 13' 30” from the following 


30° 81° sae” gar 84° 
sin 6 5 615 — 5990 зб 25600 
16. Obtain the annual premiums for ages 81 to 39 for an 
assurance of £100 from the following table : 


Age ё, а Age Е за 
20 119 6 40 8 410 
25 2 зи ‚ 45 815 B 
30 295 ю 410 7 
a5 216 2 





10. The following annual premiums for policies maturing at 
quinquennial ages being given, find the premiums for policies 
maturing at all ages from 45 to 50: 

Age 45 50 55 ° 00 65 
Premium 27871 2404 2089 1862 1712 
17. If f(z) be of tho form а + bz- cz*, shew that 


INS [so 8r) -@}. 


18. Show that frg =i {m+ 22/0) ım} 


19. Obtain the approximate formula 
LOE I {f+ (78) 
A 1 
shewing up to what order of differences it holds, 
20. Prove the approximate formula 
ffe fidem ars) +14) +10) + б), 
and show that the formula involves a small second difference 


error. 
21, If fl2)—a+be+cx*, shew that 
fife 2100+ j [o-a] 
Ex Hence find an diese EN of 


С i z а». 
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E Prove that. 
remo) (3) | dor 2)-2(G] 
А 
Hence find the value of ү уе, having given 


10) = 105, f(1) 212, f(2) = 421, /(8) 5 749 and f(4) = 1050. 
23. Establish Shovelton’s integration formula 


{еа тав [809 + fa0)) + 854/0) + fi) + J(7) + 10} 





+1542) + f) + б) + (8)  86](5)].. 


24. Find the values of одаг from the ollowing tables: 


зо 1 2 з 4 5 в 
© fü) 140 "161 — -170 -190 ‘204 7217 7280 
(ii) f(&) 1 “96154 -02456 88900 :85480 '82108 70091 
25. Find by Weddle's formula the values of 





$ ds 
Boise ө. (=. 
26. Find by (1) Weddle's and (2) Hordy'e formule the 


value of 
p de 
ovis 


27. Find the value of (7^ cos @ 40 from the table at p, 105. 
28, Use an approximate formula to evaluate 
безен 


кт dis 
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80. Find by Euler-Maclaurin formula the values of 


1 1 1 
O aoi» зо * обуз 








1 1 1 
) Tor 1088 * 105% * 
91. Calculate by (i) Lubbock'e formula, (ii) Éuler-Maelaurin 
formula. 


... ad infinitum. 





a o goa 
de 
82, Calculate by (i) Euler-Maclaurin formula, (ii) Gregory's 
formula 


7 зо а= (103 de 
s Дуу 
38. Prove that Grogory’s formula may be obtained by 


making m increase indefinitely in Lubbock's formula. 








34. Calculate {уя о Ца from the following table: 





= 1'05 106 107 108 го 110 

fiz) 125886 126090 1-28619 1:0254 181908 188505. 

85. А series consists of nine terms. The sum of the first 
and the last term is 78°84 вой the sums of other pairs of terms 
equidistant from the first and the last are respectively 73:07, 
7352 and 72°19. Find tbe middle term. 

36. The number of deaths at age 44 is calculated from a 
mortality table giving the number of deaths at quinquennial 
эе points 32, 37, 42, 47 and 52, and is estimated at 1000. 
Calculate what must bave been the tabular number of deaths 
at age 32, it being given that the number of deaths at ages 37, 
42, 47 and 52 were respectively 869, 952, 1112 and 1971. 

87, Prove hat Weddle's rule for approximate integration 
сап be obtained by the combination of a -ordinate Simpson's 
formula with а 7-ordinate Three-eighths formula in the ratio 
of 9:—4. 











© 
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88. The numbers of members in a friendly society were 

available for the following years : 

1922 1028 1924 1925 1008 

995 ` 998 юз 906 976. 
‘Pho numbers for the years 1926 and 1027 were calculated on the 
assumption of a constant fourtb difference, Subsequently the 
actual number for the year 1026 was discovered to be 1002, and 
n fresh estimate for 1927 was made. Find the original estimates 
for the years 1920 and 1027 and the revised one for 1927. 

39. The following table furnishes, at noon on the given 
dutes, the value of log D, where D is the distance of Venus 
from the Earth measured in radii of the earth's orbit: 

log D 

July 1, 1912 732 
Em 
72904358 
2304402 
2894204 
72598860 
Determine as accurately as you ean from the given data— 
(0 the value of log D nt noon on July 4; 
the time at which log D was greatest, ond the correspond- 
ing value of log D. 








EINE 





(С. U., M.A. & М.8е., 1934] 
40. From the following table of expectation of life, estimate 
By 





= 20 25 зо 85 M 
Е, 44-21 4000 8581 31-71. " 
If in addition the values for ages 10(69°08) and 15 (48-57) are — — 
ven, obtain a revised estimate for Es. 3 
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DIVIDED DIFFERENCES 


5.1. Divided Differences. In previous chapters we have 
taken the argument z to increase by some constant quantity 
во that the difference between any two successive values remain 
the same. Here we sball consider the case when the argument 
increases by unequal increments, so that the difference between 
two successive values are different. The values of the argument 
аге here unequally spaced. We introduce here another class 
of differences which goes by the name of Divided Differences, * 

‘Take a function /(z) whose values are given for z=zọ, zi, 
aser. The intervals z} 20, xg j, ту—ху, eto. need 
not be equal. The divided difference of f(z) for the arguments 
fo, X. denoted by f(zoz,), or by the symbol [жогу], and № 
defined by the relation 




















(29) а) f) (9) =. 
w Jeor) D) Le) = Е) еа) fau, 


Similarly. 





These divided differences with two arguments are called 
Divided Differences of the firat order, 

Two divided differences of the first order having a common 
argument defines в Divided Difference of the second order involving 
three arguments. Thus 


o оу) = [nons za] کا‎ fetal, 





Similarly a Divided Difference of the third order involving four 
arguments is defined by 


@) онна) m [orina] = оаа) балала), 


‘This term was first окей by De Morgan, Dif. and Int. Cale, (149), p. 060, 
and afterwards by Oppermann, Jeevn. Inst. det., 18 0889), p. 146. 
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In general, a Divided Difference of the nth order involves 
(n+1) arguments арф із defined by 





ч) fzozyza 24) [готаз] 
We ean exhibit these results in a table of divided differences 
as follows: 
= fero) . 
У Г) 
= foy) f(zozyra) 
аля») Mor zaza) 
= Iles) Гауга) Hox Fats). 
ПЕ) гугл ут) 
=, Ils) (кагу) 
Nasza). 
= Nea) 


Example. Calculate tbe divided differences of f(x) from 


the following table : 
el. % 6 17, 19 


f(z) 22, 30, 82, 106, 208, 
[0. U. M.A. 4 M.Sc. 1938] 
‘The divided differences are exhibited in the following table : 





= fin) 

1 22 
8 

2 зо 

4 82 
в 

7 108 
20 

12 206 


5.2. Relation between Divided Differences and Ordinary 
Differences. Taking tbe interval of tbe argument z to be 
constant, equal to ш, we have 

К) — fto), fzo + ө) für) = à 
س‎ fev) exi feo Е 470) 


1, we get 





If in (1) we put 
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Similarly йвз) = Prom) =I 
А/М. 


žo—žg 





= тазе. 





When w=1, we have 


котла) = А60). 
Aguin eos sas) = оза) — жугуш») 
о-в 
Ато) – АЗ) 


20-2, 


азе), 


2- ы 


7g A*f(zo), when w=1, 
In general, assumo 
(2) Moy zo + ө, 2o + Duy... tnu) = 3a "feo. 


then fizo, 20+, хо +20, 
ео, totu 


+I 


ео, zo +... mg tnu to) 
=o tnet u) 


N O E La oro 
© =A (ачаарда Mle: 


Now (2) holds for n—1, 2, З and hence from (3) it will hold 
for n=4, 5, 6,...... Thus (2) holds generally. 

Again putting w=1, = for zo, we bave 
a) Ke, 2+1, #+2,...,=+т)= A*f(). 


Thus divided" differences аге generalisations of ordinary 
differences with constant increment of the argument. 

A direct proof of (2) depends on the propertiea of $5.8, 
q.v. 








"EXT 
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53 Divided Differences are Symmetric Functions of the 
Arguments. ^ 
We have by defnition 
= dzo) , fini) 
Ment Ыы сү 
Again 


Иво” ула) = 








ES M +Де® | 


ЕЕГ 
Ка) + 


fte flea) 
MEDIE СЕИ 

То general, it can be easily proved that a divided difference 
of the nth order 






„= 











=. cr to) E 
et a)" рот) toa) . 
* __ fee) жык } 
(n ао) ra alin = 20) (=. To) ire а). ei) 





Tt is easily soon that the value of the divided difference is 
unaltered by an interchange of the arguments in any order, and 
thus the divided difference is а symmetric function of ite 
argumenta. iy 

© When the interval is constant, equal to w, we have 


О ffo otw, zo + 2w., 20400) 3 
foa Fro 49—24) + LTE 
EU eT ee ice "caras 5 De = 


4 - d 
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Yos, егт) а= 1e0 Pzd -oz + Ho) а) 
Л 
eei * TA 
and similarly for differences of higher order. 
(i) A divided difference of e.f(z), where © is constant, is 
© times the corresponding divided difference of f(e). 


For example e d= ofits) ae Шота). 
zo zy 
бй) The divided differences ‘of the nth order of к“ are 
constant. Б 
Неге Siz) 








7 Коту) 





+2129 ж... +=1-1, 
LES 


a homogeneous function of zo, ¥, of degree n—1, Again 








- all s a 
Tomita) ау i ma) ET 


SE >й 
21-2 















*+...+#3-°) +... 
which is a homogeneous function of zy, 21, z of degree n —2, 
‘Thus the operation of taking the divided difference lowers 
* tho dogreo by unity. 
1n general 


Neorize--=/)= S. 


indicating that the factor z,—2, has been omited. Ир 
вее that the divided difference of f(zo2,..-z.) is а constant, 
“The divided differences of order (м +1) of z" are zero. 
(fo) The divided differences of order n of а polynomial of the 
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5.5. Divided Difference expressed by Determinante. 
Wo bave 














а) E, 





-2 
Е 






ro difference product of ғо, 

the + indicating that z, has been omitted from the arguments. 
Now, by Cauchy's generalisation of Vandermonde's theorem, 
we have tho difference product of zy, луж. а. Пу>йе, а) 





= 3 fle.) x diferenco product of 9,2), 





Seo) fe) Даа) += пе. ‚| 


а 1 у Ij 
Eo эү Xy 


+| 28 t on 






a 


Thus divided difference of order n i» expressed as the quotient. 
of two determinants each of order (n + 1), 


eee mms 


в. 
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8.6. Newton's Interpolation Formula with Divided Differ- 
ences. Writing z for хо in §5.1(4) we bayo 


Naxx... 









Hence 
(к= уг. 
Similarly 


diy) f(z ma) Un — mg) ауз в), 
Fix) fé.) + Gr miferi). 


Therefore. 
He) = fi) (®— n frm) 4 
= (ку) + (x = )/(=уга) + (2—21) (2-х) ати) 
= flay) + (е) ng) (2—21) (2— 2) fia 795) 
+ (2-24) (023) (2— yf inn mama). 
Proceeding in this way we ultimately have 
a) Ка) = fey) + (®— fir zg) (2 — 24) (2 — rg)f(z yz 925) +... 
+ ez) (2— 23)... (0— =, - aM rng mo) 
+ (0—21) (2-2) 
+ (2-2) (2—20)... 
‘This can be put os 







(2) Дае) + Xe) (еа) (e mfra sns) Rale), 


where 
(8) Rae) = (x—21) (6 — 2)... (e — te) fr iru. e). 

This ie Newton's general interpolation formula * with the 
remainder term F(z), and no restriction is imposed on the form 
of the function f(z). This can be considered as the Balo Formula 
in Finite Differences. 

‘When the values of f(z) for ==¥,,2g,...,« are known, the 
evaluation of {(2) depends on the value of Ra(z). If Н.) 
be known or negligible, the value of f(x) is obtained from (2). 


© Newton—Principia (1687), Book 10, Lemma v. 
912268. 
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Tf f(z) be a polynomial of degree n —1, the divided difference 
of order (n —1) is constant and of order n 
Herz... 








so that in this case ^ 
and ء‎ 


O fia) = fle) + E e-a) (є—гь)...(е-=„)/(вула...е,. у). 


It (в) be not n polynomial, we seo that 
R.(x)mO, for w= 24, жа, о, 
mo that the right-hand side of (4) gives а polynomial of 
degree (»—1) as the equivalent of f(z) for the given set of 
arguments, 

A polynomial, of degree n—1 at most, whose values at 
#1, Fass. are the same ав the values of f(z) at these pointe, 
is called an Interpolation Polynomial of f(z) and is denoted by 
I,-3(2). To obtain the interpolation formula. (2) we have to 
add R,(z) toit, Thus 

f 9 1, -.(2) + В„(«) 
where Ia (=) stands for f(z) in (4). 

‘An interpolation polynomial is fixed by the values of the 
function at the given points. It does not depend on the order 
in which those points are considered. 

Ezample1, Find the interpolation polynomial from the 
following tabl 

















w 09153 5 
f) 2 8 12 147. 
‘The difference table is exhibited belos 





= fæ 
= 0 TE 
1 8 4 
1 
2 12 9 
5 147 5 
Hance M 


d) 10) + =f, 1) += —1)/ (0, 1. 3)«26—1)0-2)/0,1,5; D) 
=2 + 2+ 422-1) +2(2—-1) (2-2) 
9-24 
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Example 2. Find the expression for f(z) from thé following 
data: - 
= i ^9 4 5 „То Ру 
fe)  -- о 10 о -mo GT 


Since /(2) = О for z=2, 5 and 9, it must be of the form 


f(z) = (2—2) (2—6) («—9)o(2), 
whore ¢(2) is в polynomial of the second degree in ж. 
=. = 1.4.890) = 8300). or $0)=8 — , 


14)=  2.—1.—59(4)-109(4), dt $(4)—12 got 


Ahe 6. 2. -29(7)= —209(7), or 9(7)=80 
Hence (x) =8 + (e—1)8+ (2—1) (2-4) =z 22-4. 
дю = (@—2) (®— 5) (к—9) (82 — 92 +4), 


6.7. The Remainder Term, In $5.6(2) put 
Р.а) = а) +" 











еа) (8—24)...(e— 2, less. ол), 


Then Raz) = (а) — P«-1(2). 
Evidently P, (2) is a polynomial of (n —1)6h degree and ite 


(n —1)th derivative is 
Р! (а) = (n —1) Иа). 

Let us now suppose that the function /(0) of the real variable t 
and ite first m—1 derivatives are finite and continuous, and 
100) existe, in the interval (o, b) bounded by the greatest and 
least of the arguments. 

Wi wis Mass ae 

Then since R(t) vanishes for £—2,, =; we know by 
Rolle’s theorem that А”„(6) will vanish at (n—1) points in (s, B); 
FU (t) will vanish at n—2 pointe in (a, b). Ultimately 











ea m9 -P 
но that 


a 





| 
k 
| 
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Thus the divided difference of (n —1)th order is expressed in 
terms of the (n— 1)thr derivative of f(z) at some point 7 in the 
interval (а, 5). Hence 





e акаа) "Чё, 
where £ is some point in (а, b). Therefore we bave 


[7] Вет) =). ( 72.) LOO, 


where £ is some point in (a, 6). 


This enables us to find an upper limit of the error committed 
in omitting the remainder term. 
Interpolation Formula—Forward 
3 jory-Newton formula is а speci] саве of 
Newton's formula 15.0(1) when the arguments are equidistant. 
Let 





s, = set(r—o, T=], 2,3. 





We have by §5.2(2) eta.) = уура" Mn), 


and by $5.7(2) ежа. m.) re. 


"Hence substituting these values in Newton's Interpolation 
Formula $5,6(1), we have 


ay fte) 114) + -ajaa + Ета) дада)... 





pOr а. le 28239 араада, | 


£ lies in the interval bounded by t 


least of =, 
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1t can be put in factorial notation thus: 


{= а) 
21 


(2) 1а) =1(а) + (2a) f(a) + A?f(a) +... 


2 Ania) + Ee 





or in difference notation this : 


= 4) (2 дэ, 
(8) f)» f(a) + (e—a) 2/9). p (55). ATA) ү 






afa) p 9 кың, 
It can be easily seen that the differences ERES lie оп a 
line sloping downwards from f(a). Bee table of $1.2. 
When w—>0, we obtain Taylor's theorem 


Hea) = Па) + (= ао) + 2 ry. 





(z—a)* 
GDI 
where £ lies in the interval (a, z). 
Putting а + ты in place of z, the formula (1) can be put in 
the simplified form 


(4) flo-+ zu) = Го) + zAf(a) + Vaya) +... 





"ое + TO gy, 





= 1). ењ +2 90-1779) 4 2671) (0)... (6 +1) a cn 
ررر ر‎ EME Bene apy, 
In particular, 
(8) Хань) = а) tnaro) + D A*f(a) e... ma *"tf(o) A") 
If f(z) be а polynomial of degree n—1. and в=0 and w=1, 
have ? 


He) =f0) sajo) Das) +... 267 D Cc. 
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5.9. Newton's Interpolation Formula—Backward Differences. 
Let us construct a table of differences: 





9-3 8%) 
Af(a —З=) 
a-w  f(a—24) Ajla 85) 
a Af(a—24) A*f(a Bo) 
а f(a—2) A? f(a — 26) 
Af(a—w) 
а fe 
Let z,-a—(r-l)e r=1, 2, 8,..., 
then by §5.2 faye) = Af(a— u), 


Жгут). = дда -2), 


and generally /(шугә...=,.у)= TA 'fía ro). 
"Thus we bave 


йа) Иа) а) + @—-a)df(a—w) + Ege tw) A*[(a— 20)... 


„(E-a +7 84) ааа 1а) 





+ (29) (=-0+ а= (=-а+ a) reg, 


where £lies in the interval bounded by the greatest and the 
least of =, a, a— (п —1)ь. 
From (1) we easily get 


(a) Қа + х=) = f(a)  zAf(c 


+1) (+2)... +2) д-лек 
Ета) 


ЕЕ en „ушу, 





мдааа) + - 


bonas = эү UNE Ps am Ре a 
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and 
e Гао) = (в) - zAf (a) += 






(ape D ma s-if(ac sla) 


риа) eng, 


‘The differences employed here lie on a line sloping upwards 
from f(a). The corresponding interpolation polynomial is obtained 
by omitting the remainder term. 


5.10. Lagrangoe's Interpolation Formula without Differences. 
Let f(z) be a function of the real variable = defined in the 


interval (a, b) to which belongs the set of arguments ү, gr... жь, 
and possessing derivatives. From 15,8(2) we have 


so ی‎ 
A ULCERS MUTE TELET) 


fixe) 


(e: aj. —). (ra E | 








‘Hence 


mE 
Q J-a nt) 






(ЕЕ >). 
(razza 723). 


"EST ER YER NE 
MOSCA ee уа а 








+ (2-21) (2-23). 
This is Lagrange's formula with the remainder term. 


Putting ez) = (2—21) (2-ға)... (2—23), 
wo have (1) in the form 


(z— 2 Merz Tg.. 








© 
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where Rule) (беу ояк) 

=e) ДФ), by 50.72) 
where £ is some point in (а, b). Thus 


(2) fe X Le) £9 1. (e. 
re 2-2, J) n! 
If f(z) be an interpolation polynomial of degree n—1, ite 


divided difference of nth order is zero, i.e., 














Hez...) =O 
and tho formula becomes 
® m= $ e 


-1:z—2z,' eG 


8:11. Hermite's Expression of Divided Differences. 
Hermite has given an expression of Divided Differences by 
definite integrals in the form 


a) f f asas tr Mundt, 
ө Јав e 


whore wu, = (1—)= + (f — fg)tg + ...+ (lana жа ua 
Evidently the independent varial tis fas бу... will all 
disappear when the definite integral is evaluated. 
Suppose formula (1) is true for the т arguments z}, za 
and let us proceed to the case of n+l arguments with z,, 
and a new parameter ta. Now 


ЕЕ и {а-о БРЕ 
° o 


ata + trans} ate 





O) færi 




















+ (tent 


HP (= ty dy +... + (fag — bea ny бна) 








Fa ant 
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4 | 
Hence f 4 чь f "Yas а 
$2597 


lega n) бева e ан) 
er 





which can be putas "f Esz1sa.-- o3) fir mg Fone 
Za Zarı 
бутрал) 
[7] ить. tutu c1). 
‘Thus if (1) holds for n, it will hold for т+1. Potting n=2 
in (1) we have 


1 
f ао etna] a, o nf opu, 
" "n 
so that the formula holds for n=2. Hence by induéfion it 
holds for я = 3, 4,..., i.e., it holds generally. 


5.12. Expression of Divided Differences by Contour Inte- 
gralis. Take ап analytic function /(!) of the complex variable # 
which is holomorpbie in and along the contour C enclosing & 
simply connected region in which the points эу 8. are 
situated. ‘Then by Cauchy's residue theorem 


ч) 10. a = 
f 2-і ets 
Again the residue at (= 2, of the function 


Kt) 
Шы 


На.) 
аа A т; 





is 








[n 
Hence 
1 KO 
а) 2A жышт Я" 
е, 


mI PLU 








= 
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Thus tho divided difference of order n—1 of f(a) is expressed 
in tho form of а contour integral. 


5.13. Newton's Interpolation Formula in terms of Contour 










is 2g) 
(as) 








1 
+ ef 


qoe 








+ atn neigt 







where. 
5 ГАВЕА 


which is Newton's general formula. Here f(s) is not unrestrict. 
ed as in $5.6 but is subject to the restrictions given in $5.12. 
8.14. Confluent Divided Differences. In the definition of 
divided differences the arguments were all taken to be different 
from one another and it can be exsily seen that if two or more of 
them become equal, the values of the divided differences become 


В T» ЧЕТ: 
de абык du „АРЫ Im. ow м 
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indeterminate. A limit may however be found of fizy, ¥; +4) 
as «—>0, nnd if this limit be finite, it can bo taken a» the value 
of fiz,z,). Buch limiting values of divided differences with 
two or more equal arguments are called Confluent Divided 
Differenc. 


Wo bave from $5.7 
т) 
ftr Gopi" 
where у lies in the interval bounded by the greatest and least 
Of 24, Eanes tae 

If all these argumente become equal to x, we have the 

confluent divided difference 
РЯ 
a) fers n) m ҮР. 

‘Thus confluent divided differences oan only be formed if the 
function f(x) possesses the necessary derivatives. 

For constructing а table, arguments having the same value 
are to be placed in a single group and (1) is fo be put in place of 
the corresponding divided differences. ‘The method is exhibited 
thus: 





=.) = 











= firi) Pe) 
n fen Jj array) 
ыу РА ia Kiko. oe 
z 
3 b ела) S бшка 
=. few) Nereus) 
Wea) Hartas) 
LOI fie) fzgzors) 
Heats) Жугуу) 
EX Elo с Meer) 
fie а 


"и all the arguments are equal to тү, we easily got from 
Newton's formula Taylor's expansion 


fee) + (nm) ل‎ e+ 


+ E peti) нш), 
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where Ryle) =" гё, 


£ being some point lying in the interval (=, ху). 
If there are k groups of equal arguments, p arguments equal 





to z1, q arguments equal to xg, т arguments equal to £,,...-., Wo 
from $5.12) 
r | ко 
Ney y.-.2g%s...2925...) 3 ЕЕ LE 








S | 1 x 1 gr! а ә 
(=I G-D GH 24 ài exp 655 











Н mn 
c mt rg) (=)... 
prre 


El 1 1 
-Di Q1) (r1)! Гб Tar 








= fanny). 





Ezamplel. Find an expression for f(z) in terma of powers 
of (z—8) from the following data: 


= 6 18 27 м 43 
te) 48 1627 17335 35626 78687. 
Constructing the table of differences and extending it to 
include repeated arguments for z=3, we obtain : 





* к) 
5 18-7 
-2 
8 f3-7 us à 4108) =6 i 
з 8)=7 18, 3, 5)=8 
18, 5)=18 1 
5 48 JB, 5, 18)=18 
БТЕ 198 s 1 
18 1627 2 1 
a7 17885 D 
Ж 2018 bin 1 
м 95620 uu el 
48 78687 
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Here the third differences are constant, equal to 1. 
Now 18, 5, 18)= 42—24x1.—18 
18,5) =198-10х18=18 
19) = 48— 2x18—'7 
163, 3, 5) = 18-10x1 = 8 
163) = 18- 2х8 = 2 
A" = 86- 2х1 = 6. 
Hence substituting these values in the above table and 
applying Newton's formula we get" 
Ка) =7 + 006—3) + 6( — 8)  (2—8)*. 








ExaurLES 


1. Shew that the divided differences of fz) ure unaltered 
if the arguments be each increased by the same constant C. 
while the corresponding tabular values аге left unchanged. 
[C.U, M.A. & M.So., 1939] 
2. If the argumenta be each multiplied by the same constant 
©, while the tabular values remain unchanged, shew that the 
divided difference of the nth order is multiplied by C-*. 
B. It f(e)=2"*}, shew that 
[mo 
[C. U. M.A, ё M.8c., 1034] 
4. Find the value of {(19) from the following table: 
є п 17 21 2з a 
=) 14040 88526 104486 270540 0923526. 
B. Construct a table of divided differences from the following 
data: 











= о 5 в 2 25 
f -179 1 18 87 7221. 

‘Extend the table to include arguments == 5, repeated as 
many times as may be necessary, and thus find f(z) in powers 
Ae qae Да expression for f(z) from the following data - 

= о 1 з 5 e 9 
Мә -18 о о -u8 о 18104. 





P 


© 
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7. Find an approximate value of log, 94-01 from the follow- 
ing dota: 3 
= 40002 40104 40288 40204 
f(z) 0-6020817 “0081877 -0045824 6052404. 
8. Мп, arguments are equal to тү, ng arguments equal to 
^» т» arguments equal to zp, shew that 


وم 
soft Water,‏ „$ 


Fa 





fez 





where y= (1—,)=у + (tı — tata + ...+ о-ва я, 











9. Form в table of divided differences from the following 
values: 
= 4 5 T 10 nu 18 
do 48 100 294 900 1210 2008 
and extend it to include the values of f(z) for ж=8 and z=14, 
10. Express the confluent divided difference /(aaabbo) as 
the quotient of two six-row determinants. 
11. Prove that 


flere, ed fs бо Heat ot Fra la dig. idle, 





where the integration is extended to all positive values, including 
[Genocehi] 


oro, which satisfy ty fg +......+„=1. 








CHAPTER VI 


CENTRAL DIFFERENCES 


6.1. Central Differences. Differences have hitherto been 
used which lie on а line sloping upwards or downwards from 
f(a) in the difference table. Now we shall employ differences 
which lie on a horizontal line in the table. Instead of keeping 
these as they aro, we shall use the modified notation ss employed 
by W. F. Sheppard.* He introduces в symbol 8 which is 


equivalent to ДЕ. The operator 8 is known as the Central 
Difference Operator. 


Since è = ДЕЊЕ Е 


we have A = 883. 
Thus Aus iu, Aug Puy, Auz, mug, ote. 


With this modified notation, the difference table сап be 
put thus: 
a-8e u 
г 83 
а-ә ча 
bug 





aro u 
to 





+2» чу 
bug 
а+8ы чу 
Tt may be noted that the horizontal line through w, contains 
the differences of u, of even order, and the horizonta) line 
E between u, and wu, contains differences of u, ‚у of odd order. 


© Proc. London. Math. Soc. 31 (1699), 459 
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The srithmetic mean of two successive entries or differences in 
the same column is represented by the symbol и, and is defined as 


M I 


Thus pig Muy +u). нау о +u), 
pu, = Hupp ue gy MMe = Зы tunog). 
pbug Mn. + ду), ёзи: — (us +д5и). 


pu, = Du, +8м,.}), мди, p = Miu, дм). 


Also В=МЕ-Е-)- 424. (1+2), 


and кзз = gute Вердан erp (1+5 D 


These mean differences can be inserted to fill in the gape 
between the differences from which they are derived, ‘The 
table will then stand thus : 


а чу bu, Pus (bus Suo имо ug 


puy дар ptuj Buy رقم‎ ати, ыйы 
atuu, abu, Suy pdu, Stuy pdu, Su, 





notation of central differences. Here we are not restricted to 
lending differences only, We can use the central differences in а 
horizontal line or in zigzag lines. 

6.2. Newton-Gauss Formula of Interpolation, We have 
from Newton's formula $5.6 
O fhe) = + (zy ered + urul mara) + += 

+ (uz )(u— eg) ы m ra аа-а) + Ralu), 

where Rau) = (и — 2 1)(9 — 72) (u — e) fua ame) 
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Put P 


Let u=a + zv, then f(u) = f(a + zw) 
By $5.2 we have 





i 
ауаз) ауур Аат) бурят СЯ 


and {жүз 





1 
)" gy r2" а) 





быа) 23)... ааа) 








5) xg) (u —2,)(u— و‎ 
rer r2)... (n (ее 72)... ет 


“(Tr Janie, 
Hence (ume J(u 25). (4-2, fira ciao, cl) 
(1 дее унны. 
Similarly (u—z,)(w— ә)... (му, caf oria Ра, +) 
(FH) ate remm (arty tren: 
Substituting these values in (1) we have 


(9) fla + eu) =/(a)+2As(a)+ ie Data) 





EEE) /A*[(a— e) + ree e asa ga) 
»*. (e ayes ale — 2) дэка 9) 4... Ratu) ; 


* Theaymbol 12) stands for the binomi»! ceefücient 
ETE 
ri 
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ог, in central difference notation 


(8) И ET «(i ) suo +) Эзу 


(=), + (uy ы 


1 
«(ет Jarmo (ipee Role). 
This is Newton-Gauss Forward Formula, In ît the following 
zigzag scheme of differences have been employed: 


а мо. Buy tuo 
в SM 
ou Pu 5%, 
3 à 
ste uj 
И т=2т, we have by $5.7(8) 


Roa(u)= (ех Jerre, 
so that 4 > 
(4) fla+ zu) =u, =u + E (ОЛ) 


CX (ime (m. ауан), 





where £, u lie in the interval {а — (т — е, a+ me). 


"Hn-2m-1, we have 


=. ы „+. 
O ея +78, (фу) 





a(g jeee (IR enm, 


where £, u lie in the interval (a—mo, a + то). 
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Ezample. Given the values of log sin 6 for 8—6*, 


compute the value of log sin 9” 2%. * 
Forming the table of differences, we have. 
бе 901925 
6666 
7 908589 — 809 
5767 200 
8* 914356 — 000. -62 
5077 147 17 
0° 90-19483 — 518 —45 
4534 102 19 
10* 0-230807 441 -26 
4093 76 
11° 928000 —365 
9728 


12° 931788 
=9°, w= 1" and z=} -0:30607. 
Hence by Newton-Gavss formula 66.22) we got 





log sin 9°22 = 919493 = 919483 
+080607 х 4584 + 16625 
-O3104x  —548 + es 
—0°05280 x 102 - 54 
*002100« —-45 — 10 
40701022 x 2 








Л. log sin 9*22/—9:21152.8. 


Newton-Gauss Backward Formula. Changing the sign 
of w in the formule of 16.2, backward formule can be obtained 
where z is measured in a negative direction from f(a). These 
сай be directly obtained thus: ы 
Put # =a, za, =a ru, zg,.170-Te and usa—ze, Thep 


Marea tare) = pater 8776—79) = eg nte 
fessa. gar 


" =н" 
(+016977 








= рна" eorr 
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Also (u—z)(u— zg)... (uza) 
= (1° (7-1) (ж+т—)...ж...(к—т)ы%” 


= (P7377) еды. 
and (u—z;) (u-z)... илз) 
=(-1° (æ +r) (а). 


-— (Са) атча, 





Hence 
(u—2,) (u — 23)...(u — zg )f(zr zo...) 


= (a )arte- re (**27 nw. 
and баа) ыы ) 

=. (255 Аза тә) = -( ы узи ү 
Bubstituting these values we have 
4) faz) = (a) -Afiat (2)a*fa— 9) (* 5 Ja*fts- 20) 


>< *(* 1! )atta-22-(*2?)a*fa- so) n 








(кух, 








+ (Cx? ват) (251A tarea ў 
This is Newton-Gauss Backward Formula. 
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Ifn=2m, we have 
в.) (** 5! Jeter 
where ё, м lie in the interval (a — mw, a+ m —1e). 


It n=2m+1, вые — (т Jet iron, 


where €, ч lie in the interval (а — mo, a + mw). E 


Again putting u-a-ze with z,=0, zg-=o-re and 
Fareimócre we got 


Kor no) mus mus жабы y+ (ро (ем 


+ 1") 8o + .., + ex емо + (уина ...+ Ralu). 
Thus for n=2m 


TM жи Ne reta. 
© fases) mem удае ("мо 


AER үл 
(ЭТ yen 
while for n 22m +1 


&) fore чао (т ця (уем 





+ (251 ptm. 


It should be noted that Gauss’ forward formula §6.2(5) has 
the same remainder term os the backward formula (4) if the last 
difference in each be of the same even order, 


‘Here the following zigzag scheme bas been used : 
fae lv 
эй!» БЫЗ Bou 
ae б. 7 rhe M HN n 
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6.4. Steffensen's Interpolation Formula. We have from 
Newton-Gauss forward formula 56.2 (5) 


@ weno +S (23 ente S (m уа, 





+ (EM eram o, 


where £, м lie in the interval (a—mo, a + meo). 
Now . 


GMa ryt уе 
ECE ym (ym mmy 
-À (IS e mum 
ищи. 
Substituting these in (1) above we get Steffensen’s formule 


( т) + A (ry 





(2) w,-w— 





+ (RT oer mme. 


Tt employs the following scheme of differences : 
439. "wa 
ёи 5^ 89м... 
> 5 1 Lx 
ч Ed 

ds Apo T уз 

6.0. Newton-Stirling Formula. Rearranging the terms of 
Gauss formula §6.2(2) in the form 


fo+ zo) = (е) + ={Af(0) -3A*f(a — 9) 2 A*f(a— u) 
283 73. (азаа) фаба 724) + “р Data- Be) +. 





= amemus. Durraxaxoas из 


let us replace the differences of even order within the brackets 
by means of the following relations 

A%(a—w) =Af(e)—Afla—w), 

Atja – 25) = A*f(a—w) — A*f(a — 29), 


Деј (а — 89) = A^f(a 24) — A*f(a —8«), 






TIe) Atr- f(s — re). 
Then we have 


А+ аа) 


(0 Каа) = f(a) +20 atf) 





3 —12) Affa) + A*f(a — a 
sus ctt) Ata e АМА 22 ,: mar. 19) A «(s — 2u) 








„ 262—12)53—22) A*/(a—2«)-- A*f(a— 8o) 
а 2 


азо, Bu) +... - 


This formula was given by Newton* and afterwards studied 
by Btirling.t 


In terms of Central differences it is put thus: 


(2) «uo + зиёмо + ө + 1 چو‎ 





EUM 
IR ENS 


a(z? 





22) ешо +... 





51 


Тыв result is easily obtained by taking the Arithmetic mean 
of Gauss’ forward formula $6.2(5) and backwari formula $6.8(4). 


© Newton— Methodus Diflerentiolie (1711). 
4 Stirling — Methodes Diferenticlis (17807. 
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Thus 
f(o-- ex) — u, +4 E EI пы 


н (E) (E) jare 


"(a ) mor me. 


whore Ё, м lio in the interval.(a mu, а + mw), 





=+ у жат 
ек 2+1) = — 


ia(z — 1)... (x7) 





rel 


2۵ 1(2 عا _ 
!)1+ 2( 











ко) e (527) e See. 


„223089 — 12) ) —0®).. (т 





(ат 
Hence 
CT 1 (= 1959 92). (etr) gare 
о-ва) =, “r. 2171 PLU 
aza — 1#у® —22).. (28 -F T?) gar 
e +, "ani atras 





+ 269 19)52 29). тя) „нен 


+ 
The differences employed here lie on в horizontal line through 
мо, namely, 
в, Wo, мдш, Bugo., paulo, Stig, имо... 
Stirling's formula gives more accurate results when = Шев in 
the range —3 to +4, i.e., $ on each side of the line through ug- 
Ezample. Consider the example of p. 147. 
Here the differences employed lie on the horizontal line 
through 9°, which is аз follows :— 
9° 0190498, 4805.5, —549, 1245, —45, 18. 


F EREL C = 
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Hence by Newton-Stirling formula §6.5(2) we get 


log sin 9*22'— 9710433 = — 910498 
4-086607 x 4805.5 s 1702.0 
40700722 х — мз - 86.5 
—005289 х 1245 = 6.6 
— 000485 x - 45 + 0.2 
+0"01022 x 18 * оз 


5 = 931525 
2. log sin 9°22'=9°21152.3, 
which agrees with tbe value obtained in p. 147. 


6.6. Newton-Bessel Formula of Interpolation. By вое" 
formula (6.2 (4) wbich end with в difference of odd order, we 
have Е 
() fle+evjeug+ “= (en) "y «Stn ym 


veo zn 








-1 
«CU 





mE). 


‘Tho corresponding backward formula $ 6.8(8) with the initial 
argument а + is obtained by writing z—1 in place of =. Thus 





à) flarewmuy+ (T) 


(tip) Gr emre. 


Now taking the Arithmetic mean of (1) and (2), we get 


ЕЕГ 


+1 








6) a+ ru) = pa, +S 





E үл atre ун) везна, 


le, а+ ть). 





where £, u fie; in the interval (a 

This formule was first given by Newton and afterwards used 
‘by Bessel und is known as Newton-Besrel Formule. Here ж 
‘should lie in the range 0 to 1. 


© 


154 CALCULUS OF FINITE DIFFERENCES 


Again since 
py + (z=) ug чобы, 
the formula (8) сав be put as 





(4) s, = мо + тё. 


_=-$ (ШЗ узаты, 


a* Aversa à 


(eir уме, «(tm reme. 


А more symmetric and wonvenient form is obtained by 
Putting ay +}. Thus _ 


(5) м, рени + S Xs (9+3 1) узнела 


+x! (" 2 1) Э " ( у ru 1) )etr2me 
Writing the first few terms in full we have 
O wy snag tying + ME „зм +0070 азир 
-pyt рз 
ы ыен + uly? ® 0 su, „> 
‘The differences employed here are taken from the scheme 


а Uo 
Bay pru, ч, pite, 85% 





ae u 
Example, Take the example of p. 147. 
Here the following scheme of differences is employed: 


oe 919433 
4594, —400, 102, —85.5, 19. 
10* 
Hence by Newton-Bessel formula $6.0 (4) we have 
log sin 9°22 = 919438 = 
+0-86007 x + 
—0°11611х  — 492 + 
+000516 x + 
*002100x — 355 — 
—0:00087 х 19 - 








eias Алат d ИГ 
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6.7, Everett's Formula of Interpolation. In this formula 
only even differences on horizontal lines tbrough мо and u} are 
used. It is obtained by eliminating odd differences from Newton- 
Gauss formula, 

We have from 56.204) 


Mam uo + ади, + x (271) fru, + (Grete Jarra? 


Pol bares erem. 
Now 


(7 yere (ES "е 


NIE 
— +! 








=) fare mts + tnnt 1 





Ets Mn forenses 7) ruy 


(ал а ("Z0 perme: 


Hence 
A) wm +e 47S (Л JF (HT уне} 


ZH 





Tem 1 yere. 


This is Everett's Formula. 


7-1 1-280 
фәр кс» a )-(* 
and uo + xu чый а) + хм, 


putting 1— == 3, we bave an elegant form of Everett's formula * 


* Everett, Brit. Assoc. Rep- (1900), р. 648; J. 1. А. BS, p. 482 (1901). Bee 
also Laplace, Theorie anal. des Prob., p. 15. 
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[7] umru + v (E prm 


tre Stt yen 


(m ye. 


where € lies in the interval (а 9 — 1o, a + ть). 
A second formula of Everett employing odd differences which 
Це оп horizontal lines through uy and wg. is obtained thus. 
We have from Newton-Gauss backward formula §6,8(4) 
$ (я-а s+) sar 
wm SÍ drat) ва) “} 


+( ELM Jenong. 


Now 


CE) MA + (435) otro 
-( = zr yemas (#57 (rug — 8=) 
= (257 yen би 


= (ця) erty. 
Substituting this in the above formula we bave 


mer AEG) Can) 


a 
(£n) me. 
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Hence 
* c 
(8) му + X (rm Dye-- - 
=i (9r-1) gan . 
a SHG) ea 


С а Т] 


where p=x+} and q7]—. 

‘This ia Everett's second formule. Both formule, more 
particularly the first, are specially adapted for use in statistical 
work. 

Example, Interpolate for the age of 87 from the following 
data: Я ; 

Age 295,  80,. 96, 40, ' 45 50, , $5 
Entry 1481, 1002, 1001, 2351, 2865, 8550, 4471. 
Forming the difference table we have 

= age ч, ? 





-2 25 мз Ө $ 
201 
-1. 30 1002 вз 
200 23 
o B5 1961 v 10 
390 эз 4 
1 40 291 124 24 о 
514 т 4 
9 45 995" 171 18 
d 685 65 
в во бм 296 - 
921 2 
4 5 тї 


For age 87, we have 2="4 and y="6. 
The coeff. for the first line are '4,— 056, *01075, 
j „ for the second line are ‘064, “01165. 
The differences employed lie on tbe borizontal lines through 
` O ала 1. „Hence by £6.7(2) the required value for age 37 
| (0 UA харо — 050 x 124+ 01075 x 14 
é +6 x 1901—7064 x 91 "01165 x 10 
meis —940:4—06:04-- 15 11760 —6782 11 
(0 200455. 2 















© 
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6.8. The chief advantage of Everett's formula lies in the 
fact that, for а subdivision of intervals, the terms are such that 
they ocour both in the z-expansion, and in reverse order in the 
next у-ехрапвіоп, ‘So we need only calculate these once. If 
u, is expanded seperately for each value of z, the whole advan- 
tage of the formula is lost. We illustrate it by an example 
below, This formula is specially applicable when the values 
to be interpolated lie near the middle of the series of given values, 
Everett's formula will be-wholly unsuitable for interpolating 
values at the beginning or at the end of the given series and it 
will then be necessary to apply the general advancing difference 
formula of §4.4, which again is unsuitable for interpolating values 
near the middle unless в sufficient number of terms is given 
(вее Ex. p. 94). 3 

Ezample. Calculate for all ages from 35 to 45 from the data 
given in the example on p. 157. 

We have from Everett's formula 


мези + (*5! Jaen +( #5? ity +... 


times ($3 ome e (¥ 2 atu +... 
and for the second interval 


ма) aus +( Учи, + 


ваа) ва (722) = 


‘The second line in u,,, is the sagrieds the first line in us. 

‘The values of z are “2, “4, 6 and 8 and consequently the 
corresponding values of y are “8, `6, `4 and "2, ‘The coeficiente 
of tho terms in tho first line of v, are = 




















it 
jii i 
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The terms involving differences of мо, му and чу are put in 
the following tabular form : 


























| Brd | Sum of | Sum of 

uat term aterm | rs Int fine | 20d line Result 

2 | 302 | -20 | -o | savas] 

4 7844 -509 10 77941 

6 | 11766 | -582 | -11 | 117089 
۴ "8 | 15688 -497 08 | 1564-51 
86 | -2 | 4702 | -307 | -09 | 466-82 | 1564:51| 2080:83 
87 4 $4074 -0m 15 93361 | 117089 | 2104:50 
вв | ‘6 | 14106 | -793 | -10 | 140288 | 779-41 | 218294 
39 |в | 18808 | —5%95 | 11 | 187496 | 389-35 2204-81 
41 | cmo —5°41 11 507 64 | 1874-06 | 244260 
42 |4 |11460 | —9'57 | -19 | 1186-82) 1402-68 | 258045 
48 |-0 | 17190 |—104 | 21 | 1708-27 | 933-61 | 2041.88 
44 | 8 | 22020 | -a21 | 15 | 2283-04 | 460782 276026 








"The column for tbe “sum of 2nd line obtained by 
writing in reverse order the values in “eum of Ist line "* from 
the previous tabulation. 


6.9. So far we have deduced the central difference formule 
from Newton's general interpolation formula with Divided 
differences und other formule derived therefrom. We «ball 
now try to deduce the formule direct on the assumption that 
м. can be expressed in the form of a series ond shall seek to 
determine the coefficients. Thus to express w, in terme of 
ordinary central differences lying on the horizontal line through 
Mo, assume 

м. = Ломо + А pug + Aa uo + A suo + А Î to +... 














+ Аа, 1583 “uo + Ags?" ug +... 


‘Then 








by substituting the values of central differences (£6.1), 
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E с rote Se 


‘Multiplying both sides by Е” and equating the coefficient of 
AF" wo have 


эл», tAn, = тег, 


Tyi 
Again multiplying both sides of (1) by E'** and equating 
the coefficient of 4*7*1 we have 


VET ET PME [сс уун эши 
Чолоо by subtraction 


> = (к +(®+7—1)......(—т+1)@—т) 


rer 
19)(z* — 2?).. 








„zt 





zz? — 13) (ж®* — 22)... (2 FT) 
Wima As ی‎ er 


ES. 
ЕС wo get Newton-! et formula. 
быч: 45 че v 

= 


aiet 


в. Again to obtain Newton-Béssel formula 
t Елените 











cestaat pirrERExces 181 
Then 
M Bem Ag+ dot AAt Aa S лаа +... 


+ ФА» 





жа, eer +48277 





Multiplying by E*-! and equ 
we have 





g the coefficient of 237—1, 


BS I 
17-1! 





7+1) 





449-24 Аз, 
Again multiplying (1) by E^ and equating the coefficient of 
A?” we have 


Vg et ápr Ao, = +71 








Hence 


and 


Hence 





wo mp + (oiu + Dua, t ا‎ аз, 





а-я, (+ elel) gon, + 


6.11. To deduce Everett's formula assume 
thy = Аомо + Ag Pug + A4 uo +...... + Agrê uo + ees 


t Age BP s +...... 





+A + Азы, + Ави, +. 
Substituting the values of the 3's we have 





a, As 

E*= Ao +4) +4, ۵+ (А+ Аз) tA р + 
. 

+ وھ‎ EE + 
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Hence multiplying by E*7! and equating the coefficient of 
A377 we have 






(z—r3) | (27-1 


Leer Diet r= 
Aari- t 2r-1 


Gr 
Similarly 


(rx Der 





бапа) афт 
T$) 






Aarit Agr t Aram 


[7 














(rer 3-2)... (аат) 
whence As «Ez 
а (у +т)у+т 
TU ша ESTEE 
Thun 


uam uo + Ci ye ("ye 





+ ua (5 ima + (T2 tn 
Other formule are also derived in a similar manner. 


6.12. Relative Advantages of Different Formulae. Central 
difference formulae are generally more accurate than Newton's 
advaneing difference formula, The greater accuracy of central 
difference formulae results from placing = in the central interval 
of the series of given values. The central difference coefficients 
аге as а rule smaller and more rapidly convergent than those if 
the ndvancing difference formula, and by a suitable choice of 
origin, the arithmetical work may be reduced to а minimum. 
Further the coefficients of the principal central difference 
formulae are all tabulated and can be easily applied. 

‘Advancing difference formula is useful in interpolating 
near the beginning or the end of the given series of values, 
while central difference formulae are very sdvantugeous in 
interpolating in the central interval. The selection is n matter 
of necessity rather than of preference. Successive approxi- 
mation to the required result may be obtained by the inclusion 
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of more extensive data on both sides while the basic central 
difference formula remains unchanged, but in the case of 
advancing differences a fresh set of leading differences.and a fresh 
set of coefficients will be necessary. This is опе of the principal 
advantages of central difference formulae.* When additional 
data are taken on one side only, advancing difference formula 
is to be used. 

For the bisection of an interval Bessel's form is convenient 
since the alternate terms are zero." Everott’s formula is very 
useful for subdivision of intervals in the central position (see 
$0.8. For values of z in the neighbourhood of one half Beasel’s 
form gives more accurate results. When z is quite small 
Stirling's formula will be found more convenient. See also 
$6.16. 











6.13. The Lozenge Diagram, Many interpolation formulae 
can be easily obtained with the help of Fraser's hexagon 
diagram. 


since (PGE 
we have 
w DCE) GE 


It in evident that giving to r and n consecutive integral values 
n table can be obtained akin to the difference table obtained 
from 

(2) At, Деи, - A^ 


Putting the three elements of (2) in the usual tabular form 
and the elements of (1) in a reverse tabular form, and combining 





* For detailed investigation ioto the relative accoracy of central and 
advancing diferenm formulae, see W, Р. Sheppard—Proc. Lond. Math. Soe., 
4 (1900), pp. 90-41 and 10 1911), pp. 189-72; D. C. Fraser—Journ. Inst. Act, 
50 (1906), pp. 25-97 ; ©. J. Lidstone—Tyans Fac. Act., 9 (1998), pp. 246-57. 

+ D.C. Fraset— Journ. Inst. Act., 43 (1900), pp. 205 et aeq. 
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them together, a hexago 
elements of (1) and (22. 





obtained having at its vertices the 





هکار 

enc Зам. 
Gti 

This is known as Fraser's Hexagon Diagram. As it isin the 


form of a lozenge it is also called a Lozenge Diagram. 
From (1) and (2) we have 


(2) [2-727] (т) 80] oe 


@ or (ути. + jan 
= (ати. + (н) 


Therefore if we assume that 








Аи, — — 














(i) the oblique lines denote the product of the quantities they 
connect and horizonta! lines addition ; 

(li) n live taken in a right-hand direction gives the product 
a positive sign, and in the left-hand direction в negative sign ; 


it easily follows from (8) that in travelling from the left-band 
vortex to the right-hand vertex of the lozenge, the sum of the 
elements lying along the upper route is equal to the sum of the 
elements lying along the lower route. Hence the sum of the 
operations in travelling round any closed circuit is zero, or 
simply the circuit of а lozenge is zero. n 

614. Placing in contiguity consecutive lozenges correspond- 
“ing to consecutive integral values of r and m, we get the 
following lozenge diagram : 
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By convention we have (*5')m1 for r-0, £1, £2, 
Consider any group of lozenges, The circuit of each lozenge 
is zero ond therefore the sum of the cirouits of all the lozenges 
will be zero. ‘The operations represented by the internal lines 
will cancel out, since each i» taken once positively and once 
negatively. Thero will remain only the exteroal lines of the 
lozenges and thus the circuit of any group of lozenges ia aero. 

It follows therefore that in passing from any chowen coeffi- 
ciont to any other chosen difference, the result is identically 
the same whatever route be taken. ‘Thus the following ex- 
pressions give the same result: 

мо + (мо + G)A uo + (PAu 

uo + (Ay + (AMA + (/59)A9u-, 

uo + (PAu + (Au; + ("5)2 

uy + (Ань + (Aug + (A-1 

u, + (FTN Aig + (Atu + (5 APH 

uy + (Su, + CGN AIH + (DATA 

(au, + (Au + (APs‏ + ي 

6.18. Suppose the diagram extends to the column ol zero 
differences, and consider the entire group of lozenges. ТЬе 
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pairs leaving the two extreme lines, and those on the right-hand 
side are all zeroes. Since the circuit of tbe entire group is zero 
it follows that the upper route from left to right is equivalent 
to the lower route in the same sense, 

We can add to and take sway any number of lozenges from 
the diagram. Hence every route across the diagram, bowever 
twisted or circuitous, if it ends at zero difference, gives the 
same result. 

‘Take now the lowest route through u,. The coefficient of 
1 and following the route slong Au, ,, Aupo, Дм... 
it is found that the coefficient of each difference i» zero. Thus 
the entire route through u, gives м, os the result. Hence the 
sum of the elements along ony route proceeding from any entry 
u, whatever to the line of zero difference is equal to u,. 

Thus 

MeHg Au o (PA u + (Аи y ("Atut «.. 
which is Gauss’ formul: 

wu, ctus + Au, + ( 
and taking the mean of the above routes we get the Newton- 
Stirling formula, 
16. Ву а slight modification in the lozenge diagram, for- 
mulae involving meon differences con be directly obtained. 













Instead of taking the product of the clements connected by 
oblique lines, take the product of tbe element at a vertex and 
the mean of tbe elements at the end of the two oblique lines 
meeting at the vertex. Thus to obtain the Newton-Stirling 
formula which employs differences lying on the horizontal line 
through мо, the route begins with wo and coming horizontally, 
to (f). multiply it by the mean of Au-, and Auo, and proceeding 
horizontally to A*u_; multiply it by the mean of (^2!) and ($) 
and soon. Thus 


omi s t ttn d [es] ats 


жеу Att اھ‎ „ 
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or in central difference notation. 


uz mue + aug ¥ 4248? uo + аа? — I) nS +... 

Newton-Bessel formula is similarly obtained by following 
the horizontal route through Aug. 

‘With regard to these formulae it may be generally remarked : 

(i) Formulae which proceed to constant differences are 
exact, and are true for ali values of x whether integral or 
fractional. * 

(ii) Formulae which stop short of constent differences aro 
approximations. 

(iii) Approximate formulae which terminate with the same 
difference are identically equal Й 

(v) Approximate formulae which terminate with distinct 
differences of the same order are mot identical. The difference 
between them is expressed by the chain of lines necessary to 
complete the circuit. Thus tho difference of two routes ending 


with Atuy and A?u. , respectively, 
(Азы 4 — Ази) + ("NA uy — A?u. 9) (89 uo Ам) 
= (Aus + (3!) g + AM 
=a series of horizontal diagonal 














equal to 











Mr. Fraser bas shown tbat formulae in divided differences 
can similarly be obtained by use of the hexagon diagram. * 
6.17. Central Factorial Expression. Product of the form 
ate + {п —1)(ж 4n - 2)... (£7 4n +1) 


is called the Central Factorial end is represented by the symbol 
etl, For n=O, а!®) =1; and for n>0, z!" =0 when 2-0. 





It eim be put ss 
a) > ai zie dn - 179, 


* D.C. Frasee—Nexton end Interpelaticn, а Memorial vol 
by the Mathematical Association, 1927. 
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"The Central Factorials of even und odd order may be written 
O M esf (2.1) 4)... (22 (m-1)2). 

*3 = ax(e? — 3) (29 — 3)......{2?—3(2m—1)?}, 


which shows that z!?*! is an even function and si-i) an 
odd function of x, 











=" 


Also 10-1 = дез", 
| "Taking the central difference of (1) we have 
al = (2$ (7 wy. (e-7 


| Чень) 
2 ие 





=(=+ 











| = (^ - 1) ae 
| +4 mA, _ 
і Similarly bat ии, 
(8) Sra ат), 





Ezi = в. 


648. Expansion of a Polynomial in Central Faotoriale. 
Let p(z) be а polynomial of degree т, То express it in Contral 
| Factorials assume 

plz) mas +a + 032121 +... 40,217) +. 
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Hence putting «=0, we have 
#0)=a and 4°9(0)=a,.r1. 


D 0) =9(0)+ S 73 iO) + 2 itg)... 





e etd. 
ET arg) +. 








87900), 


In particular, putting (ж) = =”, we get 
S zy. 
XU 





(3) =. 


"The values of 5—0" are known as Central Differences of Nothing, 
"They are obtained thus : 


rele) (rà BF) un 


= aeq) (+ 3-1) + ec nnw (i -3) 


0 r=) VEG) 
*car(-2)- 





k “зу =з 
=3-0(7)G 

| ‘Taking together the first term and the last term, the second 

and the gat but one term, etc., we find 

|. 8"0".20, when m+n is odd, ET 









moras’ (л) |". when m+n is even. 


o 
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6.19. Expression of a Central Factorial аз a Polynomial. 
‘The central factorial z!*) can similarly be expressed in powers 
of = by assuming ^ 


ate Хаз" 
„то 


and then differentiating it with respect tox, Thus we get 


w 





Now from §6.17(2) it is seen that the expansion of 212%) 
contains only even powers of = snd the expansion of z!2"*1! only 
cdd powers, Therefore D'0i*)—0, if n +7 № odd, Thus 


$ ротна 
tan) З gar DON 
+ «Из. ЖОГ. 


8 perge 
tamer) S arent 
e EU AN Gre 


"The coefficients in the expansion of x'*! have alternating signs. 
D2*0t9^! and D27*1012^*11 have both the signs (—1)""". 


‘The coefficients can be calculated thus: 








Binoe zima) = (x2 — т) fami, 
wo have Бу Leibnita’s theorem 
Dru t99-9) = (2? — из) D321 + (A7) D3 7 e P291 D(a — m?) 
 ()D2 7735127) D3(22 — m?) 


(29 — таура + 4те рэ 771g [901 + (2r)? 097-219". 





Putting 2—0, we have the recurrence formula. 
(8) — Darot$e-2) = 9r(2r —1) 09720021 — таратса, 
E 


with the initial values 





ay) Эт узабыз, Pump 
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То а similar manner from zí?"*!) we ha 
formula. 5 
(5) D9**19t29* 0 (2, + 1)2rD2 7-191217 
— 49m —1)* D? 012-1, 





the recurrence 


with the initial values 


(6) ротаи (= 1)" p= 





з peine 
} (2m-1jt P 
Putting z—1 in (2), we get 

$ Drom _ 





m В =0, т>1 
g Dre o-1 gei 8 425 а 
S dre (717 gg. (040984) (0m T 





6.20. Derivatives from Newton-Stirling Formula. Arranging 
the terms of Newton-Stirling formula §6,5(3) in ‘ascending 
powers of z, we bave 


(D fla + au) eu, = мо + диди — ибо + dani uo — y] opi o +...) 


+ tuo — ditus + Hou родино + 





> 
+ jê "uo 7 pê ug + понд uo —...) 





+ чө - 18u + goo...) 
+ овна dabas +...) 


таво Bt, + РЕ; 
+. 
Again by Taylor's theorem we have 


@ Газы) = Ка) + zwf'(a) + эзе” (в) + 
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Hence comparing the coefficients of different powers of 4 
їп (1) and (2) ge get. 
@) «оао Bid uo + dosi? uo — s road? +... 
фе) = Stu uo + duo — soi uo +... 
ноч Бача + rond Mo... 





wtf" (a) = 8*u, — [ug + so "uo —. 
Wbp (a) = Bug — hab? ч, + 
" 








"The above set gives the values of derivatives of f(a). To 

| obtain the values of /'{а+ жы), formula (1) bas to be differen- 

_ tinted r times. 

^ 6.21. Derivatives from Newton-Bessel Formula. Arranging 

| Newton-Bewe! Formula 66.005) in ascending powers of з, we 
havo 

a) weg Spay бы — Pa, + 04083, — festa + 121 


тар dita, + лоли — oo) 








Suy Ly t xls a m 
+ u EDI à ^ 
+ (uy AR +) 
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These determine the derivatives at points midway between 
the tabular values of the function. 


Derivatives of (a + (x + 4)w} are obtained by differentiating (1). 
Derivatives obtained from Newton-Bessel formula §6,6(4) will be 





© afta) miu] Buy + + онбеш байби. 
sfr) = puy — Lu, — dupla, + ق‎ + 


fra) =u — itu — 





уни) = piu, — 499} 


0а) miu — . 





Example 1. Values of {(6)=cos 8 are given for 0== 20°, 30°, 
40°, 50°, 60° and 70°. Find the value of f(6) for 6-45". 


Forming the difference table, we ha’ 
20° 298969 





30° 86008 02683 
"09999 зот 
40> "T6004 — 2826 os 
— 12895 372 -2 
Boe “64270 = 1954 es 
— 14279 р 
" 750000 = 
2 — 15798 


V 70” 84202 


Here „=10°, take а=40°. Ав tbe interval i» halved, we 
take the derivative from Newton-Besse! formula. By the firat 
formulg of §6.21(2) we get 


роту = = 12825-0 х :00372— „зу, x ‘00002 
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N w=10°= 7 = 174598. 
ow =10°= 5 = 174 


Ту Ё(45°)= — 77071 correct to four places of decimals, 
Example 2. Given the values of f(n) for =0, 2, 4, 6, 8, 10. 





Find the value of f*n) for n=4°8, 
‘The table i» formed below : 
о "0427876 
* 0268490 . 

2 "6091800 —8152 
255278 =312 

4 "6046584 — 8464 12 
246814 — 800 

6 "7198898 —8764 10 
238050 —200 

8 “7431448 я — 9054 
228996 

10 "1660444 


Here both Newton-Besse| and Newton Stirling formulae are 
applicable, and both give the same result. 

Неге ь=2 and z=.4, Take а=4. Differentiating Newton- 
Stirling formula §6.20(1) twice we bave’ 








арав) =u — Bug + apo + нчө 





Б Differentiating Newton-Bessel formula twice we have 
etf pite, + (а) + (2—1 — و(‎ 
= — 70008614 + "0000080 — "0000002236. 
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6.22. Numerical Integration by Stirling's Formula. We have 
by $6.5(9), 


ams + apto + u, + ЖЕ) „узш; +t 





+ 2681 — 





Integrating this, we get 


aq Judean, + PL + itus + iS 2 z yu, 


i (2* азу "i ("ба 5 
+a s) as (s catene 


1 (57.254. Aoi, 
+ a (7 he da? itus + eC, 


whore C is the constant of integration. Taking (1) between the 
limits z= —ğ and 2=4 we have 


fen * ditte - 





In the same way 











[T 
E 275 d — 37 gen, + 8%, 
$ WEST 5760 сото Жа 
Adding these together 
еде -jo iiu + agin Bom = 
CHINE Phu, = و5 - وة‎ 





88g + Bou, + Pug +... + Bu, BI, оф 829—1 р. 
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Hence substituting these values, we get 


e нон DET 


жа (Pea a) gr (Pa n4) 





* orosa( ПКЕ a)-- 





f (a) + f(a + u) +... e f(a + Tu) 
Ааа + re ge) — аба – 3)) 
= shkollat rwt Jo) - 8? /(a — 49) 
+ оно (8? Ка + т» + fw) —8°/(a— o) 





Example. Find the value of 


1 eot өдө. 
3 


from the following values of cot 6 for 6= 22°30, 27°30... 
Forming the difference table we ba: 





в cot 6 

22°90" 24140 

— 74030 
27°30 19210 417 

—'8518 —`0569 
82°80 16007 "0848 

— ‘2665 — ‘0302 
97:90 18082 “0546 

- 19 „ 7/0177 
42°30" 10013 70869 

— 1750 . — 0100 
47°80 09168 70200 

= 71490 











CENTRAL DIFFERENCES 177. 


Here LL Take 4—52*80 and r=2 Then by (3) 
we have 2 
1 3o. Cot &d6 — 1:569 + 13080 10018 + ,(— 1750-8518) 
— 50100-70100 + 0569) 
=3°9642 + у x 1703— sig 0400 
89714. e 


As tbe coefficients converge rapidly, fifth differences have 
been neglected here. 


x e cot 646 = 39714. = 9405, correct to four places of 
decimals. 


6.23. Numerical Integration by Bessel's Formula. Inte- 
grating Bessel's formula $6.6(8) 


4 
=De- 35 " 


+++ = 1)6-2)0—3) pita s 
between the limits О and 1, we has 


fae ty наем + lars — ид bonding + 
Та the same way 











feta pig fata Montag contu, +... 
_ вой generally o 
анну dant s eds y aont, y 
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Again since 
pu + dug + 





а 
we have by addition ` 
филе uo tus tut. +My du, 


биби, — иде) + pYo(u82u, — иба) 
— abl bol nb a, — мё°ш) +... 





fiir fue + mjd 





= (а) + fla + w) +... + f(a + r7 19) + fla + rw) 
— ly {pfe + гы) — nbf(a)) + Pts {ns (a + ru) — из f(a) } 
онёо + ro) — За} +... 


"his formule wil be very useful for the calculation of 
perturbations. 


Example. Calculate the value of = сов 646 from the duta 
given on p. 105. v 
Here w=4° = 177000818. Take а=28* and т=2. 
From the table on p. 113, we calculate 
ni сов 36° = 041002, ид сов 28° = 10327485 
ph" cos 86° = 000200, oo» 28° = 000159 
ы A. (icon ож} cos 28° + oos 82° + cos 90° = 10040005 
— yl сов 36° — и8 сов 28°) — -0000878 
+ sMs(ub?cos 86° — иб?сов 28°) + 70000000 
+ тоозззв 
^ ffs. cos өдө = 176030483 x 080818 = 1182173. 
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Examples 
l. Prove that 


O Eh = (+ ie 
6) Blu. { sitne 


Gn) зи = (1 + Jue: 


(iv) Du eig [8+ (1+ Fa) ] ue 
2. Prove that 








o Scn»-(-))* =0. mL 


(ant 
ger 102911 


à Зи} I. 
® xn 297 берт 7% "20. 


a. Given ug 22854, u, =8102, м, 35144 and шуу 30909, find 
the value из by Bessel's forinula. 
4. Given  logD10224018017 log 340225814789 
log 320=2°5051500 log 350 = 2:5440650 
log2330:2:5186189 ^ log 300225508026 
Find the value of log 3375 by Newton-Bessel formula. 
5. Use Newton-Stirling formula to find usg from the data 


изо = 147035 155712794 
ug 19674 go = 12:089 
39719267 45711800 


Explain bow you would alter your method if you are required 
to find all values from ugo to tgs- 
6, Establish the formula 
uy = Met 4 Du а+ C 0 р 


iet Duy + Met HE us 
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т. Shew that 


+ E-ra] ر دوب‎ аа, 2022—1024), 
dt Beye Fate THE Dae „= Dato De. 


MET Ea p+ м ее Ре 





ME*—E-*) =zuê + rd ezet-net-) us os 


З(Е* - Е-*) zi seat ane SP DOP = Ds 3. 


B. Prove that 





1 в? per 
uyang gt ga 7 gigs Yt 


9. Sbew that 
چا د رر ر‎ anu oA RPO, 


10. Find the value of file) and f(x) for 2=116, having 
= 105, 110 115 120 125 
101 1758 3080 2454 2888 3871 
11. Show that the lozenge diagram method really derives 


all the interpolation formulae by repeated summation by parts, 
i.e., by the use of the formula. 


маде, 2(и,е,) едм, 


12. Prove that. 
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18. Prove that 470" is the coefficient of Lz* in the 








expansion ot. (4-74). - 
14. Calculate /(75) and /*(67) from the following date: 
во 409434 
45 380006 65 417499 
50 3791202 70 424860 
56 4700733 ^5 481749 
16. The equation of a curve being у=? (Fein 2, find the 


area included between the curve, the axis of x and the ordinates 


16. 


хен and rede. 


Evaluate Фоа, зан 


by numerical integration, 


17. Caloulate ("Кода from the following deta: 


E 
1 
2 
a 


fx) = 1) 
22761 4 147286 
1720480 5 18870 
1537052. 6 — YVT245 


E 
7 
в 


fte) 
1790053 
211045 
237089 


z 
18. Find by Bessel's formula the value of (3 6 sin 62е. 


19. Find the value of /'(z) and /"(x) from the data given 
in Example 2, p. 174. 


variable; the required length 4 





= G уз—«#соа#Е dE 


where ¢ is the eccentricity. 
21. Find by mechanical quadrature the length of the portion 
of n cardioid r=1+ сов 9 comprised between the initial line and 
‘a line through the pole at right angles to tbe initial line. 





CHAPTER VII 


RECIPROCAL DIFFERENCES AND OTHER 
METHODS OF INTERPOLATION 


7.1. Reciprocal Differences. Suppose tbe values of f(z) are 
given for the arguments zo, у, =,..., Za, no two of these being 
equal. We define Reciprocal Difference* of the first order of f(x) 
for the arguments zo, z,, by 








me EL. 
a) ЕВ 
Evidently (1) is thé reciprocal of the divided difference 
(вожу). But the reciprocal difference of the second. order for 


three arguments is not obtained by the repetition of the operation 
denoted by р. It is defined by 








РЕР 
[7 Рагоза) = ER] tle). 

‘The operator p does not obey the index law or distributive 
law, Similarly reciprocal difference of the third order for four 


argumenta is defined by 
re eL Pa не 
азиз) Рула), 


and generally a reciprocal difference of the nth order for n +1 
arguments is defined by 





Te 
Reciprocal differences, like divided differences, аге rymmetrio 
functions of the arguments. А 


* T. N. Thiele—Diflérences récipregues (Sterber. Akad. Kopenhagen), 
ICCB, pp- 1823:1; Feterpolationsrechnung, Leipzig, 1900, pp. 199148. 
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A table for reciprocal differences can be made in the following 
manner 


o Ито) > 
prora) 
z fle) Работа) 
работу) 
за Ја) рабунак) paro угала) 
pizza) = Par rara) 
x les) patrata) 
pisa) 
=. e0 


То obtain the value of a particular entry, way pg(F as4), 
take the lines sloping upwards and downwards through it. Take 
the difference x x, of the extreme arguments тү end #4, 
divide it by the difference of the two adjoining entries in the 
previous column, namely pa(r\7%5)—Pe(tatsz4), and add to 
the result the entry which forms a parallelogram with these 
three entries, namely р(х). Thus we get 
Ды 
) - рабата) 





+ pizgzg). 





59 





Erample, Form the reciprocal differences of /(2) = ®. 


a P Pa Ps 
о о 
1 
1 -2 
à ° 
2 <u 
1 о 
-28 
B K$ E P 
4 16 -17 
2 à 
.5 25 


Prom the above table we see that reciprocal differences of 
the third order of z? are zero. 
үл 
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7.2. Thiele's Interpolation Formula. Putting x for =7 in 
$7.1(1) we have x 


~ ры) = Dc 








or = fey) + 2221 
feme) TÊY. 


Similarly from §7.1(2), 7.143),... we have 





(ке) = p(x yxy) + — — 13 
Pede lei) + tame) For)" 





рая) PR rara) + — ч 
teeta Pals eat a) + СЕТ СПА) 


.—— و 
Pa Urz aar у) —ре(жүхөгу)‏ 





Pa(zzyzars) e )وم‎ rar ул) + 


"робити rans) palarat) 





rat yt ym.) = paz armas) 





and vo on. ï 
Hence f(x) is expressed in the form of a continued fraction 


= 
а) fe) =e) * e se 
Pal 





z-z, E 

iD anc = 

Рет) "primal ЕЕ ОЕ 
ase.) 


а = T z-z, 
иаи) pain Mr omm ET PE mm] 


‘Here we have taken n interpolation points. И the numerator 
‘of a constituent partial fraction vanish, the continued “fraction 
stops there. Putting ===, у... x. in succedsion, we get 
from (1), 








Каз). feras... Цен) 
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and, аз the last constituent partial fraction does not affect the 
result, the same can be omitted. Thus hy Thiele's Interpolation 
get а rational function agreeing in value with а 

given function at the given points. 
If f(z) be в rational function, (1) will terminate and the 
continued fraction will be identically equal to the given function. 





Example. Interpolate the value of /0275) from the table in 
P.188. Take z,—2, and by the application of §7.2(1) we get 





=4+2°25=6-25 
agreeing with the true value. 


7.3. Confluent Case : - Reciprocal Derivative. 


"When any two or more of the arguments are equal we get 
the confluent Reciprocal difference, the most interesting саве 
being when all the arguments are equal. This latter case in 
known as Reciprocal Derivative and is denoted by 


a) Tal (2)= 





m pu(zzz. ia). 
Putting n —1, we Бате 





o ЕРЛШ СЫ 


D Sa derivative of the first order is the reciprocal 
of the ordinary derivative. 





* т. н. Thiele—Interpolaticnsrechnung, Leipzig, 1909, p. 192. 
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Again ndm Sm „базы 









Temps (A 
=flz)+ lim [2 Ha) а =) 
Zn tn ese 23) рату) 
‚ =fe)+ dim [breui up 
pa (CTE Thur 


ЛЛУ {rer ] 


=fa)- mem, 
Similarly 


Talla) у э” [== сн ] 
та general 1 


al 
rete) lita | 


=a (a) 
PEE sa)” 





trennen] 


=. (xy — 2g) + (29 — 2; +...) 
Taal) sim m Tam 





жск 
a a ell ЭЕ o ERE Em 


z-s, = 
“шырай "ЛӘ 


@) Palle) ти) mr raf Ct). ^ 
which is the recurrence formula connecting successive қесіргосаї 


derivatives. 
Putting /(2)=2", we got 
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7.4. Thiele's Theorem. When all the arguments are equal, 
we get from Thicle’s interpolation formula an expression of the 
function as a continued fraction involving reciprocal derivatives 
at a point. 


"The kth divisor in Thiele’ 





interpolation form 











Раба... Za 61) = Рь- (2: a) 
and in the limit when ху. 25,..., 24.,—>=, we have 
limp, (хут. 1) 7 раа.) 





afis) — 





afr) = kr r, у](х), by $7.3(8). 
im $7.2(0) we have Thiele 





Hence putting k=1, 
theorem, 





M fle+h)=f@) Tx Je 





AT A 
ов 
rra TIT... 


the lant partial divisor being 


lim palz+h, ж, ar. e) — re aft). 
алта 52 


This continued fraction will terminate if the function is 
rational. Taking /(2) = 2, we have 





[E 





The continued fraction corresponding to f(x) can be put as 


@) Merge = 
mim. BOF Lr орача 
Example. Express e* as a continued fraction, 

We have 





Sad US easily proved by indogiion that 
rae m(-i)et © rase t e CH Dt (n De". 


ire. 
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Hence substituting these values in (2) we bave 





which can be put as 





et-l4 





i= 2+ 3- @+ б- 2+ 7- 
7.5. Aitken's method of Linear Interpolation by Iteration.* 


Let für; a, b, c, 4,...) denote the interpolation polynomial 
corresponding to the given arguments a, b, c, d,... and coinciding 
in value with f(z) at the point ‘Then by Newton's 
divided difference formula 











) = f(a) + (x — a)f(a, b). 
) = f(b) + (2 — b)fía, b), 
trical with respect to the arguments. 





Also 


since the function ів sym 
Eliminating f(a. b) we got 





т b Ка) а-а 
п $9, om 
) fiw; a,b) | dere 


‘This is a first difference approximation to f(x). It is o linear 
interpolation polynomial. 


Again f(x; a, b, e) f(a) (s —a)f(a, b) + (#—a)(x—b)f(a, b, c) 


+ (b-a). 








изв, b)+ (e) БД, b, о), 

and еза, b, defies a, 0)+(e—a)(e—cffa, b, с). 
Eliminating (a, b, c) we get 5 

CS я A РЭО 


Қа; а, с) с-а 


* А. С. Kitken —Proc. Edin. Math. Soc. (2), 1089, р. 86. 





9 
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This is a second difference approximati 
proximation to f i 
quadratic interpolation polymonial. Я MER 
Similarly о 


f(z; a, b, o, d) f( ; а, b, c) + (z—a)(x — b)(x — е) а, b, с, d) 


=f(z ; a, b, d) + (z —a)(x —b)(z— d) f(a, 
к {= — d) f(a, b, e, d) 
ге, b.e) e— 
(3) fla: a, b, с, d) = Tan tae + (4-с). 
| fat. 4) dx 


‘The &bove method can be shewn in the following table : 








Argument Function (HD 77 D Parte 
a f(a) А a-z 
b 10) j bon 
е fe) fle; a,c) fie: a, b, c) e-z 
4 кч) „4) {= a,b, d) fle; a, b, o, d) а-а 





Each entry in (1) is obtained by cross-multiplication of the 
corresponding values of the, function with the respective parts 
and divided by the- difference of the arguments. Hence it is 
called a linear eross-mean. Each entry in (2) is obtained simi- 
larly by cross-multiplicatioa of respective entries in column 
(1) and parts, and division by the difference of arguments, and 
soon. Giving tox any definite value we get the approximate 
value of the interpolation polynomial which coincides with the 
values of /(z) at the given interpolation pointe. For tho exact 
value we want the remainder term which has been omitted. 
Thus it four arguments a, b, c, d are given, we have 


fis) fo ; a, b, с. +60) 008—0) PA, 








where £'lies in the smallest interval containing a, b, c, d, £- 
Sometimes it may be necessary to replace the entries under 


Parts by numbers proportional to them. The process will 
terminate when further interpolations do not influence the result, 
Pam 
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Езатріе. Find the volue of log sin 8* from the following 
values of log sin 6. * 


Parts 
6° 901928 * -2 
T° 908589 915255 =1 
9° 919433 913596 9144255 +1 
10° 928067 012945 — 914485 — 914360 +2 


11° 928000 9123778  9:145357 9149704 914357 +8. 


The entries in the columns are calculated thus: 





901928 - 2 901928 —2 
| + 1=9°16255, + 82913506, 
908589 —1 С 910488 +1 
and во on. 
‘The successive approximations to log sin 8° are 
9715955, 9144255, 914866 and — 914857. 


‘Thus log sin 8° = 914857, Й 
the truo value to five places of decimals being 09-14360. 

7.6. Aitken's Quadratic Process. Suppose there is an even 
number of symmetrically placed argumenta 





k, 
=k; 
then 
я) де; ө, em EAI eNA 


js expression remains unaltered when —o із substituted 
for а and hence itis an even polynomial in а. We denote it 
by fiz; а). Thus 


(2) =; а) = 








[^ 
Ж-а) z+e 


Putting a—z, we have f(z ; =) =/(=). 


|+ = 
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Now we can proceed with 63, c*, d?,... asin the process of 
linear interpolation of §7.5. Thus - 
Næ; a2) аа 





(8) Ше; аз, b3) = 





| f; b?) 2-22 
[е a2, ъз) bout 
Ге; оз, езу ad и 





‘Thus we get the following table: * 

a?) 

b?) f(x ; дз, b?) 

c?) fiw; а?, с?) fix; a*, b9, сз) 

We; d?) f(x ; а?, 4%) f(z; a”, з, d?) fie Ъз, сз, 4%) 12—21 
If there are 2n arguments, the remainder term will be 











Ge? — 3)... (x? a. 


This process improves the convergence to a marked degree. 
When the arguments afe equidistant, these will be sym- 
metrical about в middle argpment or central origin, m say, and 
hence subtracting it from each urgument and dividing by the 
interval e, wo get the prepared arguments in the form +‘, 





15, 4 25,... and z is replaced by 





Aitken's methods possess particular advantages. We can 
stop at any point and thus any required degree of approximation 
can be obtained. When the intervals are unequal, the linear 
process bns considerable practical advantage. Рог equal interval 
the quadratic to be preferred. Both these methods 
are very useful for inverse interpolation. 

Example. From the following values of f(x), find the value 
ot f(61). 








58 "26011 
55 782251 
57 "0671 
59 749608 
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‘Transferring the origin to 50, we have a3, b=5, c=7, 4=® 
and =1, x 


‘The first column_is formed by {7.6(2). Thus 





9101 -2 
fr; a?) = | + б=-23119 
"6185 aj 
32251 —4 | 
fæ; b?) = + 10-7249358, and зо оп. 
‘1393 в 


‘The second column is formed by 17.0(8). Thus 





"29119 в 
«= 52) = + 105 7222106, and so on. 
240358 34 
‘Thus we get the table 
"20119 8 
"240858 “222106 a 
27898 291082 "22258 48 
"вазала ^^^ “220032 ™ > - -22200 "22255 80 


‘Thus f(z; 03, b?, e?, d?)—-22255, agreeing with the true 
vulue correct to five places of decim 





ЛЛ. Neville's method of Iteration. A slightly different 
muthod has been given by Prof. Neville, He brings the ‘Parta 
in the second column of the table. The advantage of this 
method is that repeated arguments ore considered in this scheme 
and hence the derivatives of the function whose values are 
known аге also taken into account, 
` The table for unequal arguments is represented thus: 
Argument Part 


a «-a fla) 





Қа; ab) . 
i ; a,b, 
b z-b fb fas ho) fai e, e 
е z-e fle) arena) fle; b, с, d) Е 
=-4 №) os E 
тї bi 
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To find the value of say entry in the table, we take the 
lines sloping upwards and downwards from that entry and 
terminated at the parte. Thus 











в) Me; a, b, с) |: $ 1 o leno, 
—е jiz; b,c) 
z-a f 

a Ne; a, b, o, d) 
#-4 fle; 

and so on. 


For the confluent case, let a=b, Then by $5.14(1) 
we have Ка, а) = f'(a), 
and we write 
[7] f(z a, а) = Ка) + (= —а)/ (а) = f(e ; а). 
Similarly if three arguments are equal to а, we write 
(4) гра, а, a) m f(x ; a®) + (x — a)? f(a, a) 
= Ке; а?) + Mo — arn) m Це: аз). 
"The table is represented thus : 
аа-а f(a) 








ferat „ 
E ie; 0 
mti hol? аз) fle; 0%, b) 
aso Қа 5 
i (е: a2, b9) 
bana soy"? seit 





a, b*, 0)‏ $ و 

bz-b fib) 

ET 
"The entries are calculated by the relations (1), (2), (8) and (4). 
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Example, Find tbe value of /(°84) from the values given. 


"80 04 «оваз 
-01553 
81 оз * -00035 -01613 
01598 "01618. 
B5 o1 01278 "01018 ers 
01023 
"6 —0 “02318 


‘The entries are calculated thos: 


| 04 сооззз 
| 5-701 = 01558, and ко ов. 
Î ов -00638 

The value of f(54)9'01615, which agrees with the true 
value correct to five places of decimals. 


7.8. Inverse Interpolation. The problem of inverse intere 

polation consiste in finding, from a given table of values, the 
value of the argument corresponding to a given value of the . 
function. In direct’ interpolation, however, the value of the 
function i» required for some other values of the argument. 
‘Thus suppose the values of the function у f(x) for £), жу, z, 
эге уз, Уз, Ys... direct interpolation as we bave seen consiste 
in finding the value y, for zer, while inverse interpolation 
consinté in finding the value x. corresponding to у = yn. 


By Newton's formula we have 








fi) = НО + дуо) + E дадо) EEIE a50) +... 





and there is only one value of f(x) fora given value of #, But 
sonwidered a» an equation in z, it gives multiple values of x for 
each given value of f(z). Thus, if f(z) be of the mth degree, we _ 
get п values of x. Here we are faced with difficulty of the same 
nature ва in the case of direct interpolation ($4.1), In order 
to have n result consistent with the given Cable we asume 
‘that within the range of interpolation there should be one value 
‘of = corresponding to a given value of y. That i» y=/(z) 
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be а wtrietly monotone increasing or decreasing function over 
the range of given values, and tbat the required value of = must 
lie within the given range near the centre. Then the inverse 
function z—f-!(y) will be в single-valued Tunction and can be 
represented by polynomial 


79. Vor practical purposes we must bave some knowledge 
of the nature of the graph of y=f(x) and must form a rough 
estimate of the position of ж corresponding to the given value 
of y, and consider the values in the neigbbourbood of this 
entimate. 


Alternatively the columns of arguments and fonctions can be 
interebanged, and then the problem of inverse interpolation 
becomes that of direct interpolation for, unequal intervals, and 
hence the methode of divided differences can be employel to 
Obtain any intermediate value of = corresponding to в given 
value of y. 

When the given value of f(z) is zero, the problem of inverse 
interpolation reduces to that of numerical calculation of roota 
of tbe equation f(z) =0. 


7.10. Inverse Interpolation by Successive Approximation. 
Writing Gregory-Newton advancing difference formula 
the form 
I(x) — (0) = xAf(0) + ele — 1)A2f(0) + рае 1) (а AO) +... 








we bave 





= =: E 
№ aeaea o 





-De -2)A*/(0) + 





Since x is small, » first approximation, ац, i» obtained by 
neglecting, terms involving differences of second and higher 


order. Thus , 
=t) КО. 
E 


las. 
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For а second approximation, ay, substitute the value of a, 
in (1) and neglect differences beyond second order. Thus 


28 Ка) -1(0) ч 
= ARO) + (оу —1)АЗ](0) * 
Similarly а third approximation a is given by 
Tom ff) —— 
“АО Йа» 174970) + Ma Тау 
and so on. . 





DATO” 


Ezample. Given the present value of an annuity certain for 
я specified number of years at the following rates of interest 
Rate 2 r5 8 85 4 45 
Value 16:8514 15:5802 148775 142194 18-5008 18-0070 

Find approximately the rate of interest when the present 
value will be 14°60б1. 





—`ооз9 
0008 
—'0086 
70008. 
— ‘0088 
45 180079 
Неге the interval of differences w="5. Taking the origin at 
2, we have 


1400012162514— 722, + 0%0бе-—1) _ "009021 1e 
= 
0908s 1e-me-B). 


41 - 
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A first approximation is given by 





-боб1 — 16-: e 
“= 14:0061 710-3514 -2AM -2990.. 
For a second approximation, including tbe next term we bave 
aim —r7458 AM 
27 — 7022 + (22896 —1) 0505 — 28921. 
Similarly 
dam TUE s 

27 T6: 72-8921 —2)-0089 

2 = 28000... 


"The required result is 2-424 фл. 
‘Thus the successive approximation to the required result are 
3:1449. 519005. 5-1998. 
the true rate of interest being 8:2. 


7.11. Inverse Interpolation by Central Differences. It ів 
convenient to use a central difference formula as the basic 
equation. The coefficients will then be smaller and will converge 
more quickly, In most cases a first approximation will be 
more accurate than that obtained from advancing difference 
formul; 

If æ lies between —} and ], it is best to use Stirling's formula. 
and if z lies between } and f, Beasel’s formula will be helpful. 

We have from Stirling's formula 











v, uo m лиди аазы + Lele? Iuno + date? - 1) 


Dividing throughout by põug we get 





DUM 
Lr e 
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For а second approximation аз. substitute the value of a, 
in (1). Thus we get * 


d q-e _1 -1y 9e _ 
Бете — atat Dae” 





А third approximstion as is obtained by substituting this 
value of аз in (1) and so on. 
Again writing Stirling's formula in the form 


@ == ores Torres TET 








Dp uo + 59+. 
the successivo values can be obtained in the manner of $ 7.10, 





and so on. 





Bzample. Consider the example nt p, 196, 
Take the origin at 3. Tbe differenees to be employed lio on 


the line А 
B, 148775, —'6884, 0466, —00375, 0003 


Ву (1) wo bave 
34'0001—14:8775 ,1 20406 ,. 1 -00875 ys 
ee Ae EE o ر‎ 712 e 


Thos am 2M = 3042 





чата +) gasy MDT- PRIS LIM) 1) 
= "3942 + 70052 + "0008 
= 3007. 

‘Thus the successive 


ate . 
31071, 8719085, 

‘The first value here gives в closer approximation than the 

second value obtained in p. 197 and the second value than that of 


CUTS ~ oS aa 
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the third value. Herein lies the importance of central diffar- 
ence formulae, =. 


Again using the formula (2) we have -~ 
oy = 3042 


Thus the successive approximations are 


31971, 310979, 
the truo value being 372. 


7.12. Inverse Interpolation by Divided Differences. Inter- 
changing the columns of z and y, the problem becomes one of 
direct interpolation. The table is put this 





Argument  Funetion 

Lr = 

у i НЕТУ 

z, 
3 Nyala) we fisa 
LE fs Kyavave) 
“Kyaw 
n LN 
2#, 

whore Кутуз) 711723, eto. 


By Newton's general formula $5.6(1) we have 
Q) хех, + Yaa Va) + — iy —ualfüryoya) + + 
+ (y-y) va). (yn. lvy: 
HY iy — val... (y у) yy Yar Yu) 
For numerical calculation of roots of an equation we have 
(2) жж: ~ya) + Yafa). 
Confluent divided differences present a practical method of 


solving algebraic equations especially when the degree is low. 
Thus when e 








mus‏ ووکرو 


wo have by U) Казаа) in 
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so that 
e dem fir) + (0 — frin) + & un? Pet. 
ME 
+ EU а-а) + Ru (s). 


Putting # =0 in (3) we have 
ч) Jte) = f(0) + уко) + +..., 


the coefficients being divided differences with repeated argument 


zero, 


Example 1. Following values of y= f(z) are gi 
я -42, — 39, 18, 78 
у 50, м, 54, 55. 

Find the values of = when f(z) = 03, 57. 

Interchanging the columns of x and f(z), the difference table ix 





v = 
20 -42 

3 
51 -59 4 

19 1 
54 18 9 

55 
55 73 


Hence by (1) we bave 
ze —42+3(у— 50) + 4(y—50)(y—51) + (у — 5О)(у —51)(y — 54) 
Л For у=58, webave z= —42+3.3+4.3.2+3.2,-1=—15 

and for у=57, wehave == —4248.744.7.647.6.8—278. 





Ezample 2. Find the real root of the equation. 2? +2—8=0. 
Firstly, take four arguments equal to zero. The portion of 
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The third divided difference is constant, =1, end hence the 
last column is immediately filled up by 1. Other entries are 
filled up by the relation А 





Таа) огур) + (re — хо) (oF =. 
The table stands thus = 
о -8 
1 
o -8 o 
y a 
о -8 o 
1 1 
° -3 1 
2 1 
1 -1 22 
404 1 
12 — 7072 8741 
5/8561 1 
121 = ‘018439 863 
54287 1 
122 7035848 365 
54652 1 
1'22 *035848 306 
64652 1 
122 7085848. 90534112 
4 54411284 
152184112 — -0000025 





We find the entries corresponding to 1, 
give /(z) negative sign. But х= 122 gives f(x) a positive sign, 
wo the required root lies between 1°21 and 1°22. Repeating the 
argument 1°22 thrice and calculating the entries, we bi 

085848 + 6:4052(z — 1°22) + 3-бб(ж— 1/22)? + (= 1-22) =f(z)=0. 

Neglecting the third power of (z—1:22) and solving the 
quadratic, we bave 2172134112. 

Inserting this value in the table and calculating we get 


А 1172184112) = — 0000025, 


eo that the root is slightly less than the true value, which 
is 12134116. Thus the value obtained is slightly in error in 
the last digit, the error being "0000004. 
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7.18. Inverse Interpolation by Aitken's method of Iteration. 

Tbe methods of interpolation by iteration as explained in 
$7.5—7.7 are very useful for inverse interpolation As the 
Arguments need not be equidistant, we can interchange the 
columns of argument and fonction and thence the required 
result is obtained by the method of iteration. 

We give below the example worked out by Aitken. 





Exampl& Find the positive root of the equation 
ч 27 +2824 —480=0. 
It is easily seen that the root lies between 1'0 and 2. We 


have to find the value of £ correspondiog to y=0. Calculating 
the values of y for == 150, 1°91,... the table is formed thus 


у = 


— 257140201 190 
—14:0254167 191 23189586 


— 873074639 102 2952228 | 2882804 
+ 82499435 ros 2710920 87312 84138 
«900820890 — 1:94 2483078 91702 17 68. 


Since y=0, the column y reprerente the “parts” of the 
method. The entries in the columns are calculated thus: 
190 —257140901 | 
эзїнөзөв = + 110886004, ond so on. 
191 —145254107 
Thus we obtain z-1'02988415, which is correct to nine 
places of decimals, 


7.14. Interpolation with Two Independent Variables. ^ 
function of two independent variables ж and y ін usually 
represented by flr; у) or wey. Let the values of x be zy, ту, 





= Loe, cit, p. 188. This equation js due to W. B. Davis, Educational Timen, 
3867, p. 108. Та bae been solved by Stirling's formale by Whittaker 
Robinson, Calculus of Obeerrations, p. 01. 
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755 und those of у be уу. Va, Yar---» then the function table is 





put thus: 5 
= 9ı D э» 
= Пару) fi val Ж; у») 
Ts Kezi Уз) Mees ys) 
= fürs i уз) Nes i у») . . 





From the above table it is seen that it i» not necessary for 
= may vary while y remains constant 
interpolations with one particdlar row or 
column in the table are interpolations with one variable. The 
Problem consista in obtaining the value of f(x; y) for particular 
values of = and y, or in obtaining ап approxiaiate охрана 
of j(= ; y) ns в function of z and y, 

When the arguments proceed by equidistant intervals, say 
© and w’, we denote by E, the operation of increasing the value 
of z by e while y remains constant, and by В, the incroment 
of y by uf while x remains constant. Thus 

Efe; Wafer; y) and Eyf(r; у) ун: ye). 
но that A,f(r ; у) = fire; y) — fr; у) 
"nd — Af(r; =]; ужы) fir; у). 
Again А.А, Це; у) = Ца + ө; y + m) — firmi y) feni у + m) еу), 

‘The general formula corresponding to Gregory-Newton 

formula is 


(1) fGrmo; ут) = ЕТЕДО: y)m (1+ AL"O + A)" für y) 
ы [irmans "7 аа, ] be ene Dares] fein 
efx; м) + (mA, ena, MfG ; y) 

* [nra ema, +" Эз 

















ftiv 








(n — 1)m —2) m(m—1) m.n(n —1) 
«(n atc SEIS era lla паана, 


D az) Де; gs... 
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"This formula is due to Lambert.* 

In most cases it will be convenient to transfer the origin to 
(©, 0) end when the interval of differencing for x and y are the 
same, to take the interval as unity. The general expression for 
f ; y) cun be put as 
(9) f(e;y)=f(0;0) + (2A, уА,)/(о; 0) 


до:‏ ووا 
û7 7а +луд,д, + ag 10:»...‏ 





‘Two variable functions are of great frequency in actuarial 
work. Annuity values of joint lives аге generally tabulated for 
Quinquennia! values of z and y, and for other values methods of 
interpolation are to be used. 


Example 1. Find an expression for f(z;w) from the fol. 





lowing table - 
y=0 1 
=0 7 3 1 
1 1 в 
2 17 . 
Differencing downwards we have Д,/(0;0)= 4, A2/(00) 2 
m weross AV f(0; 01 —4, Ag/(0:0) 2. 


Also A, A,/(0; 0) 6—3—11 7 —1 


no dies رو‎ = Té وع‎ + FED, а— у MD a 


=at syty? +e- 5ys T. 
Ezample 2. Caloulate the value of joint-life annuity /(28; 17) 
from the following table: i) 
y=15 20 25 
==20 6-004 4804 3825 
25 67029 446 y 
зо 6:075 =: 


* Beytráge, Part IIJ. See also Lagrange, Nour. Mem. de Berlin (1779) 
reprinted in Eueves, З, p. 441. 
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Here the interval of differencing i» 5. Change the origin 
(20, 16) to (0, 0) and the interval to 1. ‘The required value then 
is (C074) 
Differencing downwards we bave A. = 025,  A2/= 021 
and across A, f= —1700, Aff = 721. 
Also — A,A,f(0 ; 0) = 4846 —4'304—6-020 + 6-004 = 017. 
2. (C05 4) = + OA, + 4A, —12A2 «245,4, — 12A2)/(0 ; 0) 
7767004 + "015 — 168 — "00952 + 100408 — 08662 
7525404, 
Thus — f(28; 17) 525404. 





ТАВ. When the two ages concerned are such their «um in 
a multiple of 5, the origin and the interval may be wo chosen 
that 2+ y= * 


‘The general formula then becomes 
0) fle; 1—2)=/(0; 0) + (zA, + (1—2)A, 
tizia- (А -2A,4, + APO ; 0), 





or 

( [(e;1—2)7/(0; 1) + z(/( :0) /(0 ; 1)} 
AG 1(( 5 0) —2/(1 ; 1) +000 ; 2). 

Again if 2#+y=2, the formula becomes. 
f(z; 2— a) e zz —1)/(2; 0) —2(@—2)/(1; 1) 42 —1)( -2)/0; 2). 
The problem thus reduces to one of single variable inter- 
polation. We have taken only three terms and it bas been 
found that the three term formula gives а close approximation 
to the true resulte as the formulae involving differences of higher 
order. For a detailed discussion, see Spencer, J. I. A., 40, 

рр. 298-301. 


Ezample. Consider Example 2 of p. 204. 
Changing the origin to (0, 0), (28 ; 17) becomes (6 
=+у=1. Hence by (2) 
[(0; “4) = 4304 + -6(6-029 — 4-304) 
— 76 x 4(60:075—2 x 4846 + 8/825) 
= 4304 + 1085 — "08496 = 5:25404. 
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7.16. Interpolation with two variables by Divided Differences. 
When the intervals „are not equidistant, interpolation is’ 
effected by using divided differences. Divided differences may 
be formed with respekt to = and y separately or simultaneously. 
А led difference formed with respect to z with the r+ 
arguments тү, хо ‚ and with respect to y with the #+1 
arguments уз, уз... is of order r+ and is represented by 















Suppose the values of /(r;y) are tabulated for zz, 


Жз, Ж» and у= у, Уз... Yo 





By Newton's interpolation formula $ 5.6(2) we have 


O Hes y= fers v) + Site nte nn Mery tai) 


ny n ng). (n mf 2101-0 у) 
=(2—2,)(2—24).-(2—2,) and  Xo-1. 











Again D 


еле жа. ол; ye SY flr aye Parle Vea) 





Утаа, Уз, Vases 
where Y,e(y-wiy- wal. (y—s9,) and Yo=1. 


o ENY tR > 





e fiiy) 





© 2 
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Again 


Hr We Vs Was Mad È X fry eset „аз Yo ested 








Mr, т, басти р а arr) 
and hence from (4) we have 
(5) Rez X Mr, gy eges n 
XY ir s1, ta 


The remainder term admits of being pu 
For by §5,7(2) we bave 





an elegant. form, 


für, жу, gesto tg) = 





ê: 
where £ lies between the greatest and least of =, жү, Sayan, and 
D, denotes partial differentiation, Also 


Hz, жа. ваи Mae Vae = E DENE: e Was Mae a) 
=) pp.) рее 
ni m! ies 


where q lies betwoen the greatest and least of y, уу, Vac 
Hence from (5) we g 








X» p» (dg) Xm Defers) — 
O R= 75 Di Mew) Д 
and from (4) 
X» рее; ye 3. е Janie, 
© ne Xs DEE: X, Хр Xo DIS. 


"The values of £ or у need notbe the same in all the terms in R. 
"The remainder term R vanishes when zz, xg, ... x. and 
у= уу» Уз. --- Yp and (3) will then represent a polynomial of 
(n — I)th degree in = und of (m — I)tb degree in y. 
Writing tho first fow terms in full we have 
(«у= May уз) + 2) еа з) + (®— )(®— raf ens 
^ (y iin ia Yo) + (®— ry УЗИ, 
8) tr llc valens s 
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In a similar manner for three variables =, y, =, we have 





Te 
sa- ZX Zl +R. 


TT 





fe: 


TAT. The case for equidistant arguments is obtained from 
$7.16 by making the intervals equal. As a particular case, let 


mQer-1, and y,=s—1; r—1,2,.; s=], 2. 
Then by $ 5.204) we get 





1 fee =O, 1,2,...7;0, 1, 2,0). 





l AZi(0:0, 1. 








‘This formula is the same as 17.14(2) with the remainder term. 


Again when the arguments of esch set are the same, 
when zy seg. mae and yj ==; = pos oro Бата by $5,140) 





G) Ие S ES enmt У" Dg ру Цену e B. 


fro <o 


Further when all the arguments are equal to zero, we have 
C b wt 





(3) Rem X 3 T Y D&Dif(0: 0) В. 


7.48. Another case which is of special interest* is when 
m deponds on т and is connected by the relation 
m=n-r, 720, 1, 2,0—1. . 
. 


* O. Biermann—Voriesungen über Methematiache Naherungemethoden, 
pp 15144. 
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a) Ка; у)= 


0 «то 


Writing (1) in foll we have 


@) (к; v)=foo+ (ло + Yalow + (Xafso + X,Y fas + Yalow) 
tek Goals maio + Ха ewan tet Ys авон) +В 
where by 17.10(7) 


@) Rež. Ded рр f+ 


Y. р», 
[en En каг: 1 
Here we are to culculate fewer terms than in $7.10(5). 


In the limit when the intervals becomes infinitesimal, we get 
from (2) Taylor's series for two variables. 





7.19. Construction of Difference Tables. ‘Taking the func- 
tion table on p. 208 we first form the divided differences of the 
values in the first, second, ..., columns of the function table. 
‘Thus we get the following difference table in z: 


Жк; уз) Meriva) - Iriya) 
Wey. оз Ул) May. таз Уз) Иез, оз уз) 
(ел. а, аз У) (ж, Za, таза) Me. Жан аз Уз) E 


Next we form the divided differences of the fret, second, 
rows of the above table and get the following difference table in = 
and y. 

foo (у; Ул) foy =f ri Уз, Уз) ILLUM 
Тол. zai зу = (кз. zai yis Уз) freer moi V Ve al 









У) 


‘This second table gives the divided differences to be employed 
jn the interpolation formulae. 
1419068 
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Ezample. Find an expression for fiz; y) from the following 
table: 





o x 2 з 
2 70 164 
a 138 204 
5 406 746 
в 78 eso 1104 


Forming the difference tuble in 7, we have 


2 м то 164 
з 10 68 130 
5 2 22 a2 
в ° 2 з 
and thence the difference table in = and y 
2 м эв 2з 
з 10 29 и 
5 2 10 0 
6 o 1 o 


Hence by $ 7.16(3) we have 
Sle; y) 14+ 10(2 —2) + 20 — 2) (2—8) 
+ y{28 + 20(7 — 2) + 1042 —2)(z—8) + (z — 2)(a—8)—6)) 
+ y(y—2)32 + 11 —2)) 
тази det + Ney? Dany. 


120. When the value ot аг y) №. wanted st a given point 
(кун) we need not form the second difference table in z and y. 
From the difference teble in = we calculate f(z; y) for the given 
value of x and then interpolate for tbe given value of y. Thus 
[шше fue ме ыйа ишш example we 
construct tbe following table. 4 i 
y fe gus 
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Hence #4; 5) =38 + 104y + 44y(y —2) 


= 38 + 520+ 660 = 1218, 
7.21. When the arguments sre equidistant, instead of 
forming divided difference tables, ordinary difference table in ж 


may first be formed and then ordinary difference table in = and 
У, and the result is obtained by the application of $7.14(2) or 
$770). = 


Examples 


1. 
3. 


Find the reciprocal differences of 23, (1+22)-!, 2*. 
Form the reciprocal differences of а” and shew tbat 
Рана н, zn 1i. ane n) (7 1)*a*, 





тая, 





| 1 figa) 
3. Prove that ру{лулөг,) = 





+ (уге). 
| Meoxs) Казаа) 
4. Form reciprocal differences from the following data 
= 1568, . 1560, 17570 
He) 867:61100, 556760008, — 125570550, 


and hence determine the value of /(1:5085). 


5. Interpolate the values of /(7) and /(12) by "Thieie's formula 
buving given 


z 5 10 15 20 
ta) 7835 122-05 17040 21585 


6. Express zm in the form of a continued fraction. 
7. 


Form the reciprocal derivatives of (1 + =)" and shew that 


(2— m)(8— m).. 
тн +) 









Hence develop (1 + z)" in the form of a continued fraction. 





ET CALCULUS OF FINITE DIFFERENCES 


B. Prove that 
roa- log ze п", 





ran log zlog 2420 43.1) 
Hence deduce that 


log 1+2) = = 





ives of tan 2 and shew that 





9. Obtain the reciprocal deri 
тан tan zetan т 
Temes бап z= (2n +1)(n + costz) 
a tan к= —соб x 
tan а= (n +1)(2n +1 +2 sin?z). 














D 
Hence deduce that 
tan (2+h)=ton z+ Pen = + соё z) A h 
рч” 
^ 
Ie 8 
Prove also that 
AAT h h 
aes CAO а=: d. BIER 
h h h h 


1+ Booz 1— Таваф" 


10. Determine by Aitken's process the value of log I(4) from 
Abe following data : 


E 4556 5» 718482 
№ -55038 ^ “18165 
* "42796 b "08828 
$5 "32788 18 “05261 , 
11. Find by iterated metbod tbe value of /(12) from the 
following dita: . 
5 10 15 20 25 зо 


1727628 1'62880 207893 265830 828835 4"82194 
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12, Find from the following table of annuity values 
Age 37 эв E 40 
Annuity 1540 1816 1298 1268 
‘the age corresponding to an annuity of 13°00. 
18. The equation 23 — 62 — 11-0 has a root between З and 4. 
Obtain it by inverse interpolation correct to three places of 


decimals. : 
14. Obtain by inverse interpolation the real root of the 
equations E 
[7] a3 422-2020, 
z*-2:-5 =0. 
15, From the following values of log = 
= 1 2 3 4 
log = о 780108 47712 60206 


find z when log x= "30500. 
10. From the following values of /(z) 
= n 17 a 23 з 
Ja) ‘14646 — 085267 194486 279840 923520 
find x when /(z)= 130326. , 
17. Cnleulate the valuo of joint life annuity /(91;27) from 
Аве following table: * 


y=20 25 80 
2-25 582 -555 7586, 
зо 7460 E 
35 “413 


18. If /(0; 1)= (1; 0) and /(0; 2)=/(2; 0), sbew that. 


Ha; y) (0:0) (к + {Sa + M y—VA3}/(0; 0) 
+=; 1)—f(:0)). 


19. Obtain lLagrange's formula of interpolation for two 
variables 
Ву — а) 
fo Кемо) e ge сафа) 
(z— о)(у 5) (z — а) ) 
ке eet on +18: фа 
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Hence find the value of }(31; 43) from the data: 


f(80; 40) = 14:45, /(35; 40) = 14°05, f(80; 45) = 1565,85; 45)=19-31. 
20, Find the annual 


Premium Ру, for a joint life 
insurance from the following data: 
в.а в.а 


Раз: аб =2 18 10 





Pasiss=3 4 2 
Pas 








9-811 6 


Pyonom3 14 4 


7) from the following date 
зо 35 





40 

=«=30 18-080. 19:401 19807 
35 19111 1**581 
40 12:126. 


22. Given the following annuali premiums for endowment 
assurances for specified term of years : 


Age Term 15 years — 20 years 


25 years 
?0 5'947 4418 "547 
25 0046 4:580 

зо 0144 


Obtain the annual premium for вре 25 for a term of 
17 years. 


98. Calculate the value of joint life annuity ау, from the 
following data ; 
40:50 710-804 аць 996. 
9450710591 — 445,,,20:583 
50710050 





940549 78558 
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CHAPTER І 





1 d ) о ө _ 

© ое zê сов EFI ^ () дозе sin Eri ` 
i a . 
ih Бог. im ha) 








D e a ab-o8 
(хой) аап”! аай азу Fra э Зәй + a? = Jab] + әй + B* cab «bt 





49. 3(-2)*-20—5. 
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CHAPTER II 


+. (0 22-10 @г+1) (де +з) (0245) Qr 7) 21 
10 2 


NEM 
Mic ИИ az+2 
09 B ag 5 78631)(6 3) (x 








be aa G2? + 212+17 
te) D- RFT Er FEF 


п_ 4+1 ж(ж +1) (23-2) (8x +18) 
©) ada ges: C Е 








1222 + 382 +19 
i 
= Tiar 4)8r 7) ° 


ie) 21 223 + 10211 
768 ~ 2(@> + 2)(2> + r + буе +8) ` 
e àli- (8= + 20082 +5) 
(ай Apply § 2.7 zii) Apply $ 2.6 (8). 
2. () 2*2? —4z 4 6) 6. @ бе” siete 








8 2z-l gras ут, (-in* 
в) al Axis |- Рр 





PETI 1 we EMO E nd 
* Ms Ti Руа: 2 


1 1 RE 
W i- ahr p 9 VE EF VE 





(v) 23-14427 + 1624 23); 28. 


4 © ten ,کے‎ 
а Sat? 


(6j) tan”? (342 4-18) — tan- (13) ; tan” yiy. 





в Et 4 
5. 9 (1-32. узше a 
еа 
3° (tan 20) e 
Boo 26 1 Вей. 
~ (win 26)* û” («im 26)* 
в. „ 
a? -1 
* 
п. OSE . vin ze. (i) SE cos 26. 
Aton, 
(ii) eot g eot 27—16. (i tea (e D > ten 
1 -i*i os ai] = 
® ssi are ' Vir epar 





ANSWERS 


Git) ton? (74z- 81) — tn- (81); tan” s. 
(iv) tan! (262 5) — tan-! 5; башт) $. 














-1)'** 1)6 tane 
АЕ 


ar 





1 
sin û 
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к fod eat 
d 9 coo забыв BF mer eae |. 
“л үз з 
e Gas 2) ~ sin *8 





в. (0 ве+ МИ, 12z(2—1)2—2) ‚ би 1) -2)7—8) 
2i . 81 a z 





1-1. (i) 01227 + 80r? + 282 +0 
@ 1-с». Mo TEE Iles Ress) 
r —ge-a (2-1)! 
ie) 0—2 @z—1) 
| 9. (start (i) ga ини, 


18. — 1-3(n*1)(n«2)! 


CHAPTER Ш 





арар ккк hem] 
С = 10787. 


уга T 2 2. 
baio мы TEB *812313.415 ee 








a*—0:3.17246. 2. —1%-8:%-бт+1»% 8. 9-х. 
(r—1)(2—3)2—8)-?-8z*11); 47. ^ 8. 29-241, 
27714, 2:758, 2:801, 2:843, 2:894, 2:024, 2:903, 3001, 

s 3087, 37074, 3:100, 3-143, 3-175. 





one 





of 











Ae WO шо 


1625. 
Moon's latitude for May 1170 should read. 1910710750. 
f(3)— —74,f(4)— —12. — 13. 9211220. 

"5479. 3 
£9—10.—84., 2-11 
2-17-9, 





1, 2-18—8, 2-14-8, 2—16-2, 
-19-5,3-1-2,9-2-11. 











2 871, 2-76, 27603, 27570, 2484, 2404. ж. оош. 
gaa, 25. (i) 14258... @й 5629... 
(0 07858900, (i 05235990. 27. 0197451. 
46801500. 30. (î) (000833319. (Ü) 000490080326. 
00102512, 00084056. 32. 015292156, 0:04870010, 
00047242. 35. 36:04. зв. WM. 
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